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INTRODUCTION 


This  report  describes  work  performed  by  Bushnell,  Pearsall  and 
Trozzo,  Inc.,  under  subcontract  with  the  Institute  for  Defense  Analyses 
on  Task  A-1  of  IDA  Contract  No,  £MW-C-  0749  with  the  Federal  Emergency 
Management  Agency,  Task  A-1  calls  for  the  development  of  "a  model 
to  simulate  population  movement  during  an  evacuation  from  the  risk 
area  to  the  various  host  areas  over  a  transportation  network." 

This  report  describes,  documents  and  provides  a  user's  guide  to  a 
system  of  computer  routines  which  perform  the  various  computations  re¬ 
quired  to  apply  a  crisis  relocation  model  developed  jointly  by  IDA  and 
BPT,  Inc,  The  computer  routines  together  comprise  an  interactive  system 
resident  on  the  FEMA  Univac  1108  facility.  The  model  and  its  attached 
national  data  base  can  be  used  to  analyze  in  detail  the  evacuation  of 
risk  areas  anywhere  in  the  continental  United  States  under  a  wide  range 
of  diff‘”"ent  assumptions  regarding  the  assignment  of  reception  areas  and 
the  performance  of  the  transportation  system  during  the  evacuation. 

Section  1  of  the  report  Is  a  detailed  description  of  the  IDA/BPT 
crisis  relocation  model.  The  four  basic  elements  of  the  model  are:  1) 
a  transportation  network  data  base  which  can  be  accessed  up  to  10  states 
at  a  time  to  construct  regional  highway  networks  at  the  county  level,  2) 
a  crisis  relocation  planning  submodel  to  compute  (or  accept  as  input)  an 
assignment  of  evacuees  to  reception  areas,  3)  a  route  selection  and  load¬ 
ing  algorithm  which  generates  routes  and  computes  traffic  to  approximate 
the  loadings  which  vould  result  if  evacuees  were  free  to  chose  their 
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own  routes,  and  4)  an  evacuation  simulator,  which  determines  the  location 
of  evacuees  on  the  network  as  a  crisis  relocation  proceeds.  Section  2 
is  a  user's  guide  to  the  programs  and  Section  3  is  an  annotated  sample 
run  illustrating  program  usage,  input  and  output.  Supporting  material 
is  contained  in  five  appendices. 


1.0  Description  of  the  Model 


The  IDA/BPT  crisis  relocation  model  consists  of  four  major  components 
each  with  its  own  distinct  function  and  its  own  set  of  computer  routines. 
These  components  are: 

1.  A  transportation  network  in  the  form  of  a  national  data 
base  describing  the  U.S.  highway  system  and  providing 
information  at  the  county  level  on  population,  area,  and 
risk.  The  computer  routines  allow  users  to  assemble, 
store  and  retrieve  regional  networks  comprising  up  to  10 
states. 

2.  A  crisis  relocation  planning  submodel  enabling  users  to 
directly  input  an  evacuation  plan  or  to  compute  a  plan 
as  the  solution  to  a  transportation  problem  constructed 
at  the  user's  direction. 

3.  A  route  selection  and  loading  algorithm  which  generates 
routes  and  assigns  traffic  to  the  routes  in  a  way  which 
approximates  the  loadings  which  would  actually  result  if 
evacuees  were  free  to  choose  their  own  routes  and  received 
accurate  traffic  information  during  a  crisis  relocation. 

4.  An  evacuation  simulator  designed  to  show  the  location  of 
evacuees  and  the  traffic  distribution  on  the  highway  net¬ 
work  at  various  elapsed  times  from  the  start  of  an 
evacuation. 
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1.1  The  Transportation  Network 

The  BPT/IDA  model  incorporates  a  national  data  base  describing  the 
U.S.  highway  system  and  providing  other  information  at  the  county  level. 
The  data  base  was  originally  assembled  by  IDA  in  1981  under  a  FEMA 
contract.  Further  information  on  county  areas,  1980  population  and 
vulnerability  to  attack  was  added  by  BPT  under  the  subcontract. 

In  its  present  form  the  data  base  consists  of  three  files  holding 
the  following: 

Node  File  -  Information  for  a  set  of  approximately  3700  locations 
coinciding  with  population  centroids  and  interchanges  on  the 
highway  system.  At  least  one  node  is  located  within  each  county 
in  the  continental  U.S. 

Link  File  -  Information  for  a  set  of  approximately  8700  two- 
way  segments  of  highway  joining  nodes  within  the  continental 
U.S. 


Map  File  -  A  short  file  indicating  the  location  of  Che  infor¬ 
mation  for  each  state  within  the  Node  and  Link  files. 

Annotated  partial  dumps  of  the  Node  File  and  the  Link  File  are  provided 

in  Approdix  4.5.  A  complete  listing  of  the  Map  File  is  also  given  in 

the  appendix. 

The  Node  File  contains  the  following  items  for  each  location  in  the 
data  base: 


•  the  location's  5-dlgit  state/county  FIPS  code.  The  first 
2  digits  of  the  code  designate  the  state;  the  last  3 
digits  designate  the  county  within  the  state. 

•  a  single-digit  sequence  number  for  the  location  within  the 
county.  A  sequence  number  zero  (0)  indicates  the  popu¬ 
lation  centroid  for  the  county.  Every  county  in  the 
continental  U.S.  has  a  single  population  centroid  and  may 
have  additional  nodes  at  important  highway  intersection 
points. 
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•  northing  and  easting;  the  number  of  miles  north  and  east 
from  a  reference  point  set  beyond  the  southwest  boundary 
of  each  state. 

•  a  flag  indicating  whether  or  not  the  location  is  a  major 
urban  center  (1  =  urban,  0  =  nonurban) . 

•  place  name;  the  name  of  the  closest  city,  town  or  village. 

•  county  name . 

•  latitude  and  longitude  in  degrees.  Additionally,  for 
nodes  that  are  county  population  centroids,  BPT  has  added 
the  following  items: 

-  population  of  the  county  according  to  the  1980  Census. 

-  area  of  the  county  in  square  miles. 

-  a  2-dlgit  code  corresponding  to  FEMA's  classification 
of  the  county  for  crisis  relocation  planning  purposes. 

The  first  digit  classifies  counties  according  to  risk 
in  a  countervalue  attack  as  follows: 

0  -  county  not  at  risk 

1  -  (A)  county  evacuation  plan  exists 

2  -  (B)  major  metropolitan  center;  high  risk 

3  -  (C)  intermediate-sized  city;  moderate  risk 

4  -  (D)  small  city  or  military  installation;  lower  risk 

The  second  digit  is  a  flag  designating  counties  which 
have  significant  counterforce  targets  within  their 
boundaries . 

The  Link  File  contains  the  following  items  for  each  two-way  segment 
of  the  transportation  network: 

•  the  5-digit  FIPS  code  and  sequence  number  for  each  end 
of  the  segment. 

•  the  northing  and  easting  of  the  two  ends  of  the  segment. 

•  a  2-digit  code  indicating  the  characteristics  of  the 
highway  segment.  The  first  digit  designates  the  highway 
as  1  -  Interstate,  2  -  a  U.S.  Route,  or  3  -  a  State  Road. 

The  second  digit  classifies  the  highway  segment  as 
follows: 


4 


1  -  a  6-lane  Interstate 

2  -  a  4- lane  Interstate 

3  -  a  4-lane  limited  access  highway 

4  -  a  4-lane  unlimited  access  highway 

5  -  a  2- lane  primary  road 

6  -  a  2-lane  secondary  road 

7  -  a  2-lane  poor  quality  road 

•  the  Interstate,  U.S.  Route  or  State  Road  number. 

For  purposes  of  analyzing  transportation  problems  during  a  crisis 
relocation,  it  is  rarely  necessary  to  perform  the  computations  in  a  way 
which  treats  the  continental  United  States  as  a  single  inte-rated  unit. 
The  movementsof,  say,  Michigan’s  urban  residents  along  thi  state's 
highways  have  little  practical  impact  upon  similar  movemer  in  New  York 
or  California.  Consequently,  an  evacuation  of  Michigan's  cities 

can  be  effectively  studied  within  the  context  of  a  regional  network  en¬ 
compassing  Michigan  and,  possibly,  several  neighboring  states.  The  IDA/ 
BPT  model  is  designed  to  economize  on  computer  storage  and  processor 
time  by  utilizing  only  those  parts  of  the  national  transportation  data 
base  which  are  relevant  for  a  particular  application. 

The  model  operates  on  regional  networks  of  up  to  10  states  which  are 
assembled  state  by  state  from  the  information  contained  in  the  national 
data  base.  The  states  that  are  included  in  a  single  region  need  not  be 
contiguous  although  one  would  normally  expect  them  to  be  so. 

In  order  to  facilitate  later  computations,  a  certain  amount  of  pre¬ 
processing  is  performed  on  the  node  and  link  Information  culled  from  the 
national  data  base  as  a  regional  network  is  assembled.  Nodes  are  re¬ 
numbered  sequentially  and  link  endpoints  are  equated  to  the  renumbered 
nodes.  Each  bi-directional  highway  segment  of  the  data  base  is  also 
recast  as  a  pair  of  uni-directional  links.  These  links  are  then  sorted 
by  origin  so  that  the  data  for  all  of  the  links  with  a  common  origin 
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(or  destination,  since  the  network  is  symmetric)  can  be  retrieved 
sequentially. 

Finally,  the  distance,  speed  and  capacity  of  each  link  is  computed 
using  the  information  from  the  data  base  and  parameter  values  supplied 
by  the  model  user: 

Distance  -  Link  distances  are  computed  as  the  straight  line 
distance  between  the  end  points  of  each  link  using  the 
northings  and  eastings  from  the  Link  File. 

Speed  -  Speeds  are  assigned  by  road  categories  as  follows  for 
normal  (nonevacuation)  movements  along  the  network. 

Class  1-3,  Interstates  and  other  4-lane  limited  access 
roads  -  55  mph. 

Class  4,  Primary  unlimited  access  4-lane  -  45  mph. 

Class  5,  Primary  unlimited  access  2  lane  -  40  mph. 

Class  6-7,  Secondary  and  poor  quality  roads  -  30  mph. 

Speeds  for  crisis  relocation  (evacuation)  movements  are  sup¬ 
plied  by  users  for  the  same  road  categories. 

Capacity  -  For  normal  movements  the  network  is  regarded  as  un¬ 
capacitated;  for  crisis  relocation  the  user  supplies  parameter 
values  for  the  vehicles  per  lane  per  hour  using  the  same  road 
categories  as  speed.  These  values  are  multiplied  by  the  num¬ 
ber  of  lanes  to  obtain  the  capacity  in  each  direction. 

The  number  of  lanes  available  in  each  direction  is  partly  a  function 

of  traffic  control  during  an  evacuation.  Users  of  the  IDA/BPT  model  may 

choose  among  three  options  for  controlling  traffic  during  an  evacuation. 

Normal  Traffic  Control  -  Two-way  traffic  on  all  highway  seg¬ 
ments.  The  number  of  lanes  in  each  direction  is  one-half  the 
number  of  lanes  given  in  the  dat j  base. 

One-Way  Outbound  Traffic  Control  -  One-way  traffic  on  all  con¬ 
gested  highway  segments.  All  lanes  of  these  segments  art 
assigned  to  the  outbound  link.  Inbound  links  for  these  high¬ 
way  segments  are  assigned  zero  capacity. 

Variable  Lanes  Traffic  Control  -  All  lanes  are  available  to 
traffic  in  both  directions.  In  effect,  the  model  combines 
the  traffic  and  available  capacity  in  both  directions  along 
a  highway  segment. 
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Since  the  assembly  of  a  regional  network  from  the  Node,  Link  and 
Map  files  is  a  somewhat  tedious  job,  the  computer  program  which  imple¬ 
ments  the  model  provides  for  the  storage  and  retrieval  in  preprocessed 
form  of  up  to  10  regional  networks. 
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1.2  The  Crisis  Relocation  Planning  Submodel 

Within  the  IDA/BPT  evacuation  model,  a  crisis  relocation  plan  is 
functionally  defined  as  a  set  evacuation/reception  county  pairs  and  the 
numbers  of  evacuees  to  be  relocated  from  the  evacuation  counties  to  the 
reception  counties.  Users  of  the  model  have  the  option  of  specifying 
the  evacuation  plan  directly  or  supplying  other  information  from  which 
a  linear  program  is  constructed  and  solved  to  obtain  an  evacuation  plan. 
This  section  describes  the  computation  of  an  evacuation  plan  which  con¬ 
forms  to  criteria  specified  by  the  model  user. 

The  linear  program  which  is  solved  to  obtain  an  evacuation  plan  is 
known,  technically,  as  a  transportation  problem.  Formally,  the  problem 
has  the  following  mathematical  statement:  Let  i  =  1,  M  designate  origin 
counties,  j  =  1,  N  designate  destination  counties,  and  denote  the 
number  of  evacuees  to  be  relocated  from  origin  i  to  destination  j.  Let 


S^,  i  =  1,  M  be  the  numbers  of  evacuees  to  be  relocated  from  the  evacua¬ 
tion  counties,  let  D j ,  j  =1,  N  be  the  capacities  in  numbers  of  evacuees 


that  can  be  accommodated  in  the  reception  counties  and  let 


be  the  time 


it  takes  to  travel  over  the  transportation  network  from  evacuation  county 


i  to  reception  county  j.  Then,  the  linear  program  which  is  solved  to 


obtain  an  evacuation  plan  is: 
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An  efficient  algorithm  for  solving  transportation  problems  has  been 
incorporated  into  the  computer  routines  which  implement  the  IDA/BPT 
transportation  model.  This  algorithm  is  applied  to  solve  transportation 
problems  which  are  constructed  employing  information  supplied  by  users. 
The  nonzero  components  of  the  solution  to  the  transportation  problem  con¬ 
stitute  an  evacuation  plan. 

The  elements  of  the  transportation  problem  which  are  specified 
directly  or  indirectly  by  users  of  the  IDA/BPT  model  are  the  evacuation 
counties  and  the  number  of  evacuees,  Sj^,  i  =  1,  M,  and,  the  reception 
counties  and  the  reception  capacities,  Dj,  j  =  1,  N.  The  times,  are 

computed  as  the  normal  travel  times  over  the  regional  highway  network 
from  the  population  centroids  of  the  evacuation  counties  to  the  popula¬ 
tion  centroids  of  the  reception  counties.  The  times  are  determined  for 
normal  speeds  over  the  shortest  route.  The  rationale  for  inserting 
normal  travel  times  into  the  transportation  problem  is,  first,  that  the 
generally  slower  speeds  of  travel  during  an  evacuation  have  little  effect 
on  the  choice  of  the  shortest  routes  and,  second,  a  great  deal  of  traffic 
at  normal  speeds  can  be  expected  following  the  relocation  between  an 
evacuation  county  and  its  reception  counties.  This  traffic  is  necessary 
to  support  the  relocated  population  and,  possibly,  to  maintain  production 
and  services  within  the  evacuated  county. 

A  list  of  evacuation  counties  together  with  the  number  of  evacuees 
from  each  county  can  be  specified  directly  by  the  user.  Alternatively, 
evacuation  counties  in  a  region  can  be  selected  according  to  risk  cate¬ 
gories  given  in  FEMA's  Crisis  Relocation  Conglomerate  Listing.  These 


categories  are: 
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A  -  counties  for  which  crisis  relocation  plans  already  exist; 
usually  because  a  major  military  installation  is  located 
within  the  county. 

B  -  the  most  densely  populated  counties;  counties  containing  major 
metropolitan  areas  which  could  be  targeted  in  the  event  of 
an  attack. 

C  -  heavily  populated  counties;  usually  the  counties  contain 

a  major  city  with  a  population  of  several  hundred  thousand. 

D  -  less  densely  populated  counties;  counties  containing 
either  a  concentration  of  population  and  industry  or  a 
military  installation  which  could  make  the  county  a 
target  in  the  event  of  an  attack. 

The  IDA/BPT  model  allows  users  to  specify  evacuation  counties  as  1)  all 
A  counties,  2)  all  B  counties,  3)  all  B  and  C  counties,  4)  all  A,  B,  C 
and  D  counties.  In  addition,  FEMA's  Crisis  Relocation  Conglomerate 
Listing  designates  those  counties  which  would  be  threatened  in  the  event 
of  a  counterforce  attack.  These  counterforce  targets  can  be  added  to 
the  list  of  evacuation  counties  by  users  of  the  model. 

The  number  of  evacuees  from  each  evacuation  county  selected  in  this 
way  is  computed  using  a  set  of  rules  as  follows: 

1.  If  counterforce  targets  are  included  in  the  list  of 
evacuated  counties,  then  all  counties  with  counterforce 
targets  are  completely  evacuated.  The  number  of 
evacuees  is  set  equal  to  the  county's  1980  population. 

2.  All  other  counties  are  evacuated  down  to  the  point  where 
the  remaining  population  meets  both  of  the  following 
criteria: 

a.  The  population  following  evacuation  does  not  exceed 
a  percentage  of  the  county's  residents  specified  by 
the  user,  e.g.,  10  percent. 

b.  The  population  density  does  not  exceed  a  number  of 
people  per  square  mile  specified  by  the  user,  e.g., 

200  people  per  square  mile. 

The  designation  of  reception  counties  is  similar  to  the  designation 


of  evacuation  counties.  Again,  users  have  the  option  of  directly 
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specifying  these  counties  and  their  capacities  to  hold  evacuees. 
Alternatively,  the  following  rules  can  be  applied  to  establish  a  list 
of  reception  counties; 

1.  No  evacuation  county  is  simultaneously  allowed  to  serve 
as  a  reception  county. 

2.  No  county  is  permitted  to  be  a  reception  county  if  its 
resident  population  density  is  less  than  a  specified 
minimum  value  such  as  10  people/square  mile.  The  pur¬ 
pose  of  thf  minimum  density  Is  to  exclude  as  reception 
counties  areas  such  as  mountains  and  deserts  which  cannot 
support  significant  numbers  of  evacuees. 

3.  No  county  is  permitted  to  be  a  reception  county  if  its 
resident  population  exceeds  a  specified  maximum  value 
such  as  500  people/square  mile.  The  maximum  density 
limit  enables  the  model  user  to  exclude  counties  which 
are  already  heavily  populated  bv .  have  not  been  designated 
as  evacuation  counties. 

All  counties  within  a  region  which  are  not  excluded  by  1-3  above 
are  included  as  potential  reception  counties.  The  number  of  evacuees 
each  potential  reception  county  can  hold  is  computed  according  to  the 
following  criteria: 

1.  The  population  of  a  reception  county  following  an  evacuation 
cannot  exceed  a  user-specified  multiple  of  the  county's 
original  population.  For  example,  a  multiple  of  three 
would  limit  the  capacity  of  leception  counties  to  no  more 
than  two  evacuees  for  each  resident.  By  limiting  recep¬ 
tion  capacity  in  this  way  users  can  prevent  the  model  from 
assigning  more  evacuees  to  a  reception  county  than  the 
local  infrastructure  can  support. 

2.  The  population  density  following  the  evacuation  does  not 
exceed  a  number  of  people  r  square  mile  stipulated  by 
the  user,  i.e.,  500  people/square  mile.  This  maximum 
limit  lets  users  prevent  the  creation  of  potential  counter- 
value  targets  In  reception  counties. 

The  model  also  allows  users  to  establish  a  ger graphic  band  of  pre¬ 
ferred  reception  counties  around  each  evacuation  county.  The  inner  and 
outer  limits  of  these  bands  are  stipulated  by  the  model  user  in  the  form 
of  a  minimum  and  a  maximum  normal  travel  time.  For  example,  limits  of 
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one-half  hour  and  two  hours  would  have  the  effect  of  establishing  a  band 
of  preferred  reception  counties  within  commuting  distance  of  most  popu¬ 
lation  centers. 

The  travel  time  limits  are  used  to  modify  the  coefficients,  t^^,  of 
the  transportation  problem.  If  the  normal  travel  time  between  an  evacua¬ 
tion  county  and  a  reception  county  falls  outside  the  limits,  then  1,000 
hours  is  arbitrarily  added  to  the  normal  shortest  route  time  between  the 
two  counties.  The  effect  of  this  change  is  to  make  the  model  select 
evacuation  plans  in  which  the  reception  capacity  within  the  preferred 
counties  is  entirely  used  up  before  any  evacuees  are  assigned  to  reception 
counties  outside  the  travel  time  limits.  In  the  example  given,  the  solu¬ 
tion  to  the  transportation  problem  would  strongly  favor  assignments  of 
evacuees  to  reception  counties  between  one-half  hour  and  two  hours  normal 
travel  time  from  their  homes. 
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1.3  The  Route  Selection  and  Loading  Algorithm 

The  heart  of  the  IDA/BPT  crisis  relocation  model  is  an  efficient 
algorithm  for  finding  routes  along  the  regional  transportation  network 
and  then  determining  the  traffic  the  routes  would  carry  during  an  actual 
evacuation.  The  end  product  of  the  route  selection  and  loading  algorithm 
is  a  list  of  routes,  no-less-than  one  for  each  evacuation/reception 
county  pair  in  the  crisis  relocation  plan,  with  an  assignment  of  traffic 
to  the  routes  such  that  for  any  given  evacuation/reception  pair; 

1.  all  evacuees  from  the  evacuation  county  to  the  reception 
county  are  assigned  to  one  of  the  routes  connecting  the 
population  centroids  of  the  two  counties . 

2.  all  routes  assigned  traffic  are  equivalent  in  terms  of 
average  total  transit  time.  Total  transit  time  consists 
of  actual  travel  time  from  origin  to  destination  plus 
the  length  of  the  average  delay  imposed  by  the  most  con¬ 
gested  link  along  the  route. 

3.  all  other  possible  routes  have  longer  average  total 
transit  times  than  those  which  are  assigned  traffic  by 
the  model.  (NOTE:  See  comments  on  setting  the  delay 
time  termination  limit.) 

In  brief,  the  algorithm  contained  in  the  IDA/BPT  model  selects  routes  and 
loads  them  in  a  way  which  approximates  the  loading  which  would  result 
if  evacuees  were  aware  of  the  travel  and  delav  times  occurring  on  tho 
network  and  were  free  to  choose  their  »'wn  routes. 

Normally,  the  I'.S.  highwav  svsteir  Is  un.  on  ge  s  t  t-vl  exiert  ^or  certain 
links  at  well-known  iiours,  i.e.,  the  '.ei'rg*  c  Htiage  ^t  r  ish 

hour  and  U.S.  Route  >0  between  >,a>liiiigt  <■.  it>.'  'i  ■.  i;  1 1  t  *•  r 

noons  during  the  summer.  Hi’wev'M  .  <  :  ■.  t-  •  ■ 

would  place  such  e  xt  r.iot  1 1  n.ii  .  i.-m.i;.  •  .  -  ■■ 
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congestion  of  major  outbound  highways  would  become  the  rule  rather  than 
the  exception.  This  fact  has  two  Important  Implications  for  evacuation 
studies  and  plans.  First,  the  technical  characteristics  associated  with 
highway  usage  such  as  average  speed  and  lane  capacity  can  be  expected  to 
degrade  and  should  be  somewhat  influenced  by  traffic  control  during  an 
evacuation.  Second,  roundabout  routes  which  would  ordinarily  be  of 
little  interest  become  important  as  alternatives  to  the  heavily  congested 
shortest  routes . 

The  IDA/BPT  model  allows  users  to  stipulate  the  major  performance 
parameters  of  the  highway  system  during  an  evacuation.  These  performance 
parameters  are: 

1.  Vehicles/ lane/hour  for  four  highway  categories:  Class 
1-3:  Interstates  and  other  4-lane  limited  access  roads. 

Class  4:  Primary  unlimited  access  4-lane  highways. 

Class  5:  Primary  unlimited  access  2-lane  highways. 

Class  6-7:  All  other  roads. 

2.  Average  speeds  of  movement  in  miles/hour  for  the  same 
four  highway  categories. 

3.  Traffic  control;  three  options  are  available:  normal 
two-way  traffic,  one-way  outbound  traffic  on  congested 
links  and,  variable  lanes  control  in  which  the  traffic 
and  capacity  in  both  directions  along  a  segment  is  pooled. 

4.  Average  number  of  passengers  per  vehicle.  The  number  of 
passengers  per  vehicle  is  used  to  convert  evacuees  to 
vehicles  along  routes.  Normally,  the  number  of  occu¬ 
pants  per  vehicle  is  low,  between  one  and  two.  However, 
enforced  car  pooling,  intensive  busing  and  the  fact  that 
families  would  typically  evacuate  as  a  unit  should  raise 
the  average  during  an  evacuation. 

5 .  Normal  evacuation  time  in  hours .  It  would  take  time  to 
effect  an  evacuation  even  if  the  highway  system  imposed 
no  additional  delays.  Evacuees  must  be  informed  that  an 
evacuation  is  underway,  they  must  assemble  with  their 
families,  do  a  necessarily  minimal  mount  of  preparation 
and  then  begin  the  evacuation  by  travelling  on  local 
roads  to  the  major  arteries  of  the  highway  system.  The 
time  it  would  take  for  the  last  evacuees  to  reach  the 
highway  system  and  begin  moving  outbound  is  the  nominal 
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evacuation  time .  Half  the  nominal  evacuation  time  is 
the  average  delay  encountered  by  evacuees  just  reaching 
the  regional  transportation  network. 

Formally,  a  route  is  defined  as  a  sequence  of  links  v.iich  begins  at 

the  population  centroid  of  an  evacuation  county  and  ends  at  the  centroid 

of  a  reception  county.  The  total  transit  time  for  a  route  is  composed 

of  two  parts  as  follows: 

1.  Average  travel  time:  the  total  time  required  to  travel 
the  route  computed  as  the  sum  of  the  times  required  to 
cross  each  link  at  the  speed  established  for  the  links  by 
the  model  user. 

2.  Average  delay  time:  one-half  the  "'essor  of  the  nominal 
evacuation  time  or  the  length  of  the  delay  imposed  by  the 
most  congested  link  on  the  route.  The  delay  Imposed  by 

a  link  is  computed  from  the  formula: 

Average  delay  _  number  of  evacuees  traversing  the  link/passengers  per  vehicle 
in  hours  vehicles  per  lane  per  hour  x  number  of  available  lanes 

The  number  of  evacuees  traversing  a  link  is  the  sum  of 
all  evacuees  assigned  to  routes  which  include  the  link 
(in  both  directions  for  variable  lanes). 

The  route  selection  and  loading  algorithm  is  an  iterative  process 

in  which  route  generation  and  network  loading  are  alternated  until  no  new 

routes  can  be  found  to  improve  the  solution. 

Step  1  -  Initialization;  all  links  on  the  regional  transpor¬ 
tation  network  are  unblocked  and  may  be  included  in  routes. 

The  list  of  routes  is  empty. 

Step  2  -  Route  generation;  the  shortest  route  following  an  un¬ 
blocked  sequence  of  links  is  generated  for  each  origin/ 
destination  pair  of  county  centroids  in  the  relocation  plan. 

The  method  used  to  generate  the  best  route  from  an  evacuation 
county  to  a  reception  county  is  described  below.  A  route  is 
added  to  the  route  list  if  it  does  not  duplicate  a  route 
already  on  the  list  and  its  transit  time  is  shorter  than  the 
transit  times  of  the  routes  already  on  the  list  for  the  same 
origin/destination  pair. 

Step  3  -  Route  loading;  if  no  new  routes  are  added  to  the 
route  list  in  Step  2,  the  algorithm  terminates.  Otherwise, 
loadings  for  the  new  list  of  routes  are  computed  using  an 
iterative  method  for  assigning  traffic  to  a  given  set  of 


routes.  The  iterative  method  is  described  below.  It  con¬ 
verges  fairly  rapidly  on  a  route  loading  which  equalizes 
total  travel  times  for  all  routes  connecting  the  same 
origin/destination  pairs.  (If  Step  3  is  being  executed 
for  the  first  time,  then  the  route  loading  is  computed 
simply  by  assigning  the  evacuees  for  each  origin/desti¬ 
nation  pair  to  the  single  route  on  the  route  list  for 
that  pair.) 

Step  4  -  Link  blocking;  the  network  is  scanned  for  the 
unblocked  link  with  the  longest  average  delay  time  Implied 
by  the  route  loadings  of  Step  3.  This  link  and  all  other 
links  with  comparable  delay  times  are  blocked.  In  practice 
all  links  with  delay  times  of  90  percent  or  more  of  the 
longest  average  delay  time  are  blocked.  Blocking  a  link 
removes  the  link  from  the  network  only  for  the  purpose  of 
generating  routes  in  Step  2.  Blocking  links  has  the  effect 
of  forcing  the  route  generator  to  find  link  sequences  which 
avoid  those  parts  of  the  network  which  will  cause  delays 
during  an  evacuation. 

Step  5  -  Link  list  editing;  if  the  longest  average  delay 
time  found  in  Step  4  is  less  than  a  user-specified  delay 
time  termination  limit,  the  algorithm  terminates.  Other¬ 
wise,  the  link  list  is  edited  by  deleting  unloaded  or  highly 
loaded  routes  with  transit  times  which  significantly  exceed 
the  transit  times  on  loaded  routes  between  the  same  origin/ 
destination  pair.  The  algorithm  then  returns  to  Step  2. 

Routes  are  generated  in  Step  2  by  applying  an  efficiently  programmed 

version  of  the  well-known  Ford-Fulkerson  method  for  finding  the  shortest 

route  between  two  points  along  an  uncapacitated  network.  The  algorithm 

is  more-than-adequately  described  in  many  operations  research  textbooks, 

for  example.  Principles  of  Operations  Research  by  Harvey  M.  Wagner. 

(2nd  edition,  1975,  Prentice-Hall,  Englewood  Cliffs,  N.J.) 

Basically,  there  are  two  phases  to  the  computations.  First,  the  shortest 
time  to  the  destination  node  is  determined  for  every  point  on  the  network. 
These  times  are  then  used  in  the  second  phase  of  the  algorithm  which  is 
to  determine  the  actual  shortest  route  from  any  given  origin  to  the 
destination. 

The  efficient  programming  of  the  algorithm  is  made  possible  by  the 
presorting  of  the  network  links  described  earlier.  This  device  for 
speeding  computations  was  developed  by  BPT  in  a  previous  project  for  the 
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Department  of  Transportation  of  the  State  of  Michigan.  The  presorted 
network  simplifies  the  computations  in  both  phases  of  the  Ford-Fulkerson 
method  by  allowing  the  link  list  to  be  processed  sequentially  in  segments. 

In  applying  the  Ford-Fulkerson  method  to  generate  routes  in  Step  2, 
the  usual  computational  roles  of  the  origin  and  destination  are  reversed. 
Times  are  computed  from  every  point  on  the  network  to  each  evacuation 
county  and  routes  are  generated  as  though  evacuees  were  moving  from 
reception  counties  to  evacuation  counties.  Since  the  network  is  symmetric, 
routes  generated  in  this  way  are  identical  when  the  link  sequence  is 
inverted  to  those  which  would  be  generated  if  origins  and  destinations 
were  not  reversed.  However,  the  reversal  of  the  origins  and  destinations 
speeds  computations  since  the  first  phase  of  the  Ford-Fulkerson  method 
is  performed  once  for  each  evacuation  county  rather  than  once  for  each 
reception  county.  Ordinarily,  there  are  far  few  evacuation  counties 
than  reception  counties. 

Route  loading  is  performed  in  Step  3  of  the  algorithm  by  employing 
a  simple  iterative  process.  The  process  begins  with  the  route  loading 
from  the  previous  execution  of  Step  3.  Transit  times  for  each  route  on 
the  route  list  are  computed  using  this  initial  loading.  The  route  with 
the  shortest  transit  time  for  each  origin/destination  pair  is  determined 
and  a  provisional  network  loading  is  constructed  by  assigning  all  of  the 
evacuees  for  each  origin/destination  pair  to  the  route  with  the  shortest 
transit  time.  A  new  network  loading  is  now  obtained  by  forming  a 
weighted  average  of  the  original  network  loading  and  the  provisional 
network  loading. 

The  process  is  now  repeated  with  the  new  network  loading  replacing 
the  network  loading  used  to  initiate  Step  3.  A  new  provisional  loading 


is  constructed,  a  new  weighted  average  Is  formed,  and  the  new  network 
loading  is  used  to  initiate  another  cycle  of  the  process. 

The  process  is  continued  for  a  fixed  number  of  iterations.  The 
weights  which  are  used  to  form  the  weighted  average  are  a  function  of 
the  number  of  iterations.  Let  N  be  the  number  of  iterations,  then  the 
weight  applied  to  the  provisional  loading  at  iteration  N  is  given  by 
W  =  2/(N  +  19)  and  the  weight  for  the  original  loading  is  1-W.  One 
hundred  iterations  are  performed  for  each  execution  of  Step  3.  At  the 
termination  of  the  route  selection  and  loading  algorithm  an  additional 
four  hundred  iterations  are  performed  to  refine  the  final  network  loading 

The  overall  route  selection  and  loading  algorithm  terminates  if 
either  one  of  two  conditions  exist.  These  are: 

1.  In  Step  1,  if  the  route  generator  fails  to  find  any  new 
routes  to  add  to  the  existing  route  list.  After  several 
Iterations  through  Step  4  the  regional  highway  network 
tends  to  get  reduced  and  disconnected.  At  some  point  so 
many  links  may  be  blocked  that  there  exist  no  new  routes 
meeting  the  criteria  for  inclusion  on  the  route  list  in 
Step  1. 

2.  In  Step  4,  if  the  delay  times  on  all  of  the  remaining 
unblocked  links  in  the  network  are  less  than  the  delay 
time  termination  limit.  If  the  delay  time  termination 
limit  is  set  less  than  or  equal  to  half  the  nominal  de¬ 
lay  time  then  the  algorithm  will  not  stop  until  there 
are  no  more  routes  with  longer  average  total  transit 
times  than  those  which  are  assigned  traffic  by  the  model. 

If  the  delay  time  termination  limit  is  set  higher  than 
half  the  nominal  delay  time  then  the  algorithm  will 
terminate  earlier  with  fewer  routes.  However,  since 

the  algorithm  generates  the  best  routes  first  there  is 
frequently  little  to  be  gained  by  carrying  computations 
to  the  point  where  every  possible  route  has  been  examined. 
Usually,  a  good  computational  strategy  is  to  set  the 
delay  time  termination  limit  above  half  the  nominal 
delay  time  but  at  a  level  which  is  still  small  in  compari¬ 
son  to  the  delays  likely  to  be  encountered  by  evacuees 
on  the  most  congested  links  of  the  network. 
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1.4  The  Evacuation  Simulator 

The  IDA/BPT  model  enables  users  to  perform  a  computer  simulation  of 
an  evacuation  conducted  according  to  a  previously  generated  set  of  routes 
and  route  loadings.  The  simulation  provides  a  picture  of  the  progress 
of  the  crisis  relocation  at  intervals  of  2,  4,  6,  8,  12,  16,  20,  24,  32, 

40,  48,  60,  72,  and  84  hours  after  the  start  of  the  evacuation. 

The  simulation  is  performed  route  by  route  by  determining  the  loca¬ 
tion  of  evacuees  at  each  of  the  elapsed  times  given  above.  Evacuees  are 
assumed  to  embark  on  each  route  in  a  steady  stream  commencing  at  time  zero 
and  lasting  for  twice  the  length  of  the  average  delay  time  for  the  route. 
The  number  of  evacuees  embarking  per  hour  is  the  total  number  of  evacuees 
assigned  to  the  route  by  the  route  loading  algorithm  divided  by  twice  the 
average  delay  time.  The  rate  at  which  vehicles  are  fed  onto  the  route  is 
obtained  by  dividing  the  number  of  evacuees  per  hour  by  the  average  number 
of  passengers  per  vehicle.  In  effect  the  delay  which  attaches  to  the  use 
of  a  route  is  imposed  on  evacuees  at  the  start  of  the  route  and  not  at 
the  location  of  the  most  congested  link  which  would  actually  be  responsible 
for  the  delay. 

The  stream  of  evacuees  using  a  particular  route  is  of  a  fixed  length 
in  time  and  volume  in  vehicles  per  hour.  As  this  stream  proceeds  out 
over  the  route  towards  its  destination  in  a  reception  county  the  links 
along  the  route  are  successively  entered,  covered  and  exited  by  the  stream 
of  vehicles.  The  situation  for  a  typical  time  in  mid-evacuation  is  shown 
below.  The  arrow  represents  the  stream  of  vehicles  using  the  route.  Its 
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twice  average  delay 

progress  along  the  route  is  governed  by  the  speeds  of  movement  on  the 
sequence  of  links  that  form  the  route.  When  the  stream  completely  covers 
a  link,  such  as  line  2,  the  volume  of  the  stream  is  added  to  the  traffic 
using  the  link.  Where  a  link  is  only  partly  covered,  such  as  links  1  and 
3,  a  proportion  of  the  volume  equal  to  the  fraction  of  the  link  which  is 
covered  by  the  stream,  is  added  to  traffic  on  the  link.  A  record  is  also 
kept  of  the  numbers  of  evacuees  not  yet  departed  from  the  evacuation 
county  and  already  arrived  at  the  reception  county. 

The  status  of  the  crisis  relocation  at  each  interval  of  time  is 
summarized  for  model  users  in  a  series  of  tables  which  display: 

•  The  status  of  all  evacuees  as  at  risk  (unevacuated),  en 
route  or  at  distination. 

•  The  distribution  of  departures  by  evacuation  counties. 

•  The  distribution  of  arrivals  by  reception  counties. 

•  The  distribution  of  traffic  on  the  regional  highway 
network. 


2.0  User's  Guide  to  the  Computer  Program 

The  IDA/BPT  crisis  relocation  model  has  been  programmed  in  FORTRAN 
and  made  available,  together  with  the  data  files  comprising  its  data 
base,  on  FEMA's  Univac  1108  computer. 

All  of  the  model's  source  data  and  source  programs  are  available  as 
public  files  on  the  Univac  1108.  Procedure  files  which  may  be  added  to 
the  run  stream  to  compile,  assemble  and  execute  the  source  programs  are 
listed  in  the  appendix. 

The  computer  program  is  designed  for  interactive  usage  from  a  suit¬ 
able  computer  terminal  signed  on  to  the  Univac  1108  in  the  demand  process¬ 
ing  mode.  A  suitable  terminal  should  have  a  BAUD  rate  of  300  or  greater 
and  be  capable  of  printing  single  lines  of  132  characters.  After  signon, 
users  should  set  the  line  length  to  132  characters  with  an  @@TTY  or  @@DCT 
command.  FEMA  users  may  run  the  IDA/BPT  computer  program  by  adding  three 
prepared  files  of  commands  to  the  run  stream  with  the  commands: 

(3(aADD  BPT  *  EXECUTE.  COMPILE 
i3(3ADD  BPT  *  EXECUTE.  ASSEMBLE 
@(3ADD  BPT  *  EXECUTE.  RUN 

The  commands  must  be  executed  in  the  order  given.  The  first  command 
compiles  the  ASCII  FORTRAN  source  code.  The  relocatable  object  code  is 
stored  iu  a  temporary  file  named  OBJECT.  The  second  command  assembles 
the  relocatable  object  code  in  an  element  named  OBJECT.  MAIN.  The  third 
command  links  the  source  data  files  to  the  channels  8,  9  and  10;  links 
channel  7  to  a  temporary  file  named  "ROUTES"  and  links  channels  11 
through  20  to  temporary  files  named  FILEll,  FILE12,  etc.  which  can  be 
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conveniently  referenced  by  their  channel  numbers.  The  program  is  then 
run. 

We  have  attempted  to  make  the  operation  of  the  model  as  self-evident 
and  forgiving  to  a  terminal  user  as  possible.  Thus,  users  are  always 
prompted  with  an  appropriate  question  or  request  at  all  points  where  the 
program  requires  user-supplied  information.  In  addition,  a  sufficient 
amount  of  terminal  output  is  supplied  to  enable  users  to  verify  input 
values  and  to  confirm  the  completion  of  major  computational  segments. 
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2.1  Master  Control  Sequence 

Subject  to  several  sequencing  rules,  the  order  in  which  the  major 
component?  of  the  IDA/BPT  model  are  exercised  is  under  the  control  of 
the  user.  The  computer  program  always  begins  by  reproducing  the  follow¬ 
ing  listing  of  control  options  at  the  terminal: 

0  -  Repeat  List  of  Control  Options 

1  -  Input /Compute  Crisis  Relocation  Plan 

2  -  Run  Route  Finding/ Loading  Algorithm 

3  -  Simulate  the  Evacuation 

4  -  Assemble  Network  from  Source  Files 

5  -  Read  Stored  Network 

6  -  Save  Network 

7  -  Reassign  Input  Channel  for  Parameter  Values 

8  -  Reassign  Output  Channel  for  Summary  Tables 

9  -  Reassign  Output  Channel  for  Route  Listing 
99  -  Terminate  Program 

After  execution  of  major  program  segments,  the  program  always  returns  to 
the  master  control  sequence  with  the  following  printed  request: 

INPUT  CONTROL  OPTION  NUMBER 

The  user's  response  at  the  terminal  determines  the  next  major  program 
operation  subject  to  the  following  sequencing  restrictions: 

1.  Options  1  and  6  cannot  be  selected  until  after  either 
Option  4  or  5  has  been  executed. 

2.  Option  2  cannot  be  selected  until  after  Option  1  has  been 
executed  and  a  crisis  relocation  plan  input  or  computed. 

3.  Option  3  cannot  be  selected  until  after  Option  2  has 
been  executed  and  a  route  loading  has  been  computed. 

Since  exercising  control  options  often  causes  the  loss  or  destruction 

of  internally  stored  computational  results,  the  program  requires  the 

verification  of  control  option  numbers  with  the  request: 
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CONTROL  OPTION  n  TYPE  "OK"  TO  EXECUTE 
The  response  "OK"  leads  directly  to  execution  of  the  option;  any  other 
response  recycle  the  master  control  sequence. 
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2.2  Input/Output  Channel  Options 

The  program  receives  parameter  values,  prints  summary  tables 
following  most  major  computational  segments  and  prints  detailed  list¬ 
ings  of  routes  on  completion  of  the  route  finding/loading  algorithm. 
Initially,  parameter  values  are  input  via  the  terminal  (channel  5), 
summary  tables  are  printed  at  the  terminal  (channel  6)  and  route  list¬ 
ings  are  sent  to  a  temporary  file  named  "ROUTES"  which  is  attached  to 
the  channel  7  and  can  be  referenced  by  users  as  "7." 

These  assignments  can  be  altered  by  selecting  Options  7,  8  or  9. 
Following  selection  of  each  of  these  options,  a  request  for  a  channel 
number  is  printed  at  the  terminal: 

INPUT  NEW  CHANNEL  NUMBER  (5  or  6  =  TERMINAL) 

Assignments  of  new  channel  numbers  must  obey  the  following  rules: 

1.  Parameter  values  may  be  input  through  channel  5  and 
channels  11  to  20. 

2.  Summary  tables  and  route  listings  may  be  output  through 
channel  6  and  channels  11  to  20.  The  route  listings 
may  also  be  output  through  channel  7. 

3.  Different  channels  must  be  used  for  parameter  values, 
summary  tables  and  route  listings. 

If  a  channel  other  than  6  (the  terminal)  is  designated  for  the  sum¬ 
mary  tables,  copies  of  the  very  short  tables  which  are  printed  at  the 
terminal  are  also  output  through  the  summary  table  channel.  Summary 
tables  which  are  output  through  channels  other  than  6  are  also  numbered. 
The  numbering  system  follows  the  numbering  of  the  control  options.  Fol¬ 
lowing  is  a  complete  list  of  the  tables  produced  by  the  program: 
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Table  1.0 
Table  1.1 
Table  1.2 
Table  2.0 
Table  2.1 
Table  2.2 
Table  2.3 
Table  3.0 
Table  3.1 
Table  3.2 
Table  3.3 
Table  4.0 
Table  4.1 
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Crisis  Relocation  Input  Parameter  Summary 

Evacuation/Reception  County  Lists 

Evacuation/Reception  County  Pairs 

Network  Finding/ Loading  Parameter  Summary 

Algorithm  Iteration  Record 

Route  Characteristics  and  Loadings 

Route  Listings 

Simulation  Summary 

Evacuation  Departure  Summary 

Reception  Arrival  Summary 

Network  Link  Characteristics  and  Loadings 

Regional  Network  Assembly 

Population  Centroid  List 


2.3  Network  Assembly,  Storage  and  Retrieval 

Control  Options  4,  5  and  6  provide  users  with  the  ahilitv  to 
assemble  regional  transportation  networks  from  the  source  data  files, 
save  assembled  networks  and  retrieve  previously-saved  networks. 

Under  Option  4  regional  networks  of  up  to  ID  states  can  be  assembled 
from  the  source  files  comprising  the  model's  data  base.  The  size  of 
the  regional  network  is  limited  to  no  more  than  500  nodes  and  2,000 
one-way  links.  This  is  sufficient  to  accommodate  most  regional  networks 
of  interest  and  the  limits  can  readily  be  expanded  if  necessary. 

To  assemble  a  regional  network,  the  user  nust  specify  the  number 
of  states  to  be  included  in  the  network  and  the  FIPS  numbers  for  the  states 
in  response  to  the  requests: 

INPUT  THE  NUMBER  OF  STATES  IN  THE  REGION 

and 

INPUT  n  STATE  NUMBERS 

A  directory  of  states  listed  by  FIPS  number  is  displayed  in  the  table  on 
the  following  page.  The  state's  FEMA  region  number  is  also  shown  as  an 
aid  to  users .  The  computer  program  echos  the  FIPS  numbers  and  names  of 
states  to  be  included  in  an  assemoled  regional  network.  Users  must  verify 
the  list  of  included  states  by  replying  "OK”  to  the  request: 

TYPE  "OK"  TO  ASSEMBLE  NETWORK 

As  nodes  and  links  are  assembled  for  each  state  from  the  Node  File 
and  the  Link  File,  a  short  message  is  printed  at  the  terminal.  If  more 
than  500  nodes  or  2,000  links  are  found,  an  error  message  is  printed  and 
the  network  assembly  is  aborted.  Occasionally,  a  link  is  read  with 
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TABLE:  State  Directory 


FIPS 

State/District 

FEMA 

FIPS 

State/District 

FEMA 

No. 

Name 

Reg ion 

No. 

Name 

Region 

1 

Alabama 

4 

29 

Missouri 

7 

2 

Alaska 

10 

30 

Montana 

8 

3 

31 

Nebraska 

7 

4 

Arizona 

9 

32 

Nevada 

9 

5 

Arkansas 

6 

33 

New  Hampshire 

1 

6 

California 

9 

34 

New  Jersey 

2 

7 

35 

New  Mexico 

6 

8 

Colorado 

8 

36 

New  York 

2 

9 

Connecticut 

1 

37 

North  Carolina 

4 

10 

Delaware 

3 

38 

North  Dakota 

8 

11 

Dist.  of  Columbia 

3 

39 

Ohio 

5 

12 

Florida 

4 

40 

Oklahoma 

6 

13 

Georgia 

4 

41 

Oregon 

10 

14 

Guam 

9 

42 

Pennsylvania 

3 

15 

Hawaii 

9 

43 

Puerto  Rico 

2 

16 

Idaho 

10 

44 

Rhode  Island 

1 

17 

Illinois 

5 

45 

South  Carolina 

4 

18 

Indiana 

5 

46 

South  Dakota 

8 

19 

Iowa 

7 

47 

Tennesee 

4 

20 

Kansas 

7 

48 

Texas 

6 

21 

Kentucky 

4 

49 

Utah 

8 

22 

Louisiana 

6 

50 

Vermont 

1 

23 

Maine 

1 

51 

Virginia 

3 

24 

Maryland 

3 

52 

Virgin  Islands 

2 

25 

Massachusetts 

1 

53 

Washington 

10 

26 

Michigan 

5 

54 

West  Virginia 

3 

27 

Minnesota 

5 

55 

Wisconsin 

5 

28 

Mississippi 

4 

56 

Wyoming 

8 

1 
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endpoints  which  cannot  be  matched  to  any  location  on  the  node  list. 

When  this  occurs  a  one-line  message  is  printed  and  the  network  assembly 
is  allowed  to  continue.  When  the  data  assembly  is  completed  the  total 
numbers  of  nodes  and  links  in  the  regional  network  are  printed  and  the 
network  is  sorted. 

Up  to  ten  regional  networks  can  be  sorted  in  ready- to-use  form  on 
files  attached  to  channels  11  through  20.  Option  6  stores  an  assembled 
network  and  Option  5  retrieves  one  which  was  previously  stored.  Users 
must  supply  the  program  with  a  channel  number.  When  a  network  is  re¬ 
trieved,  the  program  lists  the  states  included  in  the  region  and  the 
number  of  nodes  and  links  in  the  network. 

Following  the  assembly  or  retrieval  of  a  regional  network,  users 
are  given  an  opportunity  to  list  the  information  pertaining  to  the  popu¬ 
lation  centroids  in  the  region.  The  response  "YES"  or  just  "Y"  to  the 
question: 

LIST  THE  POPULATION  CENTROIDS  (YES  OR  NO)? 
results  in  the  production  of  Table  4.1  presenting  the  data  from  the  node 
file  for  each  centroid. 
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2.4  The  Crisis  Relocation  Plan 

Control  Option  1  directs  the  program  to  obtain  a  new  crisis 
relocation  plan  either  directly  from  the  user  or  by  constructing  and 
solving  a  linear  program  (the  transportation  problem)  according  to 
directions.  A  simplified  flow  chart  for  this  segment  of  the  program  is 
presented  on  the  next  page. 

The  Interactive  sequence  begins  with: 

INPUT  EVACUATION  PLAN  (YES  OR  NO)? 

To  directly  input  an  evacuation  plan  a  user  responds  "YES"  or  just  "Y." 
Any  other  response  is  interpreted  as  a  "NO."  If  the  response  is  "YES," 
the  program  then  requests  from  the  user  a  list  of  numbers,  three  num¬ 
bers  for  each  evacuation/reception  county  pair  in  the  crisis  relocation 
plan: 

FOR  EACH  PAIR  INPUT:  EVACUATION  COUNTY  NO., 

RECEPTION  COUNTY  NO.,  NUMBER  OF  PEOPLE 
The  three  numbers  for  each  pair  must  be  typed  on  a  single  line  and 
separated  by  commas.  The  county  numbers  are  5-digit  state/county  FIPS 
codes.  Example:  The  line 

26163,  26115,  100000 

calls  for  the  relocation  of  100,000  people  from  Wayne  County  (26163) 
Michigan  to  Monroe  County  (26115).  The  program  echo  checks  information 
as  it  is  input.  If  a  nonpositive  number  is  input  for  the  evacuation 
county,  the  result  is  to  delete  the  immediately  preceding  evacuation/ 
reception  pair.  To  terminate  the  input  sequence  for  the  crisis  reloca¬ 
tion  plan,  the  user  types  the  line: 
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999,  999,  999 

The  program  prints  the  whole  plan  and  returns  to  the  master  control 
sequence . 

If  the  crisis  relocation  plan  is  not  input  directly,  the  program 
computes  one  by  constructing  and  solving  the  transportation  problem 
described  in  1.2.  First,  users  are  given  the  opportunity  to  directly 
input  the  evacuation  counties  and  the  number  of  evacuees  from  each 
county : 

SPECIFY  EVACUATION  COUNTIES  (YES  OR  NO)? 

If  the  user  answers  "YES"  or  "Y,"  the  program  interrogates  the  user  as 
follows : 

FOR  EACH  EVACUATION  COUNTY  INPUT: 

COUNTY  NO.,  NUMBER  OF  PEOPLE 
For  each  evacuation  county  the  user  types  the  county's  5-digit  FIPS  code 
and  the  number  of  evacuees  from  the  county  separated  by  a  comma.  For 
example: 

26163,  2000000 

instructs  the  program  to  treat  Wayne  County  (26163)  as  an  evacuation 
county  with  2,000,000  residents  to  be  relocated.  The  data  is  echo 
checked  on  input,  a  nonpositive  county  number  deletes  the  previously 
input  evacuation  county,  and  the  input  sequence  is  terminated  by  typing 
the  line : 

999,  999 

If  the  user  does  not  directly  specify  the  evacuation  counties,  then 
they  will  be  indirectly  specified  by  a  set  of  controls  and  parameter 
values  which  must  be  provided  by  the  user  in  response  to  a  series  of 
information  requests  appearing  at  the  terminal: 


INPUT  VULNERABILITY  GROUP  CODE  (1  TO  4) 

The  user  supplies  the  first  digit  of  FEMA’s  crisis  relocation  county 
classification  system.  All  counties  with  a  first  digit  which  is  less 
than  the  supplied  parameter  value  (except  for  A  counties)  are  considered 
evacuation  counties,  i.e.,  the  response  "3"  has  the  effect  of  treating 
all  B  and  C  counties  as  evacuation  counties. 

INPUT  MINIMUM  PERCENT  TO  BE  EVACUATED 
The  user  supplies  a  number  in  the  range  0-100.  This  percentage  is  ap¬ 
plied  to  each  evacuation  county's  1980  population  to  determine  the 
minimum  number  of  evacuees . 

EVACUATE  COUNTERFORCE  TARGETS  (YES  OR  NO)? 

If  the  user  responds  "YES"  or  "Y"  then  all  counties  which  have  signifi¬ 
cant  counterforce  targets,  as  indicated  by  the  second  digit  of  the  FEMA 
classification  code,  are  added  to  the  list  of  evacuation  counties.  All 
of  the  populations  of  counties  with  counterforce  targets  are  evacuated. 

INPUT  POST  EVACUATION  MAXIMUM  DENSITY  IN 

EVACUATED  COUNTIES  IN  PEOPLE/SQUARE  MILE 

The  number  supplied  is  compared  to  the  normal  population  density  of  each 
f’vacuation  county.  If  the  normal  density  exceeds  the  maximum,  then  the 
number  of  evacuees  required  to  reach  the  maximum  density  is  computed. 
This  figure  is  then  compared  to  the  number  of  evacuees  previously  com¬ 
puted  using  the  minimum  percentage.  The  larger  of  the  two  numbers  is 
retained  as  the  number  of  eva  uees  from  the  county. 

A  similar  sequence  is  followed  for  reception  counties  and  their 
capacity  to  hold  evacuees.  To  directly  input  reception  counties  and 
capacities  the  users  answers  "YES"  or  "Y"  to: 

SPECIFY  THE  RECEPTION  COUNTIES  (YES  OR  NO)? 


33 


For  each  reception  county,  the  user  supplies  the  county's  5-digit  FIPS 
code  and  the  number  of  places  for  evacuees  in  the  county.  For  example: 

26115,  100000 

indicates  that  Monroe  County  (26115)  is  a  reception  county  capable 
of  accommodating  100,000  evacuees.  The  data  is  echo  checked  on  input,  a 
nonpositive  county  number  deletes  the  previously  input  evacuation  county, 
and  the  input  sequence  is  terminated  by  typing  the  line: 

999,  999 

If  the  reception  counties  are  not  specified,  then  a  set  of  controls 
and  parameter  values  are  used  to  identify  reception  counties  within  the 
region. 

INPUT  MINIMUM,  MAXIMUM  PRE- EVACUATION  RECEPTION 

AREA  DENSITY  IN  PEOPLE/ SQUARE  MILE 

Two  numbers,  separated  by  a  comma,  are  typed  at  the  terminal.  The  list 
of  non-evacuation  counties  is  then  scanned  for  counties  whose  1980  popu¬ 
lation  density  falls  within  the  stipulated  limits.  Such  counties  are 
included  on  the  list  of  reception  counties.  For  example,  the  response: 

20,  500 

would  limit  the  selection  of  reception  counties  to  those  with  1980  popu¬ 
lation  densities  between  20  and  500  people/square  mile.  The  lower  limit 
serves  to  exclude  areas  such  as  mountains  and  deserts  which  cannot  sustain 
evacuees  while  the  upper  limit  excludes  counties  which  are  densely  popu¬ 
lated  but  have  not  been  selected  for  evacuation. 

INPUT  MINIMUM,  MAXIMUM  TRAVEL  TIME  IN  HOURS  X  100 


The  response  to  this  request  establishes  a  travel  time  band  around  each 
evacuation  county.  Travel  times  are  measured  along  the  fastest  route  at 
normal  highway  speeds.  As  the  transportation  problem  described  in  1.2 
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Is  constructed  the  travel  times  assigned  to  movements  from 

evacuation  counties  (i)  to  reception  counties  (j)  are  the  normal  travel 
times  between  the  counties.  However,  when  these  normal  travel  times 
fall  outside  the  travel  time  band,  the  travel  times  are  artificially 
increased  by  1,000  hours.  In  effect,  reception  counties  within  the 
travel  time  bands  of  evacuation  counties  are  filled  before  any  use  is 
made  of  reception  counties  which  lie  outside  the  time  bands.  For  example, 
the  values: 


50,  400 

establish  a  band  from  half  an  hour  to  four  hours  around  each  evacuation 
county.  The  program  will  relocate  evacuees  from  an  evacuation  county  to 
reception  counties  within  this  band  until  the  reception  capacity  of  these 
counties  is  expended  before  using  any  other  reception  counties.  Tn  ■ 
lower  limit  (one-half  hour)  pushes  evacuees  beyond  the  Immediate  suburbs 
of  most  cities  while  the  upper  limit  (four  hours)  holds  evacuees  within 
the  range  of  a  one-day  round  trip  from  their  homes  and  places  of  work. 

INPUT  MAXIMUM  REStDENT  POPULATION  MULTIPLE  X  100 

The  maximum  resident  population  multiple  is  applied  to  each  reception 
county's  1980  population  to  determine  the  number  of  evacuees  which  the 
county  can  accommodate.  For  instance,  a  response  of  "300"  is  interpreted 
to  mean  that  reception  counties  can  hold  two  evacuees  for  every  pre¬ 
evacuation  resident.  Thus,  the  capacities  of  the  reception  counties  would 
initially  be  set  at  two  tiroes  each  county's  1980  population. 

INPUT  POST  EVACUATION  MAXIMUM  DENSITY  IN 

RECEPTION  COUNTIES  IN  PEOPLE /SQUARE  MILE 

The  area  of  each  reception  county  is  multiplied  by  the  maximum  density 
supplied  by  the  user  to  obtain  the  maximum  number  of  people  which  the 
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the  county  can  hold  following  an  evacuation.  The  1980  population  is 
then  subtracted  to  yield  the  maximum  number  of  evacuees.  This  maximum 
is  compared  to  the  reception  capacity  computed  using  the  maximum  resi¬ 
dent  population  multiple  and,  if  smaller,  it  becomes  the  county's 
reception  capacity. 

At  this  point  in  the  program,  all  of  the  information  required  to 
formulate  the  transportation  problem  described  in  1.2  has  been  input  and 
processed.  A  short  table  is  printed  at  the  terminal  displaying  the 
total  number  of  evacuation  counties,  evacuees,  reception  counties, 
reception  capacity  and  the  values  of  all  parameters  supplied  by  the 
user.  The  table  listing  is  followed  by  the  line: 

TYPE  "OK"  TO  COMPUTE  EVACUATION  PLAN 
If  the  user  responds  with  an  "OK,"  a  standard  algorithm  is  applied  to 
either  solve  the  transportation  problem  or  to  determine  that  it  is  in¬ 
feasible.  A  short  message  is  printed  at  the  terminal  reporting  the 
outcome  of  the  attempt  to  solve  the  problem.  Any  response  other  than 
"OK"  returns  the  program  to  the  master  control  sequence  without  solving 


the  transportation  problem. 
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2.5  Route  Selection  and  Loading 

Control  Option  2  results  in  the  selection  and  loading  of  routes  for 
a  crisis  relocation  plan  as  described  in  1,3.  A  simplified  flow  chart 
appears  on  the  following  page. 

To  select  and  load  routes  a  certain  amount  of  additional  informa¬ 
tion  about  the  operation  of  the  highway  network  during  an  evacuation  is 
needed.  This  information  must  be  supplied  by  the  user  when  requested 
by  the  program. 

INPUT  VEHICLES /LANE/HOUR  FOR  1)  LIMITED  ACCESS, 

2)  PRIMARY  4-LANE,  3)  PRIMARY  2-LANE,  4)  SECONDARY 
HIGHWAYS 

Four  numbers  are  input  giving  the  emergency  carrying  capacity  in  vehicles/ 
lane/hour  of  four  general  classes  of  highways .  The  capacity  of  each  link 
is  found  by  multiplying  the  number  of  available  lanes  by  the  vehicles/ 
lane/hour  for  the  type  of  highway  comprising  the  link.  Example: 

2000,  1500,  1400,  1200 

Speeds  of  movement  during  an  evacuation  are  input  for  the  same  highway 
classes  in  response  to  the  following  request: 

INPUT  AVERAGE  SPEED  IN  MILES /HOUR  FOR  1)  LIMITED 

ACCESS,  2)  PRIMARY  4-LAME,  3)  PRIMARY  2-LANE, 

4)  SECONDARY  HIGHWAYS 

Normal  nonevacuation  speeds  along  the  highway  network  are  taken  to  be 
55,  45,  40  and  30  miles  per  hour.  Slower  average  speeds  would  likely 
prevail  during  an  evacuation.  For  example: 


40,  35,  35,  30 
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As  discussed  in  1.3  three  options  exist  for  controlling  the  use  of 
highway  lanes  during  an  evacuation.  One  of  these  options  must  be 
selected  by  the  user  in  response  to: 

INPUT  HIGHWAY  CONFIGURATION  COPE  1)  NORMAL, 

2)  ONE  WAY,  3)  VARIABLE  LANES 

The  responses  "1,"  ”2"  and  "3"  correspond  to  the  three  options  described 
in  1.3. 

INPUT  PASSENGERS/VEHICLE  X  100 

The  user-supplied  number  is  needed  to  convert  evacuees  to  vehicles  along 
routes . 

INPUT  NOMINAL  EVACUATION  TIME  IN  HOURS  X  100 
Half  the  nominal  evacuation  time  is  imposed  as  an  initial  minimum 
average  delay  for  all  evacuees  on  all  routes. 

INPUT  DELAY  TIME  TERMINATION  LIMIT  IN  HOURS  X  100 
The  route  selection  and  loading  algorithm  terminates  if  the  delay  times 
on  unblocked  links  are  all  less  than  the  termination  limit  supplied  by 
the  user  or  half  the  nominal  evacuation  time. 

A  short  table  is  printed  at  the  terminal  echoing  the  input  parameter 
values  followed  by: 

TYPE  "OK"  TO  GENERATE  ROUTES 
The  response  "OK"  is  immediately  followed  by  execution  of  the  route 
selection  and  loading  algorithm  described  in  1.3.  Any  other  response 
returns  the  program  to  the  master  control  sequence.  As  the  algorithm 
proceeds  a  short  listing  appears  at  the  terminal  describing  the  progress 
of  the  computations.  The  listing  describes  the  status  of  the  link  list 
after  each  pass  and  the  links  which  are  blocked  prior  to  the  next  pass 
through  the  Ford-Fulkerson  route  selection  routines. 
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2.6  Evacuation  Simulation 

Control  Option  3  directs  the  program  to  simulate  an  evacuation 
along  the  regional  highway  network  according  to  a  previously  computed 
set  of  routes  and  loadings.  No  additional  input  is  required  from  the 
user.  At  the  conclusion  of  the  simulation,  a  short  table  is  printed  at 
the  terminal  showing  the  percentile  distribution  of  evacuees  as  at  risk 
(unevacuated),  en  route,  or  at  their  destination. 


3.0  Sample  Run 

A  sample  run  using  FEMA's  UNIVAC  1108  computer  illustrates  the 
interactive  usage  of  the  IDA/BPT  model  and  provides  examples  of  the 
variety  of  tables  that  can  be  generated  by  model  users.  There  are 
three  parts  to  the  sample  run: 

1.  The  on-line  (terminal)  input/output  stream. 

2.  Summary  tables  output  through  channel  20. 

3.  Route  listings  output  through  channel  19. 
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3.1  Terminal  Input/Output  Stream 

The  sample  run  Is  a  straightforward  application  of  the  model  to  a 
small  but  typical  problem,  the  evacuation  of  the  major  population.  In¬ 
dustry  and  military  centers  of  New  England.  The  terminal  listing  has 
been  annotated  to  show  the  major  elements  of  the  computations: 

1.  The  control  options  are  listed  as  the  program  begins 
execution. 

2.  Route  listings  are  routed  to  channel  19  which  Is  linked 
to  the  file  FILE  19. 

3.  Summary  tables  are  routed  to  channel  20  which  Is  linked 
to  the  file  FILE  20. 

4a.  The  New  England  regional  transportation  network  is  retrieved 
in  preprocessed  form  through  channel  11.  Previously  it 
had  been  stored  on  FILE  11,  or, 

4b.  The  New  England  network  is  assembled  from  the  transpor¬ 
tation  network  data  base. 

5.  Table  4.0  is  printed. 

6.  The  controls  and  parameter  values  for  the  crisis  reloca¬ 
tion  planning  submodel  are  input. 

7.  Table  1.0  is  printed;  the  transportation  problem  is  solved. 

8.  The  controls  and  parameter  values  for  the  route  selection 
and  loading  submodel  are  Input. 

9.  Table  2.0  is  printed;  the  route  selection  and  loading 
algorithm  is  begun. 

10.  Table  2.1  Is  printed  as  the  algorithm  progresses.  Each 
line  of  the  table  shows  a  link  to  be  blocked  at  the  next 
iteration  of  the  route  selection  procedure.  Also  shown 
is  the  size  of  the  route  list  after  each  pass. 

11.  The  evacuation  Is  simulated;  Table  3.0  is  printed. 


12.  End. 


42 


LIST  OF  CONTROL  OPTIONS 
NO.  CONTROL  OPTION 

0  REPEAT  LIST  OF  CONTROL  OPTIONS 

1  INPUT/COMPUTE  CRISIS  RELOCATION  PLAN 

2  RON  RObTE  FINDINij/LOADING  ALGORITHM 

3  SIM'.jOATE  the  EVACUATION 

4  ASSEMBLE  NETNCRK  FROM  DATA  FILES 

5  R.EAO  _3T;:?.ED  NETWORK 

i  Ni-NERATEO  NETWORK 

'  ;  LrANNiiL  rvR  PARAMETER  VALLES 

3  REASSIGN  OLTRLT  FOR  i:--.>IMARY  TAElES 

;Sass:on_cutf';t  i-annel  for  route  _;st;;.o 

1  tiRMilNATE  FRuGRAM 

2  INPUT  CONTROL  CFTION  NUMSER 

9 

CONTROL  OPTION  9  TYPE  "OK"  TO  EXECUTE 
OK 

INPJT  NEW  CHANNEL  NUMBER  (6=T£RMINAL) 

19 

■1  INPLIT  CCNTRCX.  OPTION  NUMBER 

3 

CONTROL  OPTION  8  TYPE  "OK"  TO  EXECUTE 

INPUT  NEW  CHANNEL  NUMBER  (6=TERMINAL) 

20 


4a  INPtJT  CONTROL  OPTION  NUMBER 

COfTTROL  OPTION  5  TYPE  "CK"  TO  EXECUTE 
iDK 

ASSIGN  INPUT  CHANNEL  FOR  SAVED  NETWORK 
11 


NUMBER 

STATE  NAME 

9 

CONNECTICUT 

23 

MAINE 

y’l 

MASSACHUSEHS 

33 

NEW  HAMPSHIRE 

44 

RHODE  ISLAND 

50 

VER:M.ONT 

NUMBER  OF  NETWORK  NODES 

170 

NUMBER  OF  NETWORK  LINK,S 

788 

LIST  THE  POPULATION  CENTROIDS 

(YES  OR  NOP 

YES 

INPUT  D'JNTROL  OPTION  NUMBER 

4 

CONTROL  OPTION  4  TYPE  "OK*  TO  EXECUTE 
OK 

INPUT  THE  NUMBER  OF  STATES  IN  THE  REGION  (UP  TO  101 

5 

INPUT  6  STATE  NUMBERS 
9,23.25.33.44,50 


5  NUMBER  STATE  NAME 
9  CONNECTICUT 

23  MAINE 
25  MASSACHUSETTS 
33  NEW  HAMPSHIRE 

44  RHODE  ISLAND 

50  VERMONT 

TYPE  "CK*  TO  ASSEMBLE  NETWORK 
OK 

READING  37  NODES  FOR  STATE  9 
READING  25  NODES  FOR  STATE  23 
READING  39  NODES  FOR  STATE  25 
READING  33  NODES  FOR  STATE  33 
READING  15  NODES  FOR  STATE  44 
READING  21  NODES  FijR  STATE  50 
READING  200  LINKS  FOR  STATE  9 
READING  U4  LINKS  FOR  STATE  23 
READING  224  LINKS  FOR  STATE  25 
READING  184  LINKS  FOR  STATE  33 
READING  74  LINKS  FOR  STATE  44 
READING  11.8  LINKS  FOR  STATE  50 

NUMBER  OF  NETWORK  NtDDES  170 

NUMBER  OF  NETHORIC  LINKS  733 
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IMFUT  CONTROL  OPTION  NUMBER 
ICNTRCl  option  1  TYPE  'OK”  TO  EXEiCUTE 
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INPUT  EvACi.ftTICN  FLAN  Ut'ES  NO)^ 

NO 

-El  IF.  E.AC.ATICN  C  CONTI  EE  i'^S  CR  140 
.  ’  •  _'i£'AEI^ITY  Iri.t.ir  COliE  U  TO  41 

^E^'iENT  ee  evacuated 
{.(jljti’E  .^vNTEFF,jRf£  TAPvETE  'VES  OR  NO)"' 

:  A'.'  -Ur  E.hOwATIGN  .'WXl.’IJfl  IOISITY  in  EVACliATEi  COIWTIES  IN  PECPLE/EO.  MILE 

c  ■  ^ 

IfECIF'  FECEFTION  counties  (YES  OR  MO)’ 
nC 

IMPJ  NINIHUM,  MAXIMUM  PRE-EVACUATION  RECEFTICN  AREA  DENSITY  IN  PEOFtE/EW.  MILE 

■  , '  ^  j-y 

isPj'MINirUI.  MAXIMUM  TRAVEL  TIME  IN  HOURS  X  100 
Ilf^T  ■WXIMOM  RESIDENT  POPULATION  MULTIPLE  X  100 

'nfj-'  ’losT  p;(n:uation  maximum  density  in  reception  colinties  in  people/sq.  mile 

soo 

niimHT  if  evacuated  counties  11  evacuated  population  26395*116. 

NLIMEER  OF  reception  COUNTIES  41  RECEPTION  AREA  CAPACITY  4251722. 

.ul.‘£RABIl:ty  group  coce;  3 

MINIMUM  EVACUATED  Pi'jPULATIONt  90.  PERCENT 

COUNTEPfORCE  TARGETS  ARE  EVACUATED 

MINIMUM  FCCEF'TION  fKlP,  DETfSITY:  10.  PEOFlE/SQ  MILE 

MAXIMUM  f€CEFTION  AREA  DENSITY:  300.  PEQFLE/SQ  MILE 

MINIMUM  TRAVEL  TIME  TO  RECEPTION  AREAS'.  .50  't+DURS 

MAXIMLW  TRAVEL  TIME  TO  FfiCEPTION  AREAS:  8.00  HOURS 

MAXIMUM  RESIDENT  r-iOPULATION  MULTIPLE:  3.00 

POST  EVACUATION  MAXIMUM  DENSITY  IN  EVADJATED  AREAS:  500. 

POST  PVAIUATION  MAXIMUM  DENSITY  IN  RECEPTION  AREAS:  500. 


TYFE  "OF"  TC  COHPlfTE  EV.ROJATION  PLAN 
OK 

EVACUATION  PROBLEM  IS  FEASIBLE 

NUMBER  OF  EVACUATION/RECEPTION  COUNTY  PAIRS 
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8  INPUT  CONTROL  OPTION  NUMBER 

CONTROL  OPTION  2  TYPE  "DK"  TO  EXECUTE 
OK 

INPUT  VEHICLES/LANE/WOLR  FOR  DLIMITED  ADOESS,  2)PRIMARY  4-LANE,  3)PRIMAfiY  2-LANE,  4)SEKNDARY  HIGHWAYS 
2000, 1600. 1400. 1’iXl 

INPtJT  AVERAGE  'SPEED  IN  MILES/HOUR  FOR  DLIMITED  ACCESS-  2iPRlMARY  4-LANE,  3)FRIMARY  'Z-LANE,  4)SED:NCARY  HI 
45 . 40  T  5 j . UO 

INPUT  HIGHWAY  CCNFlGURATUjN  CODE  DNORHAL,  2)0NE  WAY.  3) VARIABLE  LANES 

iNrUT  FiASSENGERS/VEHiCLE  X  100 
400 

INPiJT  NOMINAL  EVAaWTION  TIME  IN  HOURS  X  100 
300 

INPIJT  DELAY  TIME  TERMINATION  LIMIT  IN  HOURS  X  100 
600 

9  HIGHWAY  CLASS  VCLS/LANE/W  MILES/HOUR 

LIMITED  ACCESS  2000.  45. 

PRIMARY  4-LANE  1600.  40. 

PRIMFiRY  2-LANE  1400.  35. 

SECONDARY  1200.  30. 

HIGHWAY  CONFIGURATION:  ONE  WAY  OUTBOUND  TRAFFIC 
AVERAGE  PASSENGERS  FER  VEHICLE:  4.00 
NOMINAL  EVACUATION  TIME:  8.00  HOU«S 

DELAY  TIME  TERMINATION  LIMIT:  6.00  HOURS 


TYPE  ‘OK"  TO  GENERATE  ROUTES 
OK 


19 

:c 


41 

41 


44 


47 


IQ  47 


JACITY 

'VEHICLES 

DELAY 

1400. 

7S053. 

27.33 

‘  4C'<T. 

7xi 

27.33 

140'j. 

'T-iCcrc 

15.91 

ir  r»’i . 

147  10. 

13,23 

IRIS'). 

41422. 

14.7’’ 

l4r;0. 

41422. 

14.79 

1400. 

1 1425. 

4000. 

12.  w4 

.  - 

Ai'uO . 

9. 

iiin 

'  <*,  •  **  1' 

10.42 

4i'iuO. 

7s;'2'-i. 

•9.34 

\4i'ii). 

2'^ ’7'^. 

10.42 

4000. 

73729. 

9.34 

4000. 

737. ;9. 

9.34 

411OO. 

/ 

9.25 

4000. 

7J-T36. 

9.25 

1400. 

242)2. 

3.64 

1400. 

232^4. 

0 

1400. 

2'3'2'-'4. 

V. 

HOC. 

23294. 

3.32 

1400. 

21442. 

7.66 

1400. 

21524. 

7.69 

4000. 

611'21. 

7.64 

1400. 

23294. 

3.32 

400IJ. 

57642. 

7.21 

4000. 

57'642. 

7.21 

14IX). 

23415. 

7.29 

4iyXN 

57642. 

7.21 

4000. 

56932. 

7.12 

4000. 

5;35.'ii. 

7.32 

4000. 

535  A 1. 

7.32 

41:100. 

56'932. 

7,12 

1100. 

171-5. 

7.14 

4000. 

56932. 

7.12 

4000. 

569:32. 

7.12 

140'i. 

17V5. 

Y  “» 

C  t  1  4. 

1400. 

17135. 

6.1; 

1400. 

I:S0:S1. 

6.44 

4T)(K). 

51665. 

6.46 

1400. 

I'SOSl. 

6.44 

140C1. 

13311. 

6.54 

40iM. 

4444.2. 

5.56 

5.56 

4000. 

44462. 

4000. 

5.56 

4i::0';i. 

444;,;. 

t,ZQ 

1200. 

13'?:  3. 

5.30 

1200. 

13323. 

5.80 

1200, 

132611. 

5.53 

1200. 

l.SMi. 

5.30 

400(J. 

4354.6. 

5.4+ 

1400. 

14'’03. 

5.25 

FROM 

AVON 

riASTFORD 

LEE 

•!Ai9CHEbrEF 
JiXiTH  FOSTER 
.  I'HNSrON 
1A..F.:40 
5£  j'r.EKrC'*I'l 

sos*;« 

■:pc  :>;r;F  IELQ 

a.0.vL-£ 

LKLRE*>tf^E 

:.RAw:€ 

SLATERSVILLE 

WORCHESTES 

ESMOND 

ASTON 

BIDDLEFORD 

PCiRTSMOMTH 

NDRTHHAMPL3N 

AS;HFCRO 

PUTNAM 

CHEF'AiSHET 

NORTHFIELO 

MANCHESTER 

NASHUA 

ESMOND 

MARLBORO 

NORWICH 

MERIDAtJ 

COLCHESTER 

BOSTON 

CHICOPEE 

SPRINfSFIELD 

DANVERS 

WICrtOiOON 

HAMPTON 

AMESBURT 

C'ARCNcR 

(•ARC  HES  TEH 

NEWPORT 

CCICCRD 

NEW  LONDON 

NEWPORT 

TORRINGTON 

MIL:.DAL£ 

WATERBORY 

TKCMASTOWN 

FETERBOHO;,iGH 

NASHUA 

MANCHESTER 

MIlFCRD 

HOFKINTOWN 

BLRLINGrw 


CC4JMTY 

HARFCRD 

HARTFCRD 

BERK.iHIRE 

HILLSoCRCLiOH 

FPuVICENCE 

PROVIDENCE 

''CRL 

HAMPSHIRE 

SoP-'Clc  PSD 

-AMF'[En_, 

LCRCc'jrLR 

cSSEi 

-RAN'-'LIN 

rPOVIDENCE 

WGRCESTSR 

PROVIOEMI^E 

PRiOVlDEfCE 

ycrk: 

RUSKINGWYI 

hamp-shire 

1.1NQHAH 

UINOHAM 

PROVIDENCE 

MERRIMACK 

HILLSBCWCUGH 

HILLSSOROOiSH 

PROVIDENCE 

HARTFORD 

NEW  LONDON 

NEW  rtAVEN 

NEW  LONDON 

Suffolk:  psd 

HAMPDEN 

HAMPDEN 

ESSEX 

WORCESTER 

ROCKINGHAM 

ESSEX 

WORCESTER 

ViCF.CESTER 

iLlilVAN 

MERRIMACK 

MERRIMACK 

NEWFTiRT 

LITCHFIELD 

rART'ORC 

NBJ  ^A'.'EN 

LITCHFIELD 

HILLSBOROUGH 

H1LL^P0'''0‘'C'H 

HIUHPOSCri'SH 

HILLSeOiSOiJGH 

MERRIMACK 

CHITTENDEN 


TO 

WINSTEAD 

AVVN 

SCTTSFIELD 

CCNCCPC 

DANIEL-ON 

SCiOTH  FOSTER 

r'C'RTLAND 

NCSTHHAMPTG‘1 

LhARENLE 

C-ICCPEE 

‘  tWN 

kanche;*er 

:-EENFISi.C 

WL'HiCHEiTER 

GRANGE 

ASTON 

SLATERSVILLE 

FCRTlAND 

BIDDLEFORD 

PITTSFIELD 

BCiLTON 

ASHFORD 

PUTNAM 

LACONIA 

CONCORD 

MANCHESTER 

CHEPACHET 

HARTFCiRO 

CCLCHESTER 

M1U.DALE 

MARLBORO 

DANVERS 

iNORTHHAMP'TCN 

CHICiPEE 

AMESBClRY 

KEENE 

PC^TSMiQUTH 

HAMPTON 

W  iCHEHOON 

oARTNER 

CLAPiM'CN’’ 

NORTHFIELD 

NEWPORT 

SAUNOERSTOWN 

WINSTEAD 

WATESSLRV 

THCMASTCWN 

TCiRRlNOTON 

KEENE 

MILFORD 

HENN’.'oR 

PETERBORCiiX'H 

NEW  LONDON 

VERGENNES 


county 

.I'CrFiaD 

.-ihH’T-CSD 

SEFTSHiFE 

MKRIMACK 

W.vEHAM 

SSOVIOENCE 

CLMSERLhI.D 

r-AMPSHIRE 

ii:EX 

-h.”F0£n 

rrM^CrCSE  _ 
hI_L:S'.P.LUL’H 

(iC'RCESTER 

FRANKLIN 

PROVUENCE 

FRCiVICENCE 

CUMBERLAND 

YORK 

BERKSHIRE 

TOLLAND 

WINDHAM 

WINDHAM 

BELKNAP 

MrKRlMPCK 

hILLSBCROUOH 

PROVIDENCE 

HARTFORD 

NEW  l['NC':in 

hartfurd 

HARTFORD 

ESSEX 

HAMPSHIRE 

HAMPOEN 

ESSEX 

uHEEHi.RE 

ROCKINGHAM 

ROCKINGHAM 

WORCESTER 

«CRCES  ^E' 

SOL:.  Ivan 

MERRIMACK 

SULLIVAN 

UASHINIjTQN 

LITCHFia.O 

SEW  -A'.'EN^ 

;.ITCfFi£i-T 

LITCHFIELD 

CHESHIRE 

HILLS&CROUGH 

MERRi:''H:< 

HiLLS&i.iRi.ii.iijH 

MERRIMACK 

ADDISON 


44 

44 


ROUTE  NAME 
J«  •!>*  nUU  1 E 
^oOtE 
U.S.  ROUTE 
IirERSTA^; 
•• 

U.S.  kC'E 

-'TATE 
iTTE.i’’.: :  fM^'E 

2  '  n  Z  r  , 

*  •  ■Trii-c'T''' 

;  r, .  *  I  r  ‘  il 


:STATE  ROAD  122 
INTERSTATE  295 
STATE  RCAD  146 
■NT-R.STA'"E 
INTERSTATE 
STATE  ROAD 
U.S.  ROUTE 
iJ.S.  ROUTE 
U.S.  ROUTE 
U.S.  ROUTE 
U.S,  ROUTE 
STATE  RCAD  SO' 
U.S.  RC'TE 
STATE  ROAD 
STATE  RCAD 
5TA'’r:  RCAD 
STATE  ROAD 
IS'E'S'A'E 
IN'E-SFA'^E 
U.S.  ROUTE 
INTERSTATE 
STATE  ROAD 
INTERSTATE 
IN'I’ERS'^A'^ 


95 

9 

44 

44 

44 


44 


66 

95 

91 

e; 

12 

95 

95 


S'A'E  RcAD  140 
STaK  R'C'AC  14u 
11 
93 
11 


INTERSTATE 
■STATE  RCAD 
STATE  ROAD  133 
STAT,:  ROAD  3 
ISTEPSTATE  P4 

sthie  road  s 

STATE  ROAD  101 
■STATE  fiOA'D  102 
:5TAT£  RCAD  114 
SCA  ■.  xuAD  101 
INTERSTATE  6'9 
U.S.  R'jUTE  7 


11 
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INPUT  CONTROL  OPTION  NUMBi-l 

CONTPjl  -IPTION  S  Type  "OK'  '}  EXECUTE 
OK 

EVAiCliEEN  2(v5''504.  ELAPSED  HOURS: 
PERi'.l  Nf-i'E  I-  rxCUEES  AT 
PERCENTAGE  OF  EVACUEES  F-N 
PERCENTAGE  OF  EVACUEES  AT  Ct:-.-  ‘XT ION 

INPIT  CCiNTROL  OPTION  I•I1BE■ 

99 

TYPE  “OK*  TO  TERMOR*'-. 

OK 

PROijRAM  ENDS 


7^ 

15 

4 


4 

5P 

17 

23 


6 

38 

17 

43 


o 

13 

17 

64 


12 

16 

20 

24 

4.) 

48 

oO 

0 

0 

0 

0 

'0 

0 

0 

4 

0 

0 

■) 

0 

0 

)•) 

6 

94 

99 

■59 

99 

99 

99 

9'9 

9<? 

,■) 

;Y 

96 


'-;4 

0 

A 

99 


3.2  Summary  Tables 

The  summary  tables  1.0,  2.0,  2.1,  3.0  and  4.0  reproduce  tables 
from  the  on-line  input/output  stream  and  are  not  duplicated  below.  All 
of  the  other  summary  tables  from  the  sample  run  are  shown  below.  Variable 
definitions  are  given  in  the  notes  accompanying  each  table. 

Notes  for  Table  4.1 

FIPS  -  five-digit  state/county  FIPS  code 

PLACE  NAME  -  name  of  county  center  or  principal  city 

COUNTY  -  county  name 

POPULATION  -  residents  according  to  the  1980  U.S.  Census 
AREA  -  square  miles  of  area 
DENSITY  -  pre-evacuation  people /square  mile 
FEMA  -  FEMA  Countervalue  code  A,  B,  C  or  D 
FEMA  Counterforce  flag  X 

blank  field  Indicates  population  not  at  risk 
Notes  for  Table  1.1 

EVACUATION/RECEPTION  COUNTY  -  county  name 
RESIDENTS  -  pre-evacuation  population 
AREA  -  square  miles 

DENSITY  -  pre-evacuation  people/square  mile 
AT  RISK  POPULATION  -  number  of  evacuees 

REC.  CAPACITY  (evacuation  counties)  -  capacity  of  reception  counties  within 
preferred  travel  time  band  of  the  evacuation  county 
CRITERIA  -  criteria  for  determining  capacity  of  reception 
county 

REC.  CAPACITY  (evacuation  counties)  -  capacity  of  reception  county 
to  accommodate  additional  population 

Notes  to  Table  2.1 

PAIR  -  pair  number  ordered  by  normal  travel  time 
PEOPLE  -  number  of  evacuees 
ORIGIN/DESTINATION  -  name  f  county  center 
COUNTY  -  county  name 

HOURS  -  hours  from  origin  to  destination  at  normal  speeds 
along  fastest  route 
MILES  -  length  of  fastest  route 

DENSITY  -  post-evacuation  density  of  the  reception  county 


Notes  to  Table  2.2 


ROUTE  -  route  list  sequence  number 

PAIR  -  evacuation/reception  pair  number 

EVACUEES  -  loading  of  the  route 

CITY  -  name  of  county  center 

COUNTY  -  county  name 

DISTANCE  -  length  of  route  in  miles 

TRAVEL  -  average  travel  hours  at  evacuation  speeds 

DELAY  -  average  delay  hours 

TOTAL  -  TRAVEL  plus  DELAY  time 

Notes  to  Tables  3.1  and  3.2 

EVACUATION /RECEPTION  COUNT i  -  county  name 
EVACUEE  -  number  of  people  evacuated  from  or  to  the  county 
PERCENTAGE  -  percentage  of  evacuees  departed  from  or  arrived 
at  the  county  at  the  elapsed  hour  shovm 

Notes  to  Table  3.3 

FROM/TO  -  name  of  county  center 
ROUTE  NAME  -  highway  type  and  number 
DISTANCE  -  length  of  the  link  in  miles 
SPEED  -  evacuation  miles/hour 

CAPACITY  -  evacuation  vehicles/hour  for  all  available  lanes 
PERCENTAGE  -  percentage  of  capacity  in  use  for  the  evacuation 
plan  at  the  elapsed  hour  shown 
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1TABI.II  4.1 


F'FS 

PLACE  NAME 

COiTi:/ 

P'DROLATIilii 

AREA 

DENSITY 

FEMA 

"001 

DANBLIRV 

FAIRFIELD 

S07143. 

A  /'t' 

-  ’"’S 

4C0? 

iiARTPORD 

hA'^C-ORU 

‘:077cc. 

7^! 

j  1  LnD 

-ITIiFIElD 

15oR.'. 

*  i."* 

"007 

Ci.3  IA'EPI'Ik 

MI  !3LE‘-t( 

iiv! 

^lA'Th  ~A~EM 

NEW  HAVEN 

"  .  1  DT  * 

iC  ’ . 

'  "c:,. 

■* 

■^oil 

N0RWI3H 

:€U  LONDON 

23S40?! 

■ 

:  X 

0013 

BOLTON 

Ttl,/.A;i7 

114823. 

ii;.. 

-  7  i 

30  1j 

DANIELSON 

WiriDHAM 

C<  /, 

7?  r 

23001 

L&4ISTCN 

ANDR'jSi.;.)GGtN 

99857! 

473! 

20:4. 

D 

23003 

HOULTON 

ARiJOSrOOK 

91331. 

6305. 

13. 

D  X 

23005 

PORTLAND 

CliMfif.R,.AND 

^i^ji  'J-  « 

331. 

245. 

U 

23007 

WILTON 

FRANKLIN 

27093. 

1717. 

16. 

2300? 

ELLSWORTH 

HAriCD';:< 

41731. 

1542. 

27. 

23011 

AU335TA 

KE'lNIiiEC 

10933'9. 

365. 

127. 

23013 

ROCKLAND 

KNOK 

32-94 1. 

362. 

91. 

23015 

NEWCASTLE 

LINCOLN 

25891. 

457. 

56. 

23017 

NEWRY 

OXFORD 

43-?8;3. 

2<X35. 

23. 

23019 

BANGOR 

PENOHSCOT 

137015. 

3408. 

40, 

D 

23021 

FOXCROFT 

P^mTAiIKlIS 

17634. 

3948. 

4, 

23023 

BATH 

SAijAOAHoC 

23795. 

257. 

112. 

23025 

sncwhegan 

SOMErlSET 

45023. 

3943. 

11. 

23027 

BELFAST 

WALDO 

23414. 

7:34. 

39. 

23029 

HACHIAS 

WASHWjTON 

34963. 

2553. 

14. 

D 

23031 

SANFORD 

YORK 

139666. 

ICCO. 

140. 

D  X 

25001 

SAGfSIORE 

BARWiTABLE 

147925. 

3-99. 

:371, 

25003 

PITTSFIELD 

BERKSHIRE 

145110. 

942. 

154. 

D 

25005 

TAUNTON 

BRISf.JL 

474641. 

556. 

854. 

D 

25007 

EEGARTOUN 

DUKES 

8'?42. 

106. 

34. 

25009 

DANVERS 

ESSEX 

633632. 

500. 

1267. 

D 

25011 

GREENFIELD 

FRANKLIN 

64317. 

’?? 

89. 

25013 

SPRIM3FIELD 

HAMPIIEN 

44.3018. 

62i! 

713. 

C  X 

25015 

NORTHHAMPTON 

HAMPSHIRE 

l;38313. 

51’. 

253. 

25017 

LITTLETON 

1367034. 

;329. 

lo4>! 

D 

25019 

NAflTOCKErr 

mnKKET 

1:-. 

1  *  1 

D 

25021 

FRANKLIN 

NCfEi-; 

60653?! 

393. 

1524! 

25023 

PLYMOUTH 

PLYMOUTH 

405437. 

M. 

611. 

D 

25025 

BOSTON 

SUFFUIJI  PSD 

650142. 

55, 

1132'. 

B 

25027 

DORCHESTER 

WCROiSTER 

6463.52. 

1532. 

422. 

D 

33001 

LACONIA 

SE.'.UAP 

42384. 

401. 

107. 

33003 

CONWAY 

CARROLL 

279:31. 

933. 

30. 

33005 

KEENE 

CHE'-KllfE 

62116. 

717. 

37. 

33007 

GORIWM 

COOS 

3514;>. 

1  o  ’e 

1',". 

33009 

PLYMOUTH 

GRAFTjt; 

65306. 

1717. 

o8« 

33011 

NASHUA 

HILLSBOROUGH 

27660:3. 

390. 

311. 

C 

33013 

CONCORD 

MERRIM^::. 

98;;(-.. . 

■>31. 

ns. 

33015 

PORTSMOUTH 

ROLK INGHAM 

190345. 

691. 

275. 

D  X 

33017 

ROCHESTER 

STRAFFiJliO 

854013. 

377. 

’  ’ 

33019 

aAREMONT 

SUU.IVAN 

36063. 

537. 

67. 

44001 

BRISTOL 

BRISTOL 

46942. 

35. 

1  -1 

44003 

EAiST  GREENWICH 

KENT 

15416:3. 

172. 

896. 

44005 

NEWPORT 

NEWP'iRf 

8138;. 

115. 

703. 

C 

44007 

ESMOND 

PRWlOENi.E 

571349. 

422. 

1354, 

B 

44009 

WYOMING 

UASHI'^K  “I 

93317. 

324. 

D  X 

50001 

UERGENNES 

AOlU'  i  r-l 

29406. 

735. 

37. 

50003 

BENNINGTON 

BENNINOTO*: 

33:345. 

672. 

50. 

50005 

SAINT  .^HNSBURG 

CALEnoNI.'H 

2581X3. 

614. 

42. 

50007 

BURLINGTON 

CHITTENDEN 

1155:34. 

532. 

217. 

C 

50009 

ISLAND  POND 

ESSEX 

6313. 

664. 

10. 

50011 

SAINT  ALBANS 

FRTWIIN 

34788. 

659. 

53. 

50013 

ijRAND  ISLE 

GRAND  ISLE 

4613. 

77. 

60. 

50015 

MORRISVILLE 

LAMOILLE 

16767. 

475. 

35. 

50017 

WEST  TOPSHAM 

ORANGE 

22739. 

690. 

33. 

50019 

NEWPORT 

ORLEANS 

23440. 

715. 

3:D. 

50021 

RUTlilND 

RUTLPND 

5C3347. 

^29 

63. 

50023 

MONTPELIER 

WASHINGTON 

52393. 

703. 

74. 

50025 

BRATTLEBARE 

WINDHAM 

36933. 

793. 

47. 

50027 

WILDER 

WINDSOR 

51030. 

965. 

53. 
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irCLt 


EvACJATIi-iu  LuLfHTY 

AREA 

DELSITV 

AT  RISK  ^CP. 

RELi  i-hr'A'. 

,'JEW  LONDjN 

L'3S40’?. 

672. 

654.78 

AROOSTOOK 

91I3L. 

6:3:35. 

13.42 

YCF.I; 

139666. 

li  C'  .'. 

139,67 

125699. 

4i2057£ 

Hf-f>R3EN 

443013. 

621. 

713.39 

393716. 

4093977. 

SUFFOLK  PSD 

650142. 

55. 

11320.76 

622642. 

4223664. 

HlLLSBijRilUGH 

21hl-M. 

S'90. 

310. :30 

2^:5947. 

P226664. 

ROCKIMGHAM 

L90345. 

691. 

275.46 

171510. 

4120572. 

NiLPORT 

aim. 

:i5. 

707.68 

73245. 

4223664, 

PROVIDENCE 

571349. 

422. 

1653,91 

514214. 

4222P>.4. 

WASHINGTON 

■:'3317. 

324. 

2-38.02 

63935. 

4223664. 

CHITTENDEN 

n5S34. 

532. 

217.17 

103931. 

4144362, 

RECEPTION  OOUWTY 

RESIDENTS 

AREA 

[€NSITY 

CRITERION 

REC.  CAPACITY 

LITCHFIELD 

156769. 

933. 

167.13 

DENSITY 

312231. 

TOLLAND 

114823. 

416. 

276.02 

DENSITY 

93177. 

WINDHAM 

92312. 

516. 

17:3.90 

DENSITY 

165663. 

ANDROSCOGGIN 

9'?657, 

473. 

208.49 

DENSITY 

139343. 

CUMBERLAND 

21573'?. 

331. 

244.94 

I.LNSITY 

224711. 

FRANKLIN 

2im. 

1717. 

15.73 

MULTIPLE 

54196. 

HANCOCK 

41731. 

1542. 

27,10 

MULTIPLE 

33562. 

KENNEBEC 

109339. 

365. 

127.04 

MULTIPLE 

219773. 

KNOX 

22941. 

362. 

91,00 

MULTIPLE 

65S82. 

LINCOLN 

25691. 

457. 

56.22 

M).LTIFL£ 

513S2, 

OXFORD 

43?6:3. 

2085. 

23.49 

MULTIFtE 

97936. 

PENOBSCOT 

137015. 

3403. 

40.20 

MULTIFIE 

274<>30. 

SAGADAHOC 

23795. 

257. 

112.04 

MULTIFIE 

57590. 

SOMERSET 

45d2S. 

3943. 

11.41 

MULTIPLE 

9lY'56. 

WALDO 

2S4'4. 

734. 

-'$•  •’  * 

MULTIPLE 

5632S. 

WASHINGTON 

3496-3. 

2553. 

13.o9 

MIL'IFLE 

'  926. 

BERKSHIRE 

145110. 

942. 

154.04 

MULTIPLE 

290220. 

DUKES 

3942. 

106. 

:34.36 

MLILTIPLE 

17;::S4. 

FRANKLIN 

64317. 

722. 

39.03 

MLLTIPLE 

128634. 

HAMP!5HIRE 

I.3!3:313. 

537. 

253.50 

DENSITY 

129637. 

NANTUCiair 

5037. 

46. 

110.59 

MULTIPLE 

10174. 

BELKNAP 

42334. 

401. 

106.94 

MUL’IFLE 

35763. 

CARRQlL 

17931. 

933. 

29.78 

MULTIPLE 

55S62. 

CHESHIRE 

62L16. 

717. 

36.63 

MULTIPLE 

124232. 

CCOS 

35147. 

1825. 

19.26 

HLUIFLE 

70294. 

GRAFTON 

65:306. 

1717. 

3:3.33 

MULTIPLE 

131612. 

MERRIMACK 

93302. 

9:51. 

105.59 

MLLTIPLE 

1966C4. 

STRAFFORD 

35405. 

377. 

226.55 

DENSITY 

103092. 

SULLIVAN 

36063. 

537. 

67.16 

njLTIFLE 

72126. 

ADDISON 

29406. 

735. 

37.46 

MULTIPLE 

58:312. 

BENNINGTON 

33345. 

672. 

49.62 

MULTIPLE 

666 9'3. 

CALEDONIA 

2530:3. 

614. 

42.03 

MULTIPLE 

51616. 

FRANKLIN 

.347:33. 

659. 

52.79 

;-llt;ple 

65576. 

GRAND  ISLE 

4613. 

77. 

59.91 

MLILTIPLE 

9226. 

LAMOILLE 

16767. 

475. 

35.50 

MULTIPLE 

33534. 

ORANGE 

22739. 

690. 

'i',  0, 

"v 

MLLTIPLE 

i5473. 

ORLEANS 

23440. 

715. 

32.78 

MLLTIPLE 

46330. 

F'JTLAND 

53347. 

929. 

62.31 

MULTIPLE 

116694. 

WASHINGTON 

52393. 

708. 

74.00 

MLLTIPLE 

104736. 

WINDHAM 

36933. 

793. 

46.57 

MULTIPLE 

73:366. 

WINDSOR 

51030. 

965. 

52.83 

MULTIPLE 

102060. 
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ITABLE  1.2 


EVAC'jaTlON/RECEPTION  CQIJNTY  PAIRS 
PAIR  PPQPLF.  ORIGIN 

COUNTY 

DESTINATION 

CCyNn 

-'OURS 

MILES. 

i 

22637. 

E:3nCND 

FHO-VIDENCE 

BDlNIN&TuN 

3B«MINGT0N 

2.91 

1-jU. 

-i""'*  r“3' 

n 

575'1'0. 

BttSTQN 

SUFFCU;  P-3D 

BATH 

StYSACiAHOC 

2.44 

134. 

Cf 

13644. 

BOSTON 

aiFiXK  RiD 

LBllSTON 

ANEWIiSCOGGIN 

2.36 

ISO. 

500.00 

4 

1367S. 

NEWPORT 

NEIJ'ORT 

WINSTEAD 

LITCKFIP:  n 

2.S4 

115. 

SCO. CO 

5 

102060. 

BOSTON 

SUFFOliS  PSD 

WILBER 

WINDSOR 

2.23 

1^5. 

153.44 

6 

72126. 

BOSTON 

SUFFOLK  PSD 

CLAREMONT 

SULLIVAN 

2.14 

113. 

Iv 1 ■ 47 

7 

73366. 

ESMOND 

PROVID£NC:E 

BRATTLE3ARE 

WINDHAM 

2.12 

‘93. 

^ 

S 

6353?. 

ESMOND 

PROVIDENCE 

KEENE 

CHESHIRE 

2. 09 

:35. 

112. 

259. ’-'0 

■9 

32193. 

HOULTON 

AROOSTOOK 

BANGOR 

PBC’BSCOT 

2.03 

o4 « St. 

10 

83935. 

W'/OMINO 

WASHINGTON 

WINSTEAD 

LITCHFIELD 

2.01 

99. 

500.00 

11 

53401. 

BOSTON 

SUFFOLK  PSD 

PORTLAND 

CUMBERLAND 

1.91 

105. 

500.00 

12 

24126. 

BOSTON 

SLFFO-K  PSD 

LACCWIA 

BELKNAP 

i  .  V-T 

96. 

W..V  • 

13 

20133. 

ESMOND 

PROVIDENCE 

GREENFIELD 

FRAIVKIIN 

i.6'1 

C’li 

14 

129t>37. 

ESMOND 

PROVIDENCE 

NORTHHAMPTON 

HAMPSHIRE 

1.53 

72. 

500. CO 

15 

5569:3. 

NASHUA 

HILLSB0R01X3H 

KEENE 

iSHESHlRE 

1.42 

42. 

259.90 

16 

103092. 

BijSTON 

SUFFOLK  PSD 

ro-:hester 

STRAFFORD 

1.40 

77. 

500.00 

17 

131612. 

N^lSHJA 

HILLSBOROUGH 

PLYMOUTH 

DRfffTON 

1.35 

l>i. 

114.93 

500.(90 

13 

214563. 

M3RWICH 

NEW  LONDON 

WINSTEAD 

LITCHFIELD 

1.'31 

61. 

1? 

61642. 

NA:S'rtJA 

BOSTON 

hillsboroixh 

LACONIA 

BELKNAP 

1 . 23 

'I'L. 

■jJj* 

20 

196604. 

SUFFOLK  PSD 

CONOjRD 

MERRIMACK 

w  » 

,y2' . 

ijUo.  /  £’ 

21 

22 

93177. 

EGMCND 

PROVIDENCE 

BGLTCN 

TOLLAND 

1.  i? 

47. 

500. CO 

125699. 

SANFORD 

YORK 

LEWISTON 

ANDROSCOGGIN 

1.15 

53. 

500.00 

22 

290220. 

SPRINGFIELD 

HOTDEN 

PITTSFIELD 

BERKSHIRE 

I.IO 

56. 

462. 13 

24 

25 

59567. 

NEWP'GRT 

NEWPORT 

DANIELSON 

WINDHAM 

1.05 

43. 

500.  M 

171310. 

PORT'-SMCUTH 

ROCKINGHAM 

PORTLAND 

CUMBERLAND 

.36 

47. 

500.00 

26 

103496. 

:SPRINGFIELD 

HAMPDEN 

ORENFIELI' 

FRANKLIN 

34. 

267.25 

27 

45169. 

BURLINGTON 

CHITTENDEN 

MONTPELIER 

WASHINGTON 

khl 

34. 

lc7 

23 

106121. 

ESMOND 

PROVIDENCE 

DANIELSON 

WINDHAM 

• 

X  Ji 

500.00 

i  1...  -I'S 

29 

53812. 

BURLINGTON 

CHITTENIEN 

VERGENNES 

ADDISCN 

.55 
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TE  C-HApACTERISTICS  AND  LOADINGS 

~wf.y 


“C 

FAIR 

EvAOl'cES 

city 

COUNTY 

CITY 

'.I.LihT1 

CTST^'.CE 

TSA'.-'L 

CZ.AY 

•  .  •  .”w 

1 

1 

226S6. 

ESMOND 

PROVIDENCE 

BBYNINOTON 

BENNINGTON 

12'?.7 

2,;o 

^■.  4i 

C  * 

57590. 

BOSTON 

SUFFiXK  R3D 

BATH 

SAGADAHOC 

132.4 

2/^4 

4.62 

^1 

759. 

NEInFORT 

NEWPORT 

WilLSTEAD 

LT'CHFIE'lD 

:C6.6 

2.77 

4.43 

1  * 

00 

9 

I'VlQ. 

NEInFCRT 

r.EWPORT 

WINSTEAD 

LITCHFIEuD 

133. 9 

3.05 

4.12 

“7  ' 

4 

4 

39003. 

BOSTON 

SUFFOLK  PSD 

WILDER 

WINDSOR 

124.6 

y  “fj 

6, 12 

C-  ( 

'•'S 

4 

736. 

BOSTON 

a'FFCLK  PSD 

WILDER 

WINDSOR 

127.2 

2!?4 

5.99 

7~2 

42 

4 

62316. 

BOSTON 

SLiFFOLK.  PSD 

WILDER 

'.'iINDSuR 

133.0 

•: .  •-■2 

c  C7 

S  ■ 

c 

5 

69160. 

BOSTON 

SUFOLK  PSD 

CLAREMONT 

SULLIVAN 

111.7 

2.75 

6!i2 

3.70 

:S 

5 

2966. 

BOSTON 

SUFFOLK  PSD 

CLAREMONT 

SULLIVAN 

114.3 

2.75 

=;,99 

3.74 

0 

L 

73366. 

ESMOND 

PROVIDENtCE 

BRATTLEBARE 

WINDHPH 

98. 1 

*.»  -J'J 

5.21 

7t71 

9 

7 

63539. 

ESIMOND 

PROVIDENiCE 

KEENE 

CHESHIRE 

34.4 

«  '•j'y 

-.'i. 

4.70 

7.02 

7 

3 

32193. 

HOULTON 

AROOSTOOK 

BAdNGOR 

PeC'BSCOT 

111. 5 

•*v  .*0 

L»  ‘♦u 

4.00 

6. 43 

V 

Q 

4563. 

WVOMINij 

WASttINGTGN 

WINSTEAD 

LITCHFIELD 

90.7 

••  '•ir. 

to. 

4.43 

6.31 

31 

9 

79417. 

WYOMIMj 

WASHINGTON 

WINSTEAD 

LITCHFIELD 

113.0 

2.66 

4.12 

’4 

10 

10 

67C44. 

BOSTON 

SUFFOLK  PSD 

PORTLAND 

CUMBERLAND 

103.0 

to  .  'to*7 

4.. 62 

6.91 

n 

11 

23130. 

BOSTON 

SUFOLK  PSD 

LAiSijNIA 

BELKNAP 

94.7 

2.20 

6.12 

P.3’2 

'7 

11 

'?^6. 

BOSTON 

SUFFOLK  PSD 

LAC-iONIA 

3EU(NAP 

97.3 

2.37 

5. '99 

h'Sb 

i  •. 

12 

20133. 

ESMOND 

PPCVIDENCE 

C-REEIFIELD 

FRANKLIN 

:30.3 

"T'  1  * 

7,  ^2 

1  ■' 

13 

115722. 

ESMOND 

PROVIDEN.:^ 

NORTFiAMF'TON 

16^MF"SHIRE 

71,7 

liss 

5!l6 

7. '32 

4:i 

13 

13965. 

ESMOND 

PROVIDENCE 

NORTHHAMPTON 

HAMPSHIRE 

119.0 

2.77 

4.31 

7.03 

:o 

14 

55693. 

NASHUA 

HlLLSBORiOUSH 

KEENE 

CHESHIRE 

42.5 

1.42 

4.00 

5.42 

14 

15 

103092. 

BOSTON 

SUFOLK  PSD 

ROCHESTER 

STRAFORD 

75.6 

74.2 

1.63 

4.00 

5.. 63 

16 

43844. 

NASHUA 

HILLSBtSROUGH 

PLYMOUTH 

GRAFTON 

1.65 

6.12 

7.77 

“lO 

.'V 

16 

82763. 

NASHUA 

HILLSBOROUGH 

PLYMOUTH 

GRAFTON 

74.2 

1.76 

5.99 

7.75 

‘  7 

17 

93425. 

NCfiWiCH 

NEW  LGNDC?) 

'(*  instead 

LITCHFIELD 

1.55 

4.43 

5.97 

17 

121143. 

NORWICH 

NEW  LONDON 

WINSTEAD 

LITCHFIELD 

S3. 3 

i.:-6 

4.  CO 

* ,  'r  6 

lo 

IS 

13644. 

PORTSMOUTH 

ROaCMHAM 

LEWISTON 

ANDROSvDGGIN 

72.0 

1.60 

4.62 

6i  12 

13 

19 

22673. 

NASHUA 

HILLSBORiiei 

LACONIA 

BELKNAP 

61.3 

1.47 

6.12 

7.59 

34 

19 

38765. 

NASHUA 

HILLSBORiXHjH 

LACONIA 

BELKNAP 

61.3 

1.53 

5. '99 

7.57 

43 

19 

2iX). 

NASHUA 

HILLSBiSROLijH 

LACONIA 

BELKNAP 

79.5 

2.41 

5.57 

l.n 

19 

20 

13!3:393. 

BOSTON 

SUFFOLK  PSD 

CONCORD 

MERRIMACK 

■65.6 

1.46 

6.12 

7.53 

39 

20 

3206. 

BOSTON 

■SUFFOLK  PSD 

C0ND3RD 

MFRIMACK 

63. 2 

1..63 

5.99 

7. 62 

“  < 

21 

93177. 

ESMOND 

FRCVIDENiCE 

BOLTON 

TOLLAND 

47.5 

1.36 

4.16 

5)52 

'“■71 

22 

119039. 

SANFORD 

YORK 

LEWISTON 

ANIiRDSCOGC'IN 

52.9 

1 .  S5 

5.34 

6.6r' 

33 

22 

6660. 

SANFORD 

YC«K 

LEWISTON 

ANDROSCOGGIN 

7.S.4 

2.33 

4.35 

6. 3‘9 

nn 

23 

129561. 

SPRINGFIELD 

HAMPDEN 

PITTSFIELD 

BFKSHIRE 

56,0 

1.32 

=:  -’1 

7.03 

40 

77 

126503. 

SPRINGFIELD 

HAMPDEN 

PIFFIELD 

BERKSHIRE 

54.4 

1.43 

Vbi 

7.10 

46 

23 

34156. 

SPRINGFIELD 

'HAMPDEN 

PITTSFIELD 

2'ERKSHiRE 

’7.6 

■> 

4. 46 

24 

24 

59402. 

NEUPORT 

NEWPORT 

DANIELSON 

WINDHAM 

4;3.2 

• 

i.to> 

5.51 

1  7=; 

47 

24 

165. 

NEWPORT 

NEWPORT 

DANIELSON 

WINDHAM 

91.2 

2.51 

4.39 

'  i“.  '■ 

25 

25 

157667. 

PORTSMOUTH 

ROGOINOHAM 

PORTcAND 

CUMBF'LAND 

46.3 

1.04 

4.62 

5.66 

41 

26 

103496. 

SPRINGFIELD 

HAMPDEN 

GPEENFIELD 

FRANKLIN 

34.1 

.76 

4.46 

C  “•••> 

„/ 

27 

45169. 

B.Rl-lNOTON 

CHITTENDEN 

MONTPELIER 

'rtAcHINbiON 

■4.  '■ 

TL 

4.CC 

4. . 

23 

23 

64284. 

ESMOND 

PROVIDENCE 

DANIELSON 

'WINDHAM 

24.7 

!67 

5.51 

6.19 

23 

41S37. 

ESMOND 

PROVIDENCE 

DANIELSON 

WINDHAM 

;6.0 

l.SO 

4.3'’ 

C'  •  •  r 

29 

29 

53312. 

BURLINGTON 

iGilTTENDEN 

VFuENNES 

ADBISSN 

.14, 1 

4.30 
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I’ABLE  3. ! 


EVADJATION  DEF'lfiTliFtE  3UMKARV 

EyACUATICN  COUNTY  EVfiDJEES  PEFlCEf'iTAOE  LEFARTED  AT  ELAPSED  HOUR 

:■  A  L  .}  C  ••'A  -J  AA 


PROVIDENCE 

514214. 

2 

2'0 

4 

41 

L 

‘j 

61 

0 

0*' 

12 

99 

16 

9'? 

20 

99 

24 

Oil 

y? 

99 

7 

60 

72 

99 

S4 

99 

SUFFOLK  PSD 

S2264I. 

19 

S3 

57 

7r 

•  w 

99 

9? 

99 

99 

99 

‘’*‘9 

9'9 

99 

'j'i 

00 

NEWPORT 

73245. 

19 

S3 

57 

77 

99 

9? 

99 

99 

99 

99 

99 

99 

99 

^9 

AROOSTOOK 

32153. 

24 

4? 

74 

99 

99 

9? 

7 

■;9 

99 

99 

99 

?9 

WASHINGTON 

3.?-'85. 

24 

4:3 

% 

99 

99 

9? 

?9 

;9 

99 

99 

99 

99 

HILLSBOROUGH 

24:S"?47. 

13 

Si 

55 

73 

?? 

99 

99 

99 

99 

9  9 

9v 

99 

ROCKINGHAM 

171310. 

21 

43 

64 

36 

99 

99 

?9 

90 

99 

09 

OQ 

9'*' 

v'J 

9'? 

NE'W  LONDON 

21454B. 

22 

47 

71 

"^5 

99 

99 

99 

99 

"'C 

7 

'^9 

9v 

9? 

99 

YORK 

12Si'^'4. 

13 

37 

56 

99 

9? 

99 

99 

99 

't'O 

99 

99 

Dl’C 

9? 

HAMPDEN 

Sw371S. 

19 

S3 

77 

99 

99 

99 

99 

99 

99 

99 

v9 

99 

9‘'^ 

CHITTENDEN 

lOGSei. 

24 

4Q 

74 

99 

99 

99 

99 

99 

99 

go 

99 

99 

99 

99 

ITABLE  3.2 

RECEPTION  ARRIVAL  SUnWY 

RECEP'ION  CCLTJTY  EVAaflS  PERCENTAGE  ARRIVED  AT  ELAPSED  HOUR 


4 

6 

8 

12 

16 

20 

24 

32 

40 

48 

T’’ 

84 

BOJNINGTQN 

22636. 

0 

ej 

24 

44 

32 

99 

9? 

99 

V^' 

99 

99 

99 

9? 

QO 

SAGADAHOC, 

57540. 

) 

li 

33 

54 

97 

09 

99 

99 

95 

99 

99 

c..j 

LITChr IElD 

312231. 

47 

'i 

00 

’yV 

99 

•‘■V 

■  - 

'IQ 

9'- 

'•  7 

WiNCo'.R 

lOISMO. 

7 

24 

42 

'  V 

V 

99 

0.5 

T.b 

99 

'll 

SULLIVAN 

72126. 

f) 

li 

*7 

44 

76 

99 

99 

99 

■ - 

10 

'v’i 

-f 

r'r 

WINDHAM 

73866, 

0 

14 

33 

ff-i 

91 

99 

99 

99 

og 

CQ 

99 

99 

99 

99 

PENOBSCOT 

S2r<3. 

0 

19 

44 

69 

99 

99 

■99 

99 

gg 

99 

99 

■99 

CO 

99 

CHESHIRE 

124232. 

3 

24 

47 

70 

99 

99 

9? 

99 

99 

99 

99 

•5'5 

99 

99 

CUMBERLAND 

224711. 

7 

77 

49 

71 

69 

99 

■?9 

99 

■9? 

9'? 

9? 

99 

BELKNAP 

3576;3. 

•*> 

13 

35 

51 

34 

99 

99 

99 

?9 

OP 

QO 

cc 

C9 

FRANKLIN 

123634. 

ll 

33 

ctr 

-j'j 

/ 

99 

09 

99 

99 

■hb 

'9 

oo 

99 

00 

99 

HAMPSHIRE 

129.637. 

i 

2-0 

n 

59 

O’O 

99 

99 

99 

9? 

99 

99 

99 

'99 

99 

STRAFFORD 

103092. 

4 

29 

54 

79 

99 

99 

99 

99 

'99 

99 

>9 

99 

99 

.9.9 

GRAFTON 

131612. 

:• 

13 

52 

85 

99 

99 

99 

99 

^9 

90 

99 

!>ri 

9'9 

ANDROSCOGGIN 

139343. 

5 

24 

43 

63 

9'9 

'99 

•99 

99 

OQ 

99 

v9 

99 

99 

.99 

MERRIMACK 

196604. 

4 

20 

37 

53 

36 

■99 

■?9 

99 

99 

99 

99 

>p 

9? 

■99 

TOLLAND 

93177. 

7 

55 

79 

99 

9-? 

99 

99 

99 

99 

99 

•00 

99 

BERKSHIRE 

290220, 

4 

T-T 

40 

53 

93 

9-? 

99 

99 

■=>9 

99 

'99 

99 

99 

WINDHAM 

165633. 

7 

27 

4.6 

65 

99 

99 

99 

99 

99 

99 

OQ 

00 

"9 

O-D' 

WASHINGTON 

45169. 

15 

40 

t5 

60 

99 

99 

9? 

99 

00 

-5 

Q-y 

ob 

OQ 

'D'O 

ADDISON 

c--.* 

17 

42 

9< 

99 

00 

•7'.’.' 

■1 0 

7’’^ 

99 

99 

DO 

:  - 

52 


1  table  B.3 

';rWjRK  LINK  i>IftRfti:TERI:3TIC3  AND  LGAiaNCtS 


■RON 

TO 

ROUTE  NAME 

DIST. 

AVON 

UINSTEAD 

U.S.  ROUTE 

44 

20.2 

ASHFORD 

BOLTON 

U.'S.  ROUTE  54A 

13.3 

IlUTH  foster 

DANIESiON 

U.S.  ROUTE 

6 

10.3 

'L'TNhM 

ASHFORD 

U.S.  ROUTE 

M 

12.5 

:hef'Achet 

PUTNAM 

U.S.  ROUF 

44 

is.o 

rARFORD 

AViON 

U.S.  ROUTE 

44 

9.5 

l-ERIDAN 

MILLDALE 

iiTATE  ROAD 

66 

5.4 

,-ORTLAND 

EWISTON 

STATE  ROAD  900 

25.2 

BIDDLEFORD 

PORTLAND 

INFR3TATE 

95 

12.3 

SANFORD 

PORTLAND 

U.S.  ROUTE 

202 

27.7 

BICL’LEORD 

INTESTATE 

■95 

34.0 

-ZZ 

PITTSFIELD 

U.S.  ROUTE 

7 

12.4 

NORTHHAMPTON 

PITTSFIELD 

STAF  ROAD 

9 

35.4 

NOR'^HHAnPTON 

GREENFIELD 

INTERSTATE 

?l 

15.0 

ORANGE 

GEENFIELD 

3TATC  ROAD 

2 

15.0 

OHICOPE 

NORTHHAMPTON 

U.S.  R0JTE 

5 

14.0 

BELOHERTCUN 

NOFTHHAMPTON 

STATE  ROAD 

9 

11.4 

i.ATERSVILLE 

IwjRCHESTcR 

STATE  ROAD  146 

18.3 

V.ESFIELD 

EE 

INTERSTATE 

90 

"/! 

BOSTON 

LAWENCE 

INTESTATE 

93 

25!  3 

UORCHESTER 

BELCHERTOHN 

3TAF  ROAD 

9 

29,3 

WCRCHESTER 

ORANGE 

STAF  ROAD  122 

34.7 

OHICOPE 

UE3FIELD 

INTESTATE 

90 

7.3 

NCRTHFIELD 

PLYMOUTH 

INTERSTATE 

93 

23.0 

•‘Af. CHESTER 

CONCORD 

INTER3TAF 

93 

17.5 

'•ANCHESTER 

CONCORD 

U.S.  ROUTE 

3 

17.5 

^AHRENCE 

MANCIESFR 

INTERSTATE 

93 

22.3 

NASHUA 

MArCHESTE 

U.S.  ROUTE 

3 

15.1 

ESMCND 

ASTON 

INTESTATE 

2'?5 

4.5 

.GHNSTON 

SOUTH  FOSTER 

U.'i.  ROUTE 

6 

8,5 

HVOMING 

SOIJTH  FOSTE 

STATE  ROAD  102 

24.0 

ESMOND 

CHEACHET 

U.S.  ROUTE 

44 

•■V  c 

ASTON 

SLAFRSvIEE 

STATE  ROAD  146 

7.3 

FEED 

CAP.  FERCBaAGE 

Ah  A' 

•  T'., 

lE 

- 

2  4  6 

"j 

J 

lA 

24 

4; 

1- 

35.0 

2300.  96  liX)  ICO 

IGO 

) 

0 

■‘\ 

) 

f 

35.0 

2300.  lOi’i  100  100 

1)0 

s 

f. 

'J 

35.0 

2300.  100  100  100 

ICO 

41 

6 

') 

0 

•J 

25.0 

2300.  lOO  100  ICO 

DO 

) 

0 

(: 

0 

35.0 

2300.  10>  100  100 

DO 

\J 

) 

0 

G 

35.0 

2b vO.  v9  luO 

ICC 

h 

V 

c 

0 

0 

35.0 

iSOu*  100  lUO 

ICO 

l‘; 

V 

f, 

45.0 

4000.  31  34  84 

34 

j 

0 

0 

r* 

0 

■j 

45,0 

3000.  53  ICk)  100 

IOC 

V 

I'i 

f! 

0 

Cl 

.'1 

) 

35.0 

2SC'i.  100  ICC  ICu 

DO 

\] 

6 

0 

V 

45.0 

3000.  ICO  ICO  ICO 

DO 

'i 

0 

■') 

1* 

35.0 

2300.  101  1)1  101 

101 

101 

6 

0 

0 

0 

0 

35.0 

2600,  100  100  100 

100 

76 

0 

c 

0 

0 

0 

0 

45.0 

4000.  100  100  100  100 

ii 

0 

0 

0 

0 

c 

35.0 

2300.  74  100  100 

100 

ICC 

n 

0 

0 

6 

0 

0 

') 

35.0 

1400.  68  97  97 

97 

0 

C 

0 

n 

0 

0 

0 

0 

35.0 

230^.  I'Xi  100  100 

rc 

T: 

0 

l'; 

■> 

45.0 

SOOO.  90  93  93 

O-'v 

0 

,*1 

r\ 

45.0 

4000,  70  70  70 

43 

6 

6 

0 

I': 

A 

6 

45.0 

3000.  32  32  82 

'Jim 

32 

0 

0 

c 

6 

0 

0 

0 

35.0 

2800.  100  100  100  100 

0 

0 

0 

0 

0 

0 

0 

0 

35.0 

2800.  100  100  UX)  100 

10 

0 

0 

0 

0 

0 

0 

0 

45.0 

4000.  70  70  70 

70 

0 

n 

0 

0 

!*) 

0 

0 

0 

45,0 

4000.  63  i3  68 

63 

63 

6 

0 

0 

0 

0 

0 

T' 

45.0 

3000.  IOC  100  100 

DO 

*.  0'*'; 

•'» 

c 

0 

A 

35,0 

2Sai.  ICO  11)0  100 

100 

100 

0 

c 

0 

0 

ri 

0 

45,0 

3000,  32  32  32 

32 

04 

0 

0 

0 

0 

f! 

0 

c 

35.0 

1400.  100  m  100 

100 

S'3 

0 

0 

0 

0 

6 

0 

r. 

45,0 

3000.  ■‘•'3  93  93 

93 

0 

0 

0 

0 

j 

{) 

0 

0 

35.0 

2800.  100  U'.Kj  100 

too 

c 

0 

0 

0 

6 

0 

0 

0 

30.0 

1200.  100  100  100  100 

0 

0 

0 

0 

0 

0 

0 

c, 

35.0 

2SiX).  100  100  100  100 

0 

0 

0 

0 

0 

0 

0 

0 

45.0 

3000,  93  93  93 

93 

0 

0 

0 

0 

0 

0 

0 

0 
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3.3  Route  Listing 

Each  route  is  displayed  in  detail  in  Table  2,3.  The  routes  appear 

in  the  same  order  in  which  they  are  summarized  in  Table  2.2.  Since 

Table  2.3  is  long  and  quite  detailed,  only  the  part  of  it  corresponding 

to  the  first  nine  routes  of  Table  2.2  has  been  reproduced  below. 

The  sequence  of  nodes  proceeding  from  the  origin  to  the  destination 

is  shown  for  each  route.  For  each  node  (or  link  into  the  node)  the 

following  information  is  given: 

NODE  -  the  5-dlglt  FIPS  code  of  the  state/county  with  the 
node's  single-digit  sequence  number  added  as  the  sixth 
digit. 

CODE  -  a  code  number  embedding  the  highway  segment’s 
characteristics  and  number. 

TRAVEL  -  elapsed  travel  time  to  the  node  at  evacuation 
speeds . 

DELAY  -  elapsed  delay  time,  i.e.,  the  longest  average 
delay  time  imposed  by  a  congested  link  already  listed. 

TOTAL  -  TRAVEL  plus  DELAY  time. 

DIST  -  accumulated  distance  in  miles. 

SPEED  -  link  evacuation  speed  in  miles/hour. 

CAP  -  link  evacuation  capacity  in  vehicles/hour  for  all 
available  lanes . 

VCLS  -  number  of  vehicles  using  the  link  according  to  the 
route  loading. 

PLACE  NAME  -  name  of  county  center. 

COUNTY  -  county  name . 

ROUTE  NAME  -  highway  type  and  number. 
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ITflBLE  2.3 
ROUTE  I 


fiiODE 

core 

TRAVEL 

DELAY 

TOTAL 

BIST, 

SPEED 

CAP. 

440070 

12295 

.00 

4.00 

4.00 

.0 

,0 

0. 

440071 

:122^5 

.10 

4.70 

4.30 

4.5 

45.0 

3000. 

440077 

33146 

.27 

4.70 

4.97 

12.3 

45.0 

3C00. 

250270 

33146 

.69 

4.70 

5.39 

31.1 

45.0 

SOiXi. 

250112 

35122 

i.6;3 

5.21 

6.39 

65.3 

:35.0 

2800. 

250110 

35002 

2.1L 

5.21 

7.32 

30.3 

35,0 

2300. 

250111 

12091 

2.25 

5.21 

7.46 

87.1 

45.0 

4i»0. 

500250 

12091 

2.50 

5.21 

7.71 

93.1 

45.0 

4000. 

500030 

35009 

3.40 

5.21 

9.61 

129.7 

35.0 

14tXi. 

ROUTE 

NODE 

DUDE 

TRAVEL 

DELAY 

TOTAL 

DIST, 

SPEED 

CAP. 

250250 

12095 

.00 

4.00 

4.00 

.0 

.0 

0. 

2500‘^0 

12095 

.34 

4.00 

4.34 

15.5 

45.0 

3000. 

2500'?3 

12095 

.7? 

4.00 

4,79 

35.5 

45.0 

3000. 

330151 

12095 

1.06 

4.00 

5.06 

47.6 

45.0 

SMO. 

330150 

:i2C95 

1.25 

4.00 

5.25 

56.1 

45.0 

300<). 

230311 

12095 

2.00 

4.62 

6.63 

90.1 

45.0 

StX'O. 

230050 

12095 

■*) 

4.62 

6.91 

103,0 

45.0 

3000. 

230051 

12C95 

2., 31 

4.62 

7,43 

126.4 

45.0 

41X0. 

^vJ2oO 

23001 

2,94 

4.62 

7.57 

132.4 

45.0 

40C<'. 

ROUTE 

NODE 

CODE 

TRAVEL 

rajiY 

TOTAL 

0I3T. 

SPEED 

CAP. 

440OE0 

35138 

.lYi 

4.00 

4.00 

.0 

.0 

0. 

440.;, 42 

w  w  A  w  w 

.17 

4 . 00 

4.17 

6.  1 

35.0 

2300. 

44';091 

23001 

,63 

4.60 

4.63 

26.7 

45.0 

4«i0. 

ROilll 

12095 

.63 

4,00 

4.63 

23.9 

45.0 

4000. 

90110 

90113 

36002 

1.22 

4,00 

5.22 

45.2 

30.0 

1200. 

33002 

1.52 

4.00 

5.52 

53.6 

45.0 

3000. 

<>0037 

33002 

1.63 

4.00 

5.63 

6.3.6 

45.0 

3000. 

33002 

1.95 

4.00 

5,95 

77.8 

45.0 

3000. 

>0035 

25044 

*.«  i  >' 

4.43 

6.6i 

36.3 

35.0 

2300. 

90050 

25044 

2.77 

4.43 

7.20 

10.6.6 

35.0 

2300. 

ROUTE 

30 

NODE 

CODE 

TRAVEL 

lELAY 

TOTAL 

DIET. 

SPEED 

CAP. 

440050 

35133 

.00 

4,00 

4.00 

,0 

,0 

0. 

4400<>2 

35133 

.17 

4.00 

4.17 

6.1 

35.0 

2800. 

ii. 

''Tt**.*'  * 

.  2' t 

4.00 

4.63 

26..7 

45.0 

4000. 

"0 1 1 1 

12075 

t 

4.00 

4.63 

2"3,9 

45.0 

4000. 

>0112 

12075 

.75 

4.00 

4.95 

41.0 

45.0 

40>X. 

-'''i'll') 

1  ---'iC 

4.00 

5.32 

57.7 

i5.0 

4000. 

•■x'l 

33007 

1.70 

4,00 

5.70 

74.3 

45.0 

mo. 

>C034 

4.00 

5.91 

84.0 

45,0 

4ix«‘i. 

ZO'nO'j 

j 

11071 

2. 10 

4.  CO 

6.10 

92.6 

45.0 

6000. 

70036 

35066 

2.25 

4.12 

6.37 

98.0 

35.0 

2800. 

^0091 

12034 

2.52 

4.12 

6.64 

210.0 

45,0 

3000. 

>0052 

33003 

2.59 

4.12 

6.71 

113.2 

45.0 

3000. 

>0054 

33003 

2.66 

4.12 

6.98 

125.3 

45.0 

3000. 

90050 

33003 

3.05 

4.12 

7.17 

45.0 

3C0i). 

VCLS 

PLACE  NAME 

ccLirv 

FCC’E  NAME 

0. 

ESMOND 

PROVIDENiCE 

INTERSTATE  295 

75224. 

ASTON 

PRiDVIDENiCE 

INTERSTATE  295 

75224. 

3LATERSVILLE 

FPOvirEi'.CE 

•dTATE  .^lAD  1*^6 

75224. 

UORCHESTER 

WORCESTER 

oTATE  f^uAD  14& 

29173. 

ORANGE 

FRANKLIN 

STATE  ROAD  122 

29173. 

GREENFIELD 

FRANKLIN 

STATE  ROAD  2 

24133. 

BERNARDSTON 

FRANKLIN 

INTERSTATE  91 

24133. 

BRATTLEBARE 

wINDHAM 

INTERSTATE  91 

5672. 

BENNINGTON 

BENNli'i^.'TUN 

STATE  RiJAD  9 

vas 

PLACE  NAME 

CDUMTY 

ROUTE  ^JAME 

0. 

BOSTON 

SUFFOLK  PSD 

INTERSTAlt 

C'C; 

56931. 

DANVERS 

EiiU 

INTERSTATE 

56931. 

AMESBURY 

ESSEX 

INTERSTATE 

95 

56931. 

HAMPTON 

ROCKINGHAM 

INTERSTATE 

95 

56931. 

PCRTSMOIJTH 

ROCKINGHAM 

INTERSTATE 

95 

739;56. 

BIDDLEFOfiD 

YORK 

INTERSTATE 

■>5 

73936. 

PORTLAND 

ClMBERLAND 

INTERSTATE 

•36 

14397. 

BRUNSWICK 

'XMEtRLAND 

INTERSTATE 

95 

14397. 

BATH 

iAGADiAHOC 

U.5.  ROUTE 

1 

VIXS 

PLACE  NAME 

COUNTY 

PruiTE  NAME 

0. 

NEWPORT 

NEWPORT 

STATE  ROAD 

<  ^ 

IS270. 

SAUNDERSTOWN 

WASHINGTON 

STATE  ROAD 

135 

3419. 

WESTERLY 

WASHINGTON 

U.S.  ROUTE 

1 

24416. 

PAWCATUCK 

NEW  LCNDiON 

INTERSTATE 

95 

1337. 

NORWICH 

NEW  LONDON 

STATE  RiSAD 

2 

54979. 

COLCHESTER 

NEW  LONDON 

STATE  ROAD 

2 

54979. 

MARLBORO 

HARTFORD 

STATE  ROAD 

54979. 

HARTFORD 

HARTFORD 

STATE  ROAD 

24737. 

AVON 

HARTFORD 

U.S.  ROUTE 

44 

24787. 

WINSTEAD 

LITCHFIELD 

U.S.  ROUTE 

44 

VCLS 

PLACE  NAME 

COUNTY 

RiXiTE  rVAME 

0. 

NEWPORT 

NEWPORT 

STATE  ROAD  138 

18270. 

SAIJNCERSTOWN 

WASHUXmCN 

STATE  R‘.AD  13B 

3419. 

WESTERL! 

,.1.  PC'.’E  . 

24416. 

PAWCATUCK 

NEW  ..CNICN 

INTERSThTE  95 

23X)7>. 

NEW  LONDON 

NEW  LCNDON 

INTERSTATE  >5 

23079. 

OLD  SAYERCCK 

mi::lese( 

IN.ERS’A^E  <5 

/  9  • 

MIDDLETOWN 

MIDDLESEX 

:  ’A'E  *' 

23079. 

EAST  BERLIN 

H9RT^CRD 

;TATE  >CAD  0 

MERIDAN 

NE-  nAVEN 

•* 

23079, 

MILLDALE 

mARTlRu 

rATE  ROAD  66 

53272, 

WATEKBURY 

■.'Si  HAVEN 

INTERSTA~  S4 

53271. 

THCMASTOWN 

LI'CHFiELD 

STATE  ROAD  S 

53271. 

TCfiRlNOTON 

lITDHFIE-D 

I'A^E  ; 

53271. 

WINSTEAD 

tlX'r.rlELCi 

;  1  H  ^  fc.  r  .  nL  ; 
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RDJTE  4 


CODE 

IhMll 

DELAY 

TOTAL 

Dm. 

SPEED 

CAP. 

VCLS 

PLACE  NAME 

3:50250 

T  * 'O': 

4.00 

4.00 

.0 

Oi 

.2  - '  L  I'i 

*  ^  t  - 1 

j 

C  =.4 

-'C  “• 

45.  f 

7””^1 « 

•  -.'J 

sioo 

o.U7 

4I!! 

45.  •} 

-scco! 

79^51. 

-■•-•,•  '“■c 

ioos’ 

12039 

t  '  :o 

't'f. 

%.  * 

•  4  K. 

6. 12 

A. 

7!  73 
S.24 
3.71 

73.6 

95.6 

45.0 

45.0 

ii.-iX), 

iGOO* 

43546.! 

r,r:-iINTLt.N 
NEW  LCNDON 

Zl  ii. 

0  ‘^••’5 

■Jt  4^ 

1.12 

6.12 

•  it.  t 

45.0 

4000. 

2'5515. 

LEBANON 

500070 

12039 

2.17 

3.89 

12V.  6 

45.0 

40t)i0. 

25515. 

WILDER 

“iRR'J'rCK 
'■•rR'iy:ACK 
^'.R-iFTvN 
14 1 


RCiJTE  35 


NCDE 

CODE 

TRAVEL 

2:0250 

23C03 

.00 

250171 

23003 

.62 

3301:10 

22003 

.79 

S30111 

32900 

1.13 

330130 

,5003 

1.63 

330133 

12099 

1.31 

3301S5 

12039 

2.29 

33C091 

120S9 

2.76 

5C0270 

12039 

2.94 

DELAY 

TOTAL 

DIST. 

4.00 

4.00 

.0 

4.00 

4.62 

27.. 3 

4.00 

4.79 

35.6 

4.00 

5.13 

50.7 

5.99 

7.62 

63.2 

5.99 

7.80 

76.3 

5.9'9 

S.2S 

93.2 

5.99 

3.75 

119.2 

5.99 

3.93 

127.2 

SPKD 

CAP. 

vcls 

.0 

0. 

0, 

45.0 

4001). 

1S77S. 

45.0 

4fiO0. 

18773. 

45.0 

3000. 

51502. 

35.0 

2300. 

33550. 

45.0 

4000. 

230t)7. 

45.0 

3COO. 

43546. 

45.0 

4000. 

25515. 

45.0 

4000. 

25515. 

PLACE  rWME 
2DSTDN 
LCWELL 
fiiASHUA 
flANCHESTER 
CONCCRD 
HCFK INTOWN 
•NEW  LlNCDN 
LEBANON 
WILDER 


COUNTY'  ___ 

OLEFCLr-  FiD 

illDDLEEEX 

HILLSECROLiiOH 

HILLE'BCEOUijH 

lYERRIWACK 

^€RRI?!ACK 

YERRII'ACK 

ORAFTDN 

WINISCiR 


U.3.  ROUTE  3 

jj  C 

iTATE  RCAD  900 
:j.S.  PCL.'TE  3 
'NTEFSTATE  'o9 
INTERiTATE  i9 
INTERSTATE  39 

CO 


ROUTE  42 


NODE 

CCiDE 

TRAVEL 

250250 

23003 

.00 

250172 

23003 

.62 

330110 

22003 

.79 

330111 

25003 

1.22 

'  •“  A  4 

'36114 

2.03 

320133 

23202 

2.19 

330135 

12039 

2.67 

330091 

12039 

3.14 

500270 

12089 

3.32 

ROUTE  5 


NODE 

CODE 

TRAVEL 

250250 

12093 

.00 

2500? 1 

12093 

.56 

230111 

12093 

1.07 

330 130 

12093 

1.46 

330133 

12039 

1.64 

230135 

1203'! 

2.12 

330191 

:^5Ci  1 

2.35 

330190 

35011 

2.53 

lELAY 

TOTAL 

DIST. 

4.00 

4.00 

.0 

4.00 

4.62 

27.3 

4.00 

4.79 

25.6 

5,57 

6.79 

50.7 

5.57 

7.60 

74.9 

5,57 

7.75 

C'7  Cl 

5.57 

8.24 

103.9 

5.57 

3.71 

125.0 

5.57 

3.89 

133.0 

DELAY 

TOTAL 

DIST. 

4.00 

4.00 

.0 

5.00 

5.56 

nc  0 

5.00 

6.07 

43.1 

6.12 

7.53 

65.6 

6,12 

7.76 

73.6 

6.12 

3.24 

95.6 

6.12 

3.47 

103.6 

6.12 

3.70 

111.7 

SPEED 

CAP',  vas 

.0 

0.  0. 

45.0 

4000.  18778, 

45.0 

41X90.  13778. 

35.0 

1400.  15589. 

30.0 

24'Xi.  ISS'S?. 

45.0 

4000.  15539. 

45.0 

3000.  43546. 

45.0 

4000.  25515. 

45.0 

4000.  '25515. 

SPEED 

CAP. 

VCLS 

.0 

0. 

0. 

45.0 

3000. 

79951. 

<15.0 

3000. 

79951. 

45.0 

3000. 

97903. 

45.0 

4000. 

28007. 

45.0 

30<». 

43546. 

35.0 

2800. 

13031. 

35.0 

2800. 

13031. 

PLACE  NAME 
BCSTON 
LOWELL 
NASHUA 
MANCHESTER 
HENNIKER 
HORKINTCWN 
NEW  LCNCCN 
LEBANON 
WILDER 


PLACE  NAME 

BOSTON 

LAWRENCE 

MANCHESTER 

CONCORD 

HOFKINTCWN 

NEW  LGNCGN 

NEWPORT 

CLAREMiDNT 


OHJNTY 

SUFFCLK  PSD 

MIDDLESEX 

HILLSBCiROUGH 

HILLSB0R0UI3H 

YERR I  MA'.X 

MERRIMACK 

MERRIMACK 

ORAFTON 

WINDSOR 


CiXNTY 
SUFFOLK  PSD 

CCCC)( 

HILLSBCROUOH 

MERRIMACK 

MERRIMACK 

MERRIMACK 

SULLIVAN 

SULLIVAN 


ROUTE  NAME 
U.3.  ROUTE  3 
U.S.  ROUTE  3 
U,S.  ROUTE  3 
U,S.  RCllTE  3 
STATE  ROAD  :14 
U.S.  ROUTE  :02 
INTERSTATE  S9 
INTERSTATE  S9 
INTERSTATE  S9 


ROUTE  NAME 
INTERS'ATE  -3 

li^ERSTATE  ’3 
interstate  S=' 
;vTcq-,-ATc  S'^ 

STATE  ROAD  11 
STATE  SCAD  r. 


ROUTE  36 


NODE 

CODE 

TRAVEL 

250250 

23003 

.00 

250171 

23003 

« 

030110 

22003 

.79 

350111 

32900 

1.13 

330130 

25003 

1,63 

330133 

12039 

1.31 

6301:35 

12039 

2.2? 

330191 

35011 

2, 5t 

330190 

35011 

2.75 

DELAY 

TOTAL 

DIST. 

4.00 

4.00 

.0 

4,00 

4.62 

27.3 

4.00 

4.79 

•5.6 

4.00 

5.13 

50.7 

5,99 

7.62 

63.2 

5.'?? 

7.30 

76.3 

5.99 

3.23 

98.2 

5.99 

3.51 

106.2 

5.'?9 

3.74 

114.3 

SPEED 

CAP. 

VCLS 

.0 

0. 

0. 

45.0 

4iXifi. 

13773. 

45.0 

4iX)0. 

13773. 

45.0 

3000. 

51502. 

35.0 

2300. 

33550. 

45.0 

4000. 

23007. 

45.0 

3000. 

43546. 

35.0 

2800. 

18031. 

35.0 

2300. 

18031. 

Place  name 

BOSTON 

LL'i^ELL 

NASHUA 

MANCHESTER 

CONCCRD 

HCPKINTDWN 

NEW  LiNEcN 

NEWPORT 

CuAREMCNT 


OXNTY 

SUFFOLK  PSD 

MI'OLESEX 

nlLLSEOROLSH 

HIuLSBORlUuH 

MERRIMACK 

MERRIMACK 

MERRIMACK 

SULLIVAN 

SULlIvAN 


Y.IirE  NAME 
J.S.  PCUTE 


STATE  RCAD  =^'00 


IN'ERSTATE  S9 
INTERSTATE  S® 
STATE  RCAD  11 
S'ATE  RCmD  11 


4.1  Program/Storage  Size  Limits 

A  number  of  constraints  are  imposed  by  program  storage  limits  on 
the  size  of  problems  which  can  be  handled  by  the  model.  These  con¬ 
straints  are  shovm  in  the  following  table. 


Constraint  Item  Limit 


Input/Output  Channels  11-20  10 
States  per  Regional  Network  10 
Regional  Network  Nodes  500 
Regional  Network  Links  2000 
Evacuation  Counties  40 
Reception  Counties  500 
Evacuation/Reception  County  Pairs  200 
Routes  on  the  Route  List  500 
Links  per  Route  19 


These  constraints  limits  have  been  set  to  conserve  computer  memory 
without  excluding  most  applications.  If  necessary,  the  limits  can 
be  expanded  by  increasing  array  dimensions  and  making  other  minor 
changes  in  the  programs. 


56 


1 


4.2  Subprogram  Descriptions 

MAIN  -  A  dummy  main  program  which  only  calls  MASTER. 

MASTER  -  Operates  the  master  control  sequence. 

PASS  -  Operates  the  route  selection  and  loading  algorithm. 

ORDER  -  Reorders  the  link  list  for  a  regional  network. 

PRICE  -  Performs  the  first  stage  of  the  Ford-Fulkerson 
fastest  route  computation  for  a  given  destination. 

ROUTE  -  Finds  the  fastest  route  from  a  given  origin  to  a 
given  destination. 

COUNTY  -  Prints  a  table  of  county  population  centroids 
within  a  region. 

LOCK  -  Locates  a  rentimbered  node  from  its  FIPS  code  and 
sequence  number. 

STACK  -  Gets  the  nodes  of  a  regional  network  from  the 
Node  File. 

ROADS  -  Gets  the  links  of  a  regional  network  from  the 
Links  File. 

SAVER  -  Saves/retrieves  preprocessed  regional  networks. 

SHOW  -  Prints  summary  table  and  route  listing  on  completion 
of  the  route  selection  and  loading  algorithn. 

SORTR  -  Twig  sort  routine  called  by  ORDER. 

BUILD  -  Operates  the  assembly  of  a  regional  transportation 
network  from  the  data  base. 

SETUP  -  Operates  the  direct  input  or  computation  of  a  crisis 
relocation  plan. 

LOADL  -  Performs  an  evacuation  simulation  for  a  single 
network  link 

EVAC  -  Operates  the  simulation  of  an  evacuation  given  a  set 
of  route  loadings. 
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LOADER  -  Computes  traffic  on  all  links  given  the  number  of 
evacuees  using  each  route. 

BLOCK  -  Blocks  congested  links  at  Step  4  of  the  Route 
Selection  and  Loading  algorithm. 

IF SAME  -  Checks  a  new  route  to  see  if  it  duplicates  a 
route  on  the  route  list  in  Step  2. 

TIMER  -  Computes  average  travel,  delay  and  total  time  for 
a  route. 

HEADER  -  Prints  a  table  heading. 

TRNPLX  -  Main  calling  sequence  for  solving  the  linear  program 
transportation  problem. 


RESET 

XCK 

SWAP 

RESD 

ROW 

GET 

DELTA 


-  Subroutines  called  directly/ indirectly 
from  TRNPLX. 


4.3  Procedure  Files.  Elements 
BPT*EXECUTE . COMPILE 

Compiles  source  code  to  relocatable  elements. 
BPT*EXECUTE.  ASSEMBLE 

Assembles  relocatable  code  elements. 
BPT*EXECUTE.  RUN 

Executes  program. 
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BPT*EXECUTE*  COMPILE 
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JiiJG.T  OBcECT. 

?F[;P,FU  3PT»;i;lJRCE.CCI1KCN,BPT'fiulJRCE. 

,j£-F 

?F :  N,  ^iO  s?l*i OURCE .  .'•A  1 N. EB^ECT. I'AIN 
'•E'.F 


iF^N.SO  5FT»ECUP':E..;aEF.CB.jECr.aRBEft 

IfT N, NO  BPT*5QERCE .PRICE, OBJECT . PRICE 
JEOF 

JFTN.NO  BPT»SOURCE.RQUE,CS..IECT.RrjlJTE 
'5E0F 

3FTN,N0  3PT»SQL:RL;E.STftC.K,0B..£CT.iTACK. 
2ECF 

SFTN.NO  3PT»SCURCE.C0!JNTY,0BJECT.C0LiNTY 
?ECF 

N. NO  BPT*SOURCE. LOCK. OBJECT. LOCK 

0FTN,  NO  BPT»SOURCE.  R0AC6, OBXCT.ROADS 
?EOF 

JFTN, NO  BPT»SOURCE . SAVER , OBJECT. SAVER 
2E0F 

SFTN, NO  BPTtSOURCE. SHOW, OBJECT. SHOW 
2E0F 

gFTN, NO  BPTtSOURCE . SORTR, OBJECT. SORTR 
«EOF 

»FTN,NO  BPT»30URCE.8UILD,!jB,£CT.3U[LD 
^OF 

aFTN.NO  BPT»SOURC:E.SETUP.OBJECT.SETUP 
lEOF 

8FTN, NO  BPT»SOURCE. LOADL.OBJECT.LOADL 
©EOF 

■?FTN, NO  EPT»SCURCE. EVAC , OBJECT. EVAC 
'?£0F 

SFl N,  NO  BPT»'SOURCE. LOADER, OBJECT. LOADER 
aEOF 

8FTN.  NO  BPT*SOURCE. BLOCK, OBcJECT. BLOCK 
aEOF 

aFTN.NO  BPT*SOURCE. IFSAME.O&JECT. IFSAHE 
aEOF 

aFTN.NO  BPT*SCUR':£.Tl:iER,OBJECT.TIi>ER 
aEOF 

«FTN,  NO  8PT»S0LIRCE.  HEADER, OBJECT. I£ADER 
aEOF 

aFTN. NO  BPT»SOURCE.  TRNPLX,  OEljECT. THNPLX 
aEOF 

aFTN.NO  BPT»SOURCE.RESET.OBJECT.RESET 
aEOF 

aFTN.NO  3PT»£0IJRCE.  XCK. OBJECT.  XCK 
aEOF 

anN. NO  EPT»SCLRC£. SNAP- OBJECT. SWAP 
aEOF 

aFTN, NO  BPTtSOURCE. RESD. OBJECT. RESD 
aEOF 

aFTN, NO  BPT»SCURCE.  ROW,  jB-JECT.hOW 
aEOF 

aFTN.NO  6PT»SOORCE. r.ET,uB.jECT.OET 
■aEOF 

aFTN.NO  BPT*3CURCE. DELTA.OBJECT.DELTA 
•aEOF 

aFTN.NO  BPT»SOURCE.f1IN.OB.JECT.HN 
^CF 

aF'N.NC  9PT»S0URCE.NEW-CB.JECT.NEW 
aEOF 

aFTN.  NO  SPT»S01:RCE. PI VOT .OajECT.PIVOT 

«EOF 
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BPT*EXECUTE*  ASSEMBLE 
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, OBJECT. SMAIN 
IN  DBJEC'I.MAIN 
IN  .jBi-EL'T."m':JER 
IN  OBJECT. PASS 
IN  OBJECT.  [ PEE 'V 
IN  OEvECT. PRICE 
IN  ISoE'IT.ROJiE 
IN  IBwECT.i’^CK 
IN  C5..ECT.:riJNTV 
IN  I3..ECT...:CK 
IN  OBJECT. RCADP 
IN  0B.:En. CAVER 
IN  OBJECT. SriO»i 
IN  object. SORTR 
IN  OBJECT. BUILD 
IN  OBJECT.  SETUP 
IN  OBJECT.  LCADL 
IN  OBJECT. EVftC 
IN  OBJECT. LOADER 
IN  OBJECT. BLOCK 
IN  OBJECT. TRNPLX 
IN  OBJECT.RESET 
IN  OBJECT. XCK 
IN  0BJECT.3WAP 
IN  OBJECT.RESD 
IN  OBJECT. ROW 
IN  OBJECT. GET 
IN  OBJECT. DELTA 
IN  OBJECT. MIN 
IN  OBJECT. NEW 
IN  OBoECT. PIVOT 
IN  OBJECT. TIMER 
IN  OBJECT. IFSAME 
IN  OBJECT. HEADER 
SEOF 

BPT*EXECUTE*RUN 


SASij.A  BPT*IJSND. 
iJASG.A  BPTtUSLD, 
?ASO,A  BPT*MAP, 
0A3G-T  FILEll. 
aftSG>T  FILEli. 
SASG.T  FILE13. 
@A30.T  FILE14. 
i5A3G,T  FILE15. 
®ASG>T  FILE16. 
SASG.T  FILE17, 
aASG.T  FILE13. 
?A3G,T  FILE19. 
2A3G.T  FILE20. 
■?ASO,T  ROUTES. 
®U3E  3i3PT*U3ND. 
SUSE  9,PPT*USLD. 
3U';E  10,BFT»MAP. 
SUSE  ll.FILEl:!. 
?U3E  12,FILE12. 
®!J3E  13,F:LE13. 
•aUSE  14,FILE14. 
®U3E  15,rILE15. 
@uSE  I6,FILE16. 
■?USE  17,FILE17. 
®U3E  18.FILE13. 
aUSE  19,FILE19. 
aUSE  20,FILE20. 
5USE  7, ROUTES. 
eXOT  OBJECT. SMAIN 


4.4  Source  Code  Listings 


COMMON  Block  Insert 
Main  Program 


Subprograms : 


Transportation  Subprograms 


MASTER 

PASS 

ORDER 

PRICE 

ROUTE 

COUNTY 

STACK 

LOCKS 

ROADS 

SAVER 

SHOW 

SORTR 

BUILD 

SETUP 

LOADL 

EVAC 

LOADER 

BLOCK 

IFSAME 

TIMER 

HEADER 


TRNPLX 

RESET 

XCK 

SWAP 

RESD 

ROW 

GET 

DELTA 

MIN 

NEW 

PIVOT 
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INSERT  PROC 

LOGICAL  OMIT 

DIMENSION  IV(l),IA(l),IS3(ll 
COMMi'jN  /Bl/IJR(21.500),JT0(2000).CR<500).DR(5<X)t 
COMMiDN  /  B2/LFR 1 2000 ) ,  I FR  ( 20tXi ) ,  MFR  ( 2000 ) .  CFR  ( 2000 ) ,  BFR  ( 2000 ) .  DFR  ( 2 
1000) 


lOMMON  /B3/N'  DES  ( 3 . 500 ) .  I  OPR  ( 2 .  IiXT ) .  lENPR  ( 200 ) .  JEV  ( 100 )  >  OMIT  ( 20(Xi 
l,2),SNAMES(4,60),«AP(4,60) 

COMMON  /B4'NCL0C(2.6O),ANA«E3(7,500),XLL0C:(2.5iXi),LL0CI5O0) 

COMMON  /B5.  /  i  1 500 , 40 ) ,  A 1 8000 ) .  8  i  2000 ) 

COMMON  /  F6  '  CSTLNl;  s  2000 ) .  SPDLNK  <  2m )  r  C APLNK !  2000 ) 

EQUIVALENCE  i lY- V ) i ( lAi A ) > i 188) 8) 


C 

C 


INCLUDE  SOURCE. INSERT 


CALL  MASTER 

STOP 

END 


SUBROUTIIC  HA:jTER 
C  MASTER  CONTROL  SEQUENCE 
INCLUDE  SOURCE. INSERT 
DIMENSION  NREGISO) , NREADY 13) 

DATA  K.INl , KIN2, KOUTl , K0UT2. KCUT3/5.5.6.6, 7/ 

DATA  YES/'OK'/ 

DATA  (flEADY/0.0,0/ 

REUIND  10 

5  READ(10,501.END=10)  J.(MAP(1.J),M,4),(SNAMES(I,JI,I=1,4) 
501  F0RMAT(5Ifc,2X,4A4) 

GO  TO  5 

10  yRITEt)  CIJTl-/>6S) 

666  FORMAT (/•'  LIST  OF  CONTROL  lOPTlONS'/'  NO.  CONTROL  OPTION' 
1/'  0  REPEAT  LIST  OF  CONTRa.  OPTIONS' 

2/'  1  INPUT /COPPUTE  CRISIS  RaOCATION  PLAN' 

2/'  2  RUN  ROUTE  FINDING/LOADING  ALGORITHM'/ 

2'  3  SIMLILATE  THE  EVACUATHDN' / 

3'  4  ASSEMBLE  NETW3RK  FROM  DATA  FILES'/ 

4'  5  READ  STORED  NETW3RK'/'  6  SAVE  GENERATED  NETWORK' 
5/'  7  REASSIGN  INPUT  CHANNEL  FOR  PARAMETER  VALUES' 

6/'  8  REASSIC-N  OUTPUT  CHANNEL  FC«  SUMMARY  TABLES' 

7/'  9  REASSIGN  OUTFUT  CHAM®.  FOR  ROUTE  Ll'^'^ING' 

8/'  99  TERMINATE  PROGRAM') 

20  WRITE(KOIJT1,601) 

601  FORMAT!/'  IIPUT  CONTRa  OPTICW  NUMKR') 

READ(KIN1,»,ERR=20)  I 
IF(I..)€.99)  GO  TO  25 
URITE(K0UTl,608) 

o03  FORMATUX.  'TYPE  "OK"  TO  TERMINATE') 

R£AO(K1N1,6051  AYN 
IF (AYN.se. YES)  GO  TO  20 
WRITE(K0UT1.6O9) 

609  FGRMATdX, 'PROGRAM  ENDS') 

IF (K0UT2.NE. 6)  WRITE (K0UT2. 620) 

IF(K0UT3.NE.6)  WRITE(KiXlT3.620) 

620  FORMAT! ///U) 

1F!K0UT2.NE.6)  BffiFILE  K0Lrr2 
IF(K0UT3.HE.6)  ENDFILE  K0UT3 
RET1.IRN 

25  IF!I.LT.0.0R.1.GT.9)  GO  TO  20 
WRITE!K0UTl,6O3)  I 

603  FORMAT! IX. 'CONTRa  OPTION  '.12.'  TYPE  "OK"  TO  EXECUTE') 
READ(KIH1.605.ERR=20)  AYN 
m  FORMAT ''A2) 

IFIAYN.NE.YES)  GO  TO  20 
IF!I.GE.4.AND.I.NE.6)  GO  TO  40 

TC/T  cn  ("YY  f.ri  Tr»  flrY 
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*»  VW  >  V  »v 

IF(I.LE..3.AijD.NF:EAriYlI).EQ.l)  30  TO  40 
IFd.EQ.c.-r'lD.MREADYlD.EO.l!  GO  TO  40 
IFd.EQ.L)  i^ITE(KOUTl.sll) 

iFii.EQ.:)  ufi:e(;<outi.6.2) 

IFil.E'.i)  wRITEd'Oi.rri.olS) 

'  IT  ‘  ’  L  \  '  'Z  ’■'"T  ■’  '“i  IT  «  <  x 

30  "0  10 

6!!.  crF‘. '*  NO  ) 

FCR^'ATi  BRilR:  .'iO  E.-CGATICN  PLAN' ) 

413  F'j'',i'1i'iT(  ’  EF  P-L'Ri  iC  ROUTES  ) 

40  I1=It1 

GO  TO  (10,30, 10'0,:f0,200,300,300..3l0..320,.320). :i 
30  CALL  SETUP(K.QIJTl,K01JT2,KIN2.NN.NL.NREG,NF13,INFS) 

NRE»DY(2)=l 

IF( INFS.E5.il  NREADY(2)=0 

NREADY(3)=0 

GO  TO  20 

100  CALL  PAS3(KIN2,K01JTl,KGUT2,KOUT3,NN,NL.NPR.NRT,hR:3,PV,K.TCC, 
lAdOOD) 

NREAIiY  131=1 

IFINRT.LE.O)  NF;EA0Y(3)=O 
GO  TO  20 

200  CALL  BIJILD(K1N2,K0UT1,K0UT2,NN,NL,NREG,A.IY) 

240  NREADYill=l 

IF(NN.LE.O.OR.NL.I.£.OI  NREADY(1)=0 

NREADV(2)=0 

NREADY(3)=0 

GO  TO  20 

250  CALL  m:  ( KCUTl , IOjtJT2, NN, NL .HRS,PJ,fffi;T, A, Y, S(  1001) ) 

GO  TO  20 
300  KS=I-4 

CALL  SAVERIKS, MREADY , KINl ,KC4JT1  ,Ki3UT2, NN, NL, NREG) 

IF'I.EO.S)  GC  TO  240 
GO  TO  20 

310  WRlTE(KajTl,60b) 

READO:IIJ1,^,£RR=310)  KI 

IFIKI.EQ. KCUTl. OR. KI.£Q.K0IJT2.CiR.KI.£Q.K,CUT3)  GO  TO  330 

IFIKI.LE.IO.AND.KI.NE.S)  GO  TO  :330 

IF(KI.0T,30)  GO  TO  330 

KIN2=KI 

GO  TO  20 

320  )«ITE(K0lJT1.6O71 

READO:iNl,*.ERR=320)  KO 
IF(K0.E5.KIN1.0R.K0.EQ.KIN2)  GO  TO  330 
I?iKO.L£.10.AND.Ka.NE.4.Hl'lD.KO.i€.7)  GO  TO  530 
IF!I.EQ.5.AND.K0.EQ.KC!tJT2)  GO  TO  330 
IF(t;O.GT,30)  GO  TO  M 
IFCl.EO.S)  KCUT2=KD 
IF(I.EQ.9)  K0UT3=K0 
GO  TO  23 

330  WR1TE(KOUT1,6101 
GO  TO  310 

bob  FOF,«ATin,  INPUT  NEW  CWANilEL  NUMBER  (5=TEF:(1INAL)'' ) 

407  Fr;RMAT(lX,dNPU’  NEW  CHANNEL  NUMBER  (i.=TER!i:NAL) ') 

410  F0R«AT(1<,  ' Illegal  CHANNEL  ■iELEdlOti') 

END 


9®Ri'jUTINE  FASS1KIN1,K0UT1,K0UT2,K01JT3,NN,NL.NPR,NRT,HR4,PV,KTlC, 

ITST) 

C  IINE  PASS 

INCLUDE  SGLIRCE .  INSERT  , 

DIMENSION  NP';i41  ,riPH(4) ,XPV(4) ,XPH14) ,ToT (1) 


DATA  YES/'OK'V 
DO  3  L=1,NL 
UMIT(L,11=.FALSE. 
3  0Mn(L,2)=. FALSE. 


i 

I 
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C  INPUT  MrrW’jRK  PPfiAMETERS 
m  WRITE(KfflJ71,606) 

m  FORNATW IX,  INPUT  VEHICLE5/Lh(€/H0UR  FOR  rilniTED  ACCESS.  2)PRINA 
IRY  4-LANE,  I'^RIYARY  2-LANE,  4!SECC.'ffiAftY  HIGhUAYS' ) 
READ(KIN1,»,EPR=910)  (NFV(K1  .I>:=1,4) 

R20  ySlTE'.KCUTl.SOS) 

cOS  FGRMATdX, 'I>JPUT  AVEFAiCE  IfEED  IN  MILES/HQJR  "CR  liLIMlTEII  ACCESS, 
1  OPRINARY  4-iANE,  SlPRlS'jW  C-LANE,  RidECiNSAfiV  HIGHWAYS') 

READi. F INI  1 ». £nR=y20i  j ,f=1 

:0  2  ■<=! .  4 
dCj'  =‘irViL} 

2  YrH''!  ':=YFr-',.'J 
4  -^i.LL 

f.i=mfr(lwi:oo:) 

r  i = i  ''FR  ( L ;  - 1  1 )  / 1  oco 

V=XPVil) 

T=XPH(1) 

IFIK1.LE.3)  GO  TO  16 
'i'=XPV(2) 

T=XFM(2) 

IFIKI.EQ.4)  GO  TO  16 
V=XF'V\3) 

T=XFH(.3) 

IF(K1.EQ.5)  ijO  TO  16 
U=XPV(4) 

T=XPH(4) 


16  C=2.0 

IF'Kl.EQ.l)  C=3,0 
IF(K1.GT.4)  C=1.0 
C=V*C 
BFRIL)=C 
mNK(Ll=T 
CAF'LNk:(L!=C 
CFRa)=CSTLNK(L)/T 
!F(CFR(L1.LE.0.0I  CFRiDO.Ol 
4  CONTINUE 
360  HRITE(K0UT1,623) 

623  FCRHAK IX,  'INPUT  HIGHWAY  CONFIGURATHON  CODE  DNCSWL,  2)CNE  WAY,  3 
1)  VARIABLE  LAf£S') 

READ(KIN1,»,£RR=360)  kTOC 
IF(KTCC.LT.l.0R,IITCC.C,T.31  GO  TO  360 
IF(.KTCC.LE.21  GO  TO  370 
DO  365  L=1,NL 
361)  BFR(L)=2.0+BFR(L) 

370  WRITE(K0UT1,624) 

624  FORMAT ( 2 X.'IkPUT  PASSENGERS/VEHICLE  X  100') 

READ(KlNl,t,ERR=370)  KPV 
IF(KPV.LT.O.OR.KPV.GT. 10000)  GO  TO  370 
FV=FLOAT(kPV)/TOO.O 

m  WRITE(K0IJT1,629) 

629  FORMAT( IX, -INPUT  NOMINAL  EVACUATION  TIME  IN  HOURS  X  100') 
READ(KIN1,*,ERR=330)  NH 

HRS=FLOAT<W). '100.0 
IF(HRS.LT.0.0.cr,:.HRS.GT.24.)  CO  TO  380 
390  WRlTE(kOUTl,630) 

630  format; IX, 'INPl.IT  DELAY  TIME  TERMINATION  LIMIT  IN  HOURS  X  100') 
PEADlk:iNl,t.ERR=390)  ND 

CTL^FLOAT^NOdOO.O 
I F  ( 2 .  CHfDTL  .  LT .  HP  S 1  GTL=H.RS  /  2 . 0 
WRITE(KCiJTl,650)  (XPV!K),'XPmK),Ksl,41 
IF ( KTCC. EQ. 1 )  write; KOUTl , 651 ) 

IF i 'fTCC. E5.2)  '..RITEikCOTl , 6521 
IF;t'TCC.EQ.3)  WRITEtKCiUr: ■  653) 

_ I TE  ( K  OUT  1 , 6.5  4 )  P  V ,  HR'S  ,251 

650  FORMAT  1/ IX,  HIGHWAY  CLASS  VCiS/LANE/'HR  MILES /'HOUR  ' 6 
IIX,  LIMITED  ACCESS  ',F10.0,F12.0./ 

22X, 'PRIMARY  4-LANE  ',F10.0,F12.0,' 

SIX, 'PRIMARY  2-LRNE  ',F10.0,F12.0X 
41X,  SECONDARY', 7X,F10.0,F12.0) 

■651  format; IX, 'HIOrt^AY  CGNFIGURATICN:  NORMAL  TWO  WAY  TRAFFIC) 

652  format; IX, 'HIGHWAY  CONFIGIRATION:  CTE  WAY  OUTBOUND  TRAFFIC) 

653  format; IX, -HIGHWAY  CCNFIGIIRATION:  VARIABLE  LANES  TRAFFIC) 

654  format; IX.  AVERAGE  PASSENGERS  PER  VEHICLE:  ',F5.2/ 
nx,  NOMINAL  EYAtCUATION  TIME:  ',F5.2.  '  HOURS'/ 

21X.  DELAY  Tir€  TERMINATION  LIMIT:  ',F5.2.'  HOURS  I 
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C 


:F(K0UT2.EQ.2)  C’C  to  4<)0 
CALL  !-!EADER(KCLiT2.2.0) 
wR I TE ( KDJT2 , 650 )  ( XP V ( K 1 . XPH ( K ) , K= 1,4) 
IFiKTCC.EQ,  1 )  hRITEti!.CiJT2,651 ) 

IFvKTiX.  EQ.  2 1  iJIViT£(KCUT2, 652) 

IF ! KTDC.EC).  3)  kiRiTE(KG'JT2,653) 
W3ITEiXiJT:.65Al  P'rf.HRS.DTL 
CALL  nEADER  •  '•0LT2, 2. 1 ) 

A'.'O  .>RITEs:  vL-ltoIC) 


c55  FGFfflK/lX, 
;01  FCRMAT(A2) 
NIT=ij 
NRT=ij 
5  NRS=NRT 
NIT=NIT+1 
0=0 
7  NE=0 


■01 1  E’R-40‘)) 

''TYpr’^0K'''70  GENERATE  ROUTES') 


DO  30  KPR=l.i'FR 
IFiDENFR(KPR).LT.l.O)  GO  TO  30 
JR=IJFR(1,KFR) 

IF(NE.EQ.O)  GO  TO  20 
CO  10  JE=1,NE 
IF(JR.EQ.,£V(JE))  GO  TO  30 
10  CONTIKUE 
20  !€=NE+1 
JEU(NE)=.JR 
WR1TE(11,601)  r£,..lR 
401  F0R(V)T(1X,2I4) 

CALL  PRICE(JR,NN,Yll,IAE),LFR.IFR,tfR,CFR,OMITll,2),)AijDES,IA(1001) 
1,S) 


30  CONTINUE 
CMAX=99??.0 
40  COLD=CnRX 
CI1AX=0.0 
I'JIAX^O 

DO  70  kPR=l,NPR 
IF(CENPR(KFP:).LT.1.0)  GO  TO  70 
JR=I.JPR(1,KPR) 

I3=:i.JPR(2.)PR) 

DO  50  .JE=1.NE 
IFIJR.EQ.JEV(JE))  GO  TO  40 
50  CONTIMJE 
GO  TO  70 

40  IF(Y(IS.JE).GT.COLD)  GO  TO  70 
IF(Y(I3,.JE).L£.Cf'AX)  Gtl  TO  70 
C«AX=Y(IS.X) 
kl1AX=KPR 
jnAX=JE 
70  CONTINUE 

IFIKMAX.EQ.O)  GO  TO  200 

IF(NRT.EQ.500)  GO  TO  200 

NR7=NRTM 

a3.(NRT)=0.0 

CR(NRT)=i:KAX 

l.:R,:,NRT)=Kr’AX 

^R=IJPR(1.KNAX) 

IS=IJPF(2.I^J1AX) 

CALl  ROUTE' is-. iR.IFN0.NRT,I..J;:,CR,Y(1, UMAX), S.LFR.IFR,«FR,CFR,0!1IT( 
11,1). NODES) 


Y(l‘,i,.rAX)=1.0El0 
IFdFf-lCi.EQ.O)  NEw=Nti4+l 

IFiIF14i3.£Q.O.AND.IF3AliElNRT,NRS.lJR).EQ.O)  GO  TO  SO 
75  NRT=NRT-1 
«  TO  40 

30  IF^NRS.EO.O)  CO  TO  40 
nM=TIMERiNRT,FV,HRS,Tl,T2) 

DO  85  k;rt=i,nrs 

IF(klR(2.kRT).EQ.l<J1AX..AND.TY!1.LE.A(KRT))  GO  TO  40 
35  CONTINUE 
GO  TO  75 
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:00  IF(NEW.GT.O)  GO  TO  37 
205  U?3T=1 
DIV=100,.0 
GO  TO  125 
bT  UVGT=0 

'iFiNRT.LE.i'IRS)  GO  TO  155 
Div’=1.0 

"iNIT.GT.i;  DIV=i?.0 


Iii(K?R+i'>AJ)=0 
121  TST;KPR)-1.0E20 
00  loO  KRT^lirJRT 
f,FR=kiR(2.l'^;T) 

A(KRT)=TI«ER(KF:r.PV,HRS.Tl.T2) 

S(KRT)=0.0 

IF(T3T(KPR).LT.A(KRT))  GO  TO  120 

T3T(KPR)=A(KRTI 

ISS(KPR+1000)  =KRT 

130  CONTINUE 

DO  140  KFR=1.NPR 
IFIDENPRlKPRl.LT.l.O)  GO  TO  140 
JRT:=ISS<KFR+1000) 

S(JRT)=DENF'RIKFR) 

140  CONTINUE 
DIV=DIV+1.0 
Ul=2.0yDIV 
W2=1.0-H1 
k>!G=0 

DO  150  KRT=l,NRT 

I F ( ABS ( DR ( KRT 1 -S ( KRT 1 1 . GT , 1 00 . 0 )  KHG=KHG+ 1 
150  DR(KRT)=U1»S(1!F:T)+H2»DR(KRT) 

C  WRITE(6,602)  lDRlKRT),KRT=l,tiRT) 

602  F0RMAT(1X,1CF12.2) 

IF(DIV.liE.500.0)  GO  TO  175 
IF(KHG.GT.0.AND.D1V.NE.100.01  GO  TO  120 
CALL  LOADER(NRT,NL,KTCC) 

160  DO  131  KFR=l,NPR 

131  r5T(i<:Fnj=i.oE2o 
DO  132  KRT=1,NRT 
KFR=IJR(2.KRT) 

IF(A(KRT).LT.TSTtKFR))  T3T(KPR)=A(KRT) 

132  CONTINUE 

DO  138  fflTslrNRT 
kPR=IJR(2,KRT) 

IFiDR<KRT).GT.O.02»DENPR(KFP))  GO  TO  I® 
IFiA(K.RT).LT.1.02+T'3T(KFR))  C>0  TO  133 

133  DRK=DR(KRT| 

DR(KRT)=0.0 

m  134  JRT=1,MRT 
IFdJRC.JRTl.NE.KPR)  00  TO  134 
DR  ( JRT )  =DR  I JRT )  ♦[«  ( JRT )  tDRK  /  DENPR  ( l^PR ) 

134  CONTINUE 
138  CUNTIM.IE 

NR3=0 

DO  144  KRT=MURT 

iF':dr(krt).eq.o.O)  go  to  144 

NRS=NR:54-1 

CR<NRSI=DR(K'RT) 

CR(NRSI=':R(t;RT) 

NF=IJR(UKRT)+1 
DO  142  J=1,NF 
142  I.JF:(J.NRS)=IJR(.J,KRT) 

144  CONTI  a  JE 
?«T=NRS 

CALL  LOADERINRT.M.KTCCI 
IFiLAST.GT.O)  GO  TO  175 

155  CALL  BL0CK(KCIJTl,K0Lrr2-,NN,NL.<TCC.PV,HR3.Tl.NIT.NRT) 
IFIT1.LE.2.0*DTL+.001.0R.NRT.£Q.500)  Gij  TO  205. 

GO  TO  5 

175  CALL  SH0W(Ki:<JT2.K0UT3.NPR,NRT,‘€,PV,HRS,A) 

RETURN 

END 
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SUBROLlTli'lE  ORDER (NN.fl.LFR.  IFR,JTO,?;rR,CFR.  ISS,  !D) 

0  T  ( 1  j , ;  ( 1  i , . JTO  i  I ) .  MFR  i  1 )  -  ■-  FR  '•  i )  •  I  ‘  ‘  D  ^  > 

C  CALL  TI!^E(0) 

oa  10 

10  ID(Ki=lCOCCOrt.jTO!!  l+^IFFin 
CALL  SGRTH :  . j.uFri. 

"ll  ■■■:•  ,  M_  .  .  , 

■  LRITE'R.^CO'  K^LFi^iK),lrR^^)-T^G^K)■AFR.^/,vrrt^^)•.-^i 

'  oOC  FOR?*, AT(1X>5IL.F10. 0.1105 

23  lkIkol=i< 

DO  30  K=l,flL 
IO=iFR(K) 

J.>=jT0lK) 

W:.=«FR\K1 

C.0=CFR1K1 

K0=LFR(K1 

1FR(K1=1FR(K0) 

JTG(K)=0T0(K0) 

MR(K1=MFRvKOI 

cfr(K)=cfr(k:o) 

IFR(K0)=I0 

,rro(Koi=.JO 

MFR(K0)=«0 

CFR1K0)=C0 

K1=I'0S(K) 

ISS(K0)=K1 
30  LFR(K1)=K0 
DO  35  K=1.?1N 
25  l.FR(K)=0 
J=0 

DO  40  K=l,NL 
IF(.JT0(K).EG!.J)  go  to  40 
JS.JT01K) 

LFR(L')=K 
40  CGNTIMJE 
Nl=Nr'l+l 
N2=Nin 
LfRimi=NL+l 
DO  50  K=l.Nl 

lF(LFR(l'-:i).EQ.01  LFR(Ki)=J 
50  J=LFR(K1) 

DO  80  K=1m1L 
ii=inuK) 

I3S(K)=0 

I2=JT0iKI 

L3=LFRIIl) 

LF=LFR(I1+1)-1 

IF(lFRS?E^I2.AND.«FR,tL).EQ.KFRlK))  GO  TO  70 
60  CONTINUE 
CO  TO  30 
70  !3Siiv;=L 
;30  CON'IIUJE 
DO  RO  :<=1,?1L 
L=ISS(K) 

r  KiLFRlKlIcFRU'^iiTOiK) , ISS'Kl .NFR(K) . IFRd) .•JTOILl . ISSil) .i1FR(L) 
■  601  FORNATdX.lOU) 

90  .JTO(K)=ISS(K) 

RETURN 

END 


i. 


J 
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S'JBROIJT  I!€  PR  I CE I  .js ,  Nf) .  V  >  lFR  ,  li-s ,  ;iF'? ,  C?R . I T .  1 1 1  i .  i ; 

C  PRICE  CUT  tiODES  FOR  TRAFFIC  TD  -I.R 
uCCICAL  CMIT 

DII^SICM  •|il),LRR(l;,IFRa).KFR;il,;FR(lMC:‘!IT(l),:CC(2,l),I(i;,N 
1CTES13.1, 

lU  i') 

Ti=vv  go-;'  'j-s 


IS3i2,.JR)=l 
DO  DO  I=2.NN 
11=1-1 

IF'i'o( 1 1) . GT,??999i )  ijO  TO  60 

IT=IS3(2,I1I 

K3=LFR(IT) 

fF=lFRHT+1)-1 

IFIKS.OT.KFI  00  TO  25 

DO  20  K=K3.KF 

:iK=IFRiK) 

i)K=IS3(  1 1  IKl 

600  F0Rf1ATIi:<,5I5,3F9,3,4I6) 

IF(JK.LT.I)  C'O  TO  20 
P=S(Illi-CfR(K) 

IF(P.GE.SH<))  GO  TO  20 
IF(OmT(K))  GO  TO  20 

C  WRITE(11,600) 

C  IK,IT,MFR(K),JK,I,CFR(K),S(ll),S(JK).»]!]£3(l.IK>.NiDD£3(l.IT),N0DES(l.o 
C  K),NijIjES(l.I) 

15  S(.JK)=P 
20  CONTINUE 

25  IF(I.GE.NN)  GO  TO  60 
30  JS=I+1 
J=I 

P=3(I) 

DO  40  .J1=..1S,NN 
IF(S(Jt).GE.P)  00  TO  40 
P=S(.J1) 
v>=Jl 

40  CONTINUE 
S(J)=S(I) 

S(I)=P 

N=1SS(2.J1 

ISS(2.J)=1SS(2,I) 

133(2, I) =N 
133(1, N)=I 
N=ISS(2,J) 

ISS(1,N)=J 
50  CONTINUE 
60  DO  70  1=1, NN 
IY=IS3(2,I) 

70  Y(I't)=3;i) 

C  «RITE(il,t)01)  (Y(I),I=1,NN) 

601  FlRTWT(1X,10E12.4) 

RETlRN 

END 


JUBFOUTINE  ROI.ITE(IS,JR,IFNO,RR,IJR,CR,Y,S,LTO,JTO,!1TO,CTO,OI1IT,NOD 

C  FIND  THE  SEST  DISTINCT  ROUTE  FROM  IS  TO  UR 
LOGICAL  ONIT 


C 

C 

C 

c 

c 

c 


DIMENSION  IJRI21,l),CR(l),Y(ll,S(n 

dKiM  JOL/ii’  lloi/’  ^  >  >CT0(  1)  .OHITI 1  ),M]DES(3, 1 ) 

DO  5  L=1,4C7 
J=.JTO(L) 

N.J=1000OfH'ICDES(1,.J)+10*N0DE3(2,J)+N0DE3(3,U) 
e  L,LTO(L),LTO;L+i),.J,NJ.MTO(L'  ,CT0(L1  ,BTO(Li,Y(J) 

g  LUNTlNUt 

WRITE(6,602)  KR,IS,.JR 
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1FN0=0 

IJR(l,,<R)=l 

J=IS 

CR(KR):=y(IS) 

CTEST=YiIS)»TOL 
SiKR)=0.0 
10  ]L=IJR!1.KR) 

IF lY  i  1 L'T. L’U  lO  7C 
20  I=.J 

:fi:^N£..R.RN::.:L.L£.Ll:^;  00  TC_3p  . . . 

L-  rtRITEi'i  'iOu  j  r.RiUR'.KR) iSlkSI  iMTC(L)  >  i  IuR( liK+l^fR/'  i k=li  »L.' 
IF(IL.GT.LIM)  GO  TO  70 
RETORN 
SO  i;KR=i(KR) 

LS=LTO(I) 

LF=LT0(I+1)-1 
DO  50  L=LS.LF 
IFIOMITID)  GO  TO  50 
J=OTO(L) 

C  WRITER), 602)  I,.J,IS,JR.L.t1T0(L).CKR.CT0(L),Y(J).Y(I3)iY(JR) 

602  F0RMAT(1X,6I9,5F9.3) 

IF(a<R+CTO(L)+Y(J).LT.CTEST)  GO  TO  55 
50  CONTINUE 
70  IFNCi=l 

C  '^ITE(6,602I  1,.J.IS,JR,L.MT0(L).CKR.CT0(L).Y(J),Y(IS),V(.JR) 

RETURN 

55  IJR(l,KRl=IL+l 
I.jR(IL+2.KR)=L 
S(KR)=CKR+CTOiL) 

CR(KR)=S(KR)+Y(J) 

IL=IL+1 

0  URITE(10,600)  KR,CR(RS>>S(KR),MT0(L).(IJR(1.K+1,KR),K=1,IL) 

600  FORHATdX,  15. 2F9.3,20I5/i25X, 21)15)) 

GO  TO  20 
END 


SUBROUTINE  COUNTY(KajT2.NN,XU.X.LLOC.NODES.ANAMES) 

C  LIST  POPULATION  CENTROIDS 

OmeiSION  XLL0C(2. 1 ) .LLOCI 1) .N0DES(3. 1) . ANAMES-:?,  1) , ALPH(6) 
DATA  ALPH/'A'.  'Bd  'C'.  'DS'  d  'X'V 
lF(KOLrr2.N£.6)  CALL  HER[€R(KCUT2,4. 1) 

UR1TEIK0UT2.601) 

DO  10  Jd.NH 

IF(NODES(3,J).GT.O)  GO  TO  10 

K=1000*N0DES(1,J)+2*N0DES(2,J)-1 

P0P=XLL(1C(1,J) 

(«=XLL0C(2.J) 

DEN=POP/AR 
Li=LLOC(J)/lO 
IFILI.EQ.O)  Ll=5 
L2=LUX(.J)-lO+Ll+5 
AL=ALFH(L1) 

CF=ALPH(L2) 

WRITE(KCUT2.602)  K.,(ANA«ESII.J).]=1.7),P[f.AR.DEW.AL.CF 
10  CONTINUE 
RETURN 

602  FORMAT! IX, I5,2X,4A4,2X,3A4,F12.0,2F8.0.;3X.A1, LX, Al) 

601  FORMAT! /IX, ^POPULATION  CENTROID  LIST'/ 

IIX.'FIPS  PLACE  IWHE'.SX. 'COUNTY', 8X, 'POPULATION  AREA', 
2'  DENSITY  FEMA') 

END 


•jJft'.O'.JTIME  iTaCKANET ii'iREC',MAPiKINl.<iXITl,KOIJT2.i'iN.f»OL'jC) XLLuC 
1 LL  vC  •  hiJjDEo  1  ANAilES ) 

;  BUILD  SEL^CEK  ‘JODE  LISTS  FF.uH  r^ODELIST  FILS 

D  L'lE  IS  r  C'N  'ICLOC  ( J  •  1  ^ '  XLLUL 1 2  >  1  j  >  NCISS  i  3  - 1 )  •  rNAMES  ( 7  •  i )  i  LLLU  1 
DI'EMSICN  ^Enu.MREOUi.^^AF'i.l) 

-L" U’<E  FLl'i  ji ■‘vUM 1 1  jiJ'Eii'iRi 
I"  ;  I 
VuLUL  •  1 1 1 1 
XLLCC(2.I)=0 
10  CONTINUE 

o0=0 

K=0 

DO  40  J=1.40 

IF(t1ftP(l,.Jl,EQ.O.ijR.NR£'l(Jl.EQ.O)  GO  TO  40 
NR=I1AP(1,J) 

FIND  (3'NR) 

NF=W(2.J) 

!F(OO.GT.O)  NCiL0Cl2...i0)=K 

N0L0C(1,.J1=K+1 

J0=J 

NNS=NF-NR+1 

;»RITEIK01JT1,S01)  NNS,.J 

501  FURKATi  IX, 'READING  'Mi,  '  NODES  FOR  STAfE  M2) 

DO  30  UR=NR,NF 

K=K.+1 

IF(K.GT.500)  GO  TQ  50 

READ  (3'NR,502,ERR=2C'i  IT,IC.I0.(ANAt1ES(l.K),l=l,7),XIP, 
1XIA,IX 
IC=iIC+l)/2 
XLLCC;i.K)=XIP 
XLL0C(2,K)=XIA 
'J-OC(K)*lX 
N0DES(1,K)=IT 
Ni:iOES(2,K)=IC 
NGDES.(3,Kl=:i;i 
GO  TO  30 

10  WRlTEi3,i02)  J,LR 

502  F0RrtAmX,2I5) 

K=K-l 

30  CONTINUE 
40  CONTINIJE 
NN=K 

NCL0C(2,.J0)=NN 
RETURN 
50  NN=-l 
URITE(6,S03) 

603  FORMAT! IX,' NUMBER  OF  NODES  EXiSEDS  500') 

RETURN 

501  F0RMAT(At.I5,2X,5A4,I3) 

502  F0RMAT(I2,I3.4X.Il,l2X.4A4.26X.3A4,2X.2F3.0,4X.l2) 

503  FORMAT! 1 X, 315, J S. 0.213) 

END 


FUNCTION  LOCK(]T,IC,IO.NOLCC,NOreS) 

C  FIND  NODE  INDEX 

DIMENSION  N0L0C(2.ll,N0DES(3.1) 
KS=NOLOC(1,IT) 

KF=N0LX(2,IT) 

DO  10  K=KS,KF 

IF(N0DES12,K).NE.IC)  GO  TO  10 
IF(N0DES(3.K).NE.I0)  GO  TO  10 
LOCK=K 
RETURN 
10  CONTINUE 
L0CK=O 

WRITE(6,601)  IT,IC,IOtNOLOC(1,IT)>NOLX(2.IT) 
601  FORMAT (IX, 'NODE  UXATION  ERROR'.SIS) 

RETURN 

END 
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SUBROUTirJE  RClAlK;(NET.NREG.MftP.Ki]lJTl,WN.NL.CIJ,  IJ-NiDLCC. NODES) 
C  CONSTRUCT  ROfiD  LINKS  FRCM  R&^DNET  FIUE 

DIMENSION  CIJ(1).IJ(3,1).N0LK(2,1),N0DE3!3.U 
DIMENSION  NET ( 1 ) , NREG( I ) . MAP (4. U 
DEFINE  FILE  9< 10000, lOO.E.NR) 

NL=0 

DO  40  ..1=1, 60 

IF(MAP(3iO).EQ.O.iR.NREC'iJ).EQ.u)  A!  TO  40 
;,R=MAP'3,j; 

NR) 

NNS=2«(NF-NR+1) 
wRITElKOUTliiOl!  1^3, .J 

i.01.  FORMAT (U, 'READING  ,13,'  LINKS  FOR  STATE  ,12) 

BO  30  .JR=NR,NF 

READ  (9'NR.503,ERR=25)  IT.IC,I0.Ll,L2,JT,JC,.J0,L3,Lfi,;<l,K.2,JH 
I3=NREG(IT) 

.jS=NREG(JT) 

IFdS.EQ.O.OR.JS.EQ.O)  GO  TO  30 
Ii:=(IC.H)/2 
JC=(JC+l)/2 
10=10+1 
JO=JO+l 

n=((H-L3)»(Ll-L3)+(L2-L41»(L2H.4))«.5 
I;l=LOCK(  IT,  IC,  I0-1,N0L0C, NODES) 

J1=L0CK<  JT,  JC,  .JO-1,  NOLOC,  NODES) 

IF(H.EG.O.OR.Jl.EQ.O)  GO  TO  30 
KODE=1000*Kl+uH+10000»lK2-l) 

IF(IT.EQ.JT)  GO  TO  13 
DO  10  L=1,NL 

It‘(Il.NE.IJ(l,L))  Ki  TO  10 
IF(J1.NE.IJ(2.L))  GO  TO  10 
IFIKCDE.NE.IJC.D)  GO  TO  10 
GO  TO  30 
10  CONTINUE 
18  NL=NL+1 
IJ(1,NL1=11 
IJ(2,NL)=J1 
IJ(3.NL)=K0DE 
C:j(NL)=D 
NL=NL+1 

lJ(l,NL)=IJ(2iNL-ll 

IJ(2,NL)=IJ(l>NL-l) 

IJ(3,NL)=IJ(3,NL-1) 

CIJ(NL)=D 
GO  TO  30 

25  WRITE{4,-i02)  J..JR 
602  FORMAT (IX, 215) 

30  CONTINUE 
40  CONTINUE 

502  FCF'MAT(2X,I3,I2,2I5,3X,I3,I2,2I5,IX.2I1,15) 

503  FCRMAT(I2, 13, 12,215, IX, 12, 13. 12,215, 1X,2I1. 15) 

RETURN 

END 


SUBROUT INE  SAVE13  ( KFLAG ,  .NREADY ,  K. I Nl .  KCtiT  1 .  K0UT2 .  NN,  NL ,  DREG ) 

C  'SAVE  OR  RETRIEVE  ORDERED  ARRAYS  TO  BYPA'SS  OR'uERING 
INCLUDE  SOURCE. INSERT 
DIMENSION  NREGdI 
DATA  YES/  ' YES-' / 

NREADV=0 

IF(KPLAG.EQ.O)  RETURN 

iF':):flag.£Q.2)  go  to  loo 

5  yRITE(K0UTl,601) 

READ(K1N1,*,ERR=5)  IFILE 
IFdFILE.LT.il. OR.  1FILE.GT.30)  GO  TO  10 
RBiIND  IFILE 

READ  d F I L£ ,  ERR=200 )  NN ,  fL ,  ( ( K ) ,  K=  1 , 6*3 ) .  d  NOLOC  d ,  .J ) .  I = 1 , 2 ) 
1,.J=1,60) 

F:EADdFILE,ERR=2C0'i  (LFR1K),K=1,NN).LFR(NN+1) 
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READ(IFILE,ERR=200)  (IFR(J),J=1.NL) 

R£AD([FILE.ERR=20':))  (JTOtJ).J=l.>l) 

READ  ( IF  I LE ,  ERR=iOO )  I  tiFR ( J ) .  ..1=1 .  .ML ; 

READ !  I F ILE ,  ERR=;00 )  t  DSTL!«  ( .J ! ,  ■j=  I ,  ML  > 

READ(  IF  I  LE .  £RR=2(Xi )  ( ( XLLCC  ( .J.  K  | ,  i<;= 1  ■  .','N  I  -  -J*! ,  2 ) ,  ( LLGC  ( K ) .  K= 1 ,  NN ) 
READ(IFIL£iERR=20O)  ( iN0DE3(Iik)<K=lTl'iN)i  i=l;3) 

REABl  IFIL£iERR=2.30)  I  lANAfESi  IiKi  ,K=1)KN) .  I=li7) 

«RITEit:VjjTl,.i.5l) 

:p  W^N=l-ip . 

.!■ '.NREoit- i.cii.o)  i.'ij  Ti.i  .0 


>RITE:;:2GT1,&S2)  ;<,(SNAf€S(:,i<),  1=1,41 


451  FORMAT (/ IX, 'MUdBER  STATE  NAME') 
c52  ■0R«AT(1X,I4,4X,4V4) 

,.«iTCa.:cuTi.40C)  nn,nl 
REWIND  IFILE 
NREADY=1 

IF(K0UT2.EQ..4)  GO  TO  30 


CALL  HEADER(KaiT2,5.0) 

WRITE(K0UT2,451) 

DO  20  K=l,60 

IF(NREG(K).EQ.O)  GO  TO  20 
WRirE(K0UT2,452)  K,ISNA«ES(1,K),I=1,4) 

20  CONTINUE 

WR1TE(K0UT2,600)  NN,NL 
30  URITElK0UTl,fc04) 

404  FORMAT! /IX, 'LIST  THE  F'OPULATICW  CE.NTFOIDS  ('/ES  CR  Ni3)?') 
READ(KIN1,5011  AYN 

IF(AYN.EQ.YES)  CALL  C0UNTY(K01JT2,NN, XLLOC,LLiX, MODES, ANAMES) 
502  Fi3Rf1AT(A3) 

RETURN 

100  NR1TE(K0UT1,402) 

READ(KIH1,*,ERR=100)  IFILE 
IF(IFILE.LT,11.CR.IFILE.6T.30)  GO  TO  100 
REWIND  IFILE 


URITE(IFILE)  NN,ll.,(NREG(K),K=l,60),((Nll.tX(I,J),I=l,2),.J=l,.401 
WRITEdFILE)  (LFR(K),K=l,NN),Lni;(NN+l) 

WRITE(IFILE)  (IFR(J),J=1,NL) 

WRITEdFILE)  (JT0(J),J=1,NL) 

WRITEdFILE)  (MFR(j),J=l,M) 

WRITEdFILE)  (0STL^^!;iJ),.J=1.^4L) 

WRITEdFILE)  ((XLLCCtJ,K),K.=l,IW),J=l,2),  (LL0C(K),K=1,MN) 
WRITEdFILE)  dN0DES(I,K).K=l,NN),I=l,3) 

.  WRITEdFILE)  (CANAMESd,K),K=l,Nhl),I=l,7) 

600  FORMAT! /IX, 'NUMEER  OF  NETWORK  NODES  ', 15/ IX, 'NUMBER  OF  NETWOR 
IK  LINKS  M5) 

401  FORMAT! IX, 'ASSIGN  INPUT  CWANICL  FOR  SAVED  NETWORK') 

402  F0RMAT!1X, 'ASSIGN  OUTPUT  CHANNEL  FOR  'SAVED  NET’WijRK') 

ENDFILE  IFILE 

NREADY=l 

RETURN 

200  'WRITE!KOUT1,403) 

403  FORMAT!  IX,  FILE  ERRiOR  NETWORK  READ  CANCELLED'  ) 

RETURN 

BID 


SUBRiDUTINE  SH0W!KrjUTl,IFILE.I«>R,ll4T,NE,PV,HRS,XL) 
C  DISPLAY  A  ROUTE 

INCLUDE  SrjURCE.IN'SERT 
DIMENSION  XL!21,1I 
IF!K0UT1.EQ.4)  CX  TO  1 
CALL  HEADERiKOLrri.2.2) 

1  CALL  HEADER! IFILE, 2.3) 

WR]TE!KOUT1,601) 

DO  210  KPR=1,NPR 
DO  200  J^T=1,NRT 
IF(IJR!2,JRT).NE.KPR)  GO  TO  200 
DO  2  J=l,21 

2  XL !J, CRT) =0.0 

IF!DR!JRT).LT.100.0)  GO  TO  200 
WRITEdFILE,609)  JRT 
WRITE!  IFILE,  404) 

WRITE! IFILE, 410) 


74 


100  «:=IJR(l.JRT)+l 

C  WRITE(6.6b6)  JRT, (I..R(l.K.vjRT),K=l,NF), (LlR(2.K,vlRT),K=l,NF) 

666  FORMftTdX.lOlS) 

C  GO  TO  100 

3  L=IoR(NF,wRT) 

I=IFR(L) 

Ni'.'CE=10000*NQDE'i\l'i  I l»10*i2s-N0[£S(2i I)-1)-H'I0DES(3- 1) 

rXN=I 

NCa=NON 

CAPKi.O 

^iE=),.i} 

0 

Z'G’'=0.'} 

:rD=0.0 


c 


VCI.=0.1) 

:il1=HRS/2.0 

T0T=7II1 

MCOE=MFR(L) +10000 

WRITE! IFILE, 600)  NODE.flODE.STOT.DTOT.SPD,  (AWAMES(K,N0N),K=1,5) 
MR=0 
7=11. 0 


J=1 

5  0=0+1 
tij=NF-0+2 

IF(0.EQ.2)  GO  TO  25 
10  L=I0R(NO+1.0RT) 

IF  I  NO. EQ. 2)  GO  TO  12 
LO=IOR(NO,ORT) 

I=IFR(LO) 

GO  TO  20 
12  I=I0PR(2,KPR) 

20  DGT=DSTLNK(L) 

SPD<-.PDLNK(L) 

CAP=BFR(L) 

15  T=T+DST/ftliS(SPD) 

24  3T0T=T 
BTOT=DTOT+DST 
3PD=ABS(SF0) 

VCL=DFR(L)/PV 

TVf1sVCL/(2.0*CAP) 

IFITYM.GT.TIM)  TIf1=TY« 

T0T=S70T+TIM 

NOCE=IOCOO+NC:DE:5(  1 , 1  )+l0*(2+N0D£S(2. 1  )-l  J+NiDOESO,  I ) 

NCN=I 

25  CONTINUE 
XL(NO>ORT)=STOT 
NH=MFR(L)-10(»*(lfR(L)/100;i) 

KH=I1FR(L)/10000+1 

Wfi0E=MFRIL)+10000 

IP'KH.EQ.l)  HRITEiIFILE.604)  Ni3DE.HDDE.ST0T,TI!1,T0T.DT0T,SPD.CAP. 

IVCL.  ( fiNAKES(K.  WON)  .K=l  ,7) .  NH 

IF(KH.E0.2)  WRITEIIFILE.603)  NODE.HODE.STOT.TIM.TOT.DTOT.SPD.CAP, 
1VCL.(ANA(1ES(K.N0N),K=1.7).NH 

IFlKH.Ew.S)  HRITE(IFILE.602)  NOBE.rtjDE.STOT.TISl.TOT.DTOT.SPD^CAP, 

IWLi  ( ANAMES!  K,  NON ) ,  K=  1 , 7 ) .  NH 
30  IF(O.LT.NF)  GO  TO  5 
X1.(2,0RT)=ST0T 

60  WRITE! IFILE.iOS)  N0DE,M0DE.ST0T.TII1.T0T,DT0T,SPD.lANAnESlK,N0N),K=l,5) 
HRITE!IFILE.606) 

606  FORMAT i  IX'/ 


T3=TIHER<0RT,PV,HRS,T1.T2) 

WRITE ( KCUT 1 , 620 )  JRT , I OR ( 2. JRT ) , CR ( ORT ) , ; ANAHEi i K -  fC  1 ) >  K= 1 . 7 ) , 
iriAMES  i  F ,  NGN ) ,  S= 1 , 7 ) ,  DTOT ,  7 1 , 72.  T3 
620  FCRI*AT!lX,i5,2X.:4,2X,F10.0.3X,7A4.4X,7A4,2X,F3.1,3(2X.F3.2)) 
200  OjNTIWJE 
210  CONTINUE 
RETURN 

604  FCR«AT(1X,I6,2X,I5,3F3.2,2F5.1,2F;3.0,3X,7A4,''  INTE^OTATE  ' 
1.13) 


A 

\rt 


603  F0RMAT(lX,I6.2X,I5.3FS.2.2F8.1.2F3.O,3X.7A4.'  U.S.  ROUTE  ' 

1.13) 

602  F0RI^T(1X,I6.2X,I5.3F8.2.2FS.1.:F3.0,3X.7A4.  '  STATE  ROAD  ' 

1.13) 

609  FORMATiaX.  ROUTE  ,14' 

610  FORMAT!  IX,  NODE  CODE  TA'^L  DELAY  '3TAL  GIST.  SPEED  . 
1  CAP.  ;CLS  PLACE  NAME  .rX.  CCi;NTV  .5<.  RijUTE  NAME  ) 

601  FORMAT  I, 'IX.  ROUTE  CHARACTERISTICS  AND  LOADINGS '/CC'X.  FROM',  23X,  TO 
1.42X, 'AVERAGE  HOURS  /IX.  ROUTE  =AIR  EVAiCUEES  .  CITY  ,12X, 
2CfJlJNTY',10X, 'CITY  ,12X.  COUNTY'. 6X,-  DISTANCE  TRAVEL  DELAY 
3  TOTAL') 


END 
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SUBROUTINE  SCRTR^N. IDAT, INDX.L) 
INTEGER  INDX.L 
DINENSICN  IDAT(l),INDX(l.i 
DIMENSION  Ldl 
':i=o 
1=0 


xi=o 

■'=0 

:4=Ci 


..='■'1-1 


I /  iOi  *0  y-O 

■:1=N 

C  'LIIB  THE  ’REE 
ITO  IF'ol.LT.T)  GO  TO  32>) 

K2=K2*2 
32=31/ 2. 

S1=IF!X(B2) 

’4=T4+(B2-S1)»K2 
GO  TO  250 

C  INITIAL  CALCULATIONS 
320  T4=k:2-T4 
! "? 

C  NEXT  TWIG  “ 

350  IF(K1.EG.K2)  GO  TO  940 
Tl=Kly-l 
ki=k;i+i 
01=32 
T3=T2 

C  ADD  1  TO  REFLECTED  SINARV  DOUNTER  AND  CARRY 
400  ri=Tl/2. 

r.F(IFIX(Tl).LT.Tl)  50  TO  470 

Ml=Mld 

T2=T2-B1 

81=31/2. 

GO  TO  400 

C  TWIG  CALCULATIONS 
470  T2=T2+B1 

IF(31.Eg.2.)  GO  TO  550 
C  3-TWIGS  i  4-7WiGS 

IF(T3.LT.T4')  GO  TO  5.S0 
C  4-TWIG 

M1=-1»IH 
GO  TO  b30 

550  IF(T3.LT.T4)  GO  TO  410 
C  ...3-TWIG 
560  I11=MU1 
1=1*1 
Ld)=I 
L(J)=I 
J=.J*1 
C  ...2-TWIG 
610  M1=M1*1 
630  r=i+i 
L1=I 
L(:i)=I 
Lt.J>=I 
LO=J 
J=J*  1 
1=1+1 
l2=I 
L(I)=I 
L(J)=I 
iGO  TO  750 

C  MERGE  LEAVES  AND  BRANCHES 
700  .XI-,1 
L0=.J-1 
L1=L(L0) 

L2=LIJ) 

750  IF(;iDAT(Ll).LE.IDAT(L2))  03  TO  320 
L(L0i=L2  ■ 

770  LCf=L2 
L2=L1L0) 

IF;L2.EQ.L0)  GO  ^0  370 
IF(IDAT(Ll).G7.IDATIL2n  Gi3  TO  770 
L(LO)=Ll 
320  LO=Ll 
L1=L(L0) 

IF(Ll.NE.LO)  GO  TO  750 
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HL0)=L2 
GO  TO  880 
870  LlL01=Ll 
£30  M1=M1-1 

IF'Tll.GT.O)  GO  TO  'iX^ 
iFiMl.Ea.O)  GO  TO  350 
C  GEJIESATE  2^0  HftlF  OF  A  4-TUIG 

GO  T'i  oGil 


;',3a(::=u 
30  'CMTIfJJE 
?E" 


.RN 


iLiEF;CUT;r€  8UILD-;k;:n1-KjXIT1.K0UT2.Mi!,NL.NREG.CIJ. IJ) 

C  ftSSE’IBLE  NETWORK  OOMK-GNENTS  AND  iORT  LINKS 
INCLUl'E  SOLRCE.IMEFiT 
DIi"E.''i£IGN  CIjiD.lJlIM) 

B •  IIENS I LN  NET ( 60 ) > NREb ( I) 

BATA  OK.YES,YS/'OK'.  ’YES',  'V/ 

TO  URITEiK.CUTlicOS) 

603  FCRNATt,"'  INPUT  THE  W.'HBER  OF  STATES  IN  THE  REGION  (LF  TO  10)  ') 
:'EAB(rINli»>Er.R=70'  No 
IFiNb.LT.  1 . GR.NS.S’T.  10)  SO  TO  TO 
WRITEtKOiJTl.SO:'  vS 

no;  fcrh at  i  1  y  .  t  nput  ' .  i  ; ,  s  tate  .•x’laiF-s  ■ ) 

REAB(KINl,*.ERR=TO)  '...EV'K, fiSi 
DO  SO  K>1,NS 

IFi  JSV(K ) .LT.  l.OR.^EY' ■) ! .GT. GO  TO  TO 
SO  GGNTINI.E 
;iO  DO  120  1=1.60 
120  NREG(I)=<i 

■i.=:ITE(KCvTl,6;il) 

DO  130  K=1.N3 
I=.J£V(K) 

.iR  I  T£ : )  GUT  1 1 606 )  I .  ■  SNA.'’'E5  ( .  I ) .  •.'=  1 . 4 ' 

ISO  NRF.Gti;=K 

aRITE(KCUT1.607) 

READ(KIMl,5«;il)  AYN 
IF(AYN.NE.OK)  GO  TO  70 

170  GALL  3TAC.K(NET.NREG.MAP,KIN1,KC1JT1.K0UT2.NN,N0L0C.XLL0C.LLDC, 
INOBiES.AfKAMES) 

GFtNN.LE.Ol  RETURN 

CALL  ROADS (NET, NREG. MAP .kCfiJTl . NN. NL. CIJ, I J. IWLCC, NODES ) 
IF(NL..E.01  RETURN 
WRITE(KQIJT1,604!  NN 
WRITE(KQUT1,605)  ML 

«04  FORMAT!/  NUMBER  NETWORK  NODES  MS) 

605  FORMAT!  NUMBER  GF  NETWORK  LINKS  M5) 

IFiK0Ln’2.E0.6)  ijO  TO  200 

CALL  r£ADERIKCUT2.4.0) 

WRITE!r.CiJT2.601 ) 

DO  ISO  ;-l,NS 
I=vEv!K) 

ISO  WRITE'Ti'ATC.SOo)  I,  !SNAMES(J.  I),  J=l,4) 

.,FrE.K:JT2.SC4)  NN 
»Ri TEL  gut:, 605)  NL 
ZOf  33  210  i.=  l./l 
IrR  1 L  '=I'J!  1  tL) 

■jT0<L;=IUi2,L) 

MFR!L-=IJI3,L) 

210  DSTLNK(L)=CIJiL) 

CTLL  CRC£R!NN, NL.LFR, IFR, JTO.MFR.DSTLNK. IA, lYl 

WRITE(KCUT1,.';.03) 

kEAlI(KIN1.504)  AYN 

IFiAYN.EQ.YES.OR.AYN.EQ.YS)  CALu  CCUNTY 
1 ' KOUTS , NN , XLLOC . LLOC, NODES , ANAME3 ! 

504  FORMAT  I  AS) 

505  FORMAT (/ll.  LIST  THE  POPULATION  CEN'FOIDS  (YES  OR  NO)'  : 

RETuRti 

601  FORMAT! /IX.’MJMBER  STATE  NAME’) 

606  Fi:RMAT(1X,I4,4X.4A4) 

607  FOFMATi/lX, 'T’t’FE  "OK."  TO  ASSEMBLE  NETWijRK  ) 

501  FORMAT 1A2) 

END 


SUBROUTINE  SETUP(K.OUT  1 ,  K0UV2.  KINl ,  NN,  NL.  NREG>  .  INFS  I 

cOGICSL  B 
REfiUB  :,P,R,)(0 

DIMENSION  SV(  100) ,DV(500) .1(600) .IV(500).JV(500) ..JR(500) ,.jE(  100) .P 
ltli')0),R(500). 8(5000) 

INCLUDE  SOURCE. INSERT 
DIMENSION  NREOd) 

DATi)  VES.  CK.  VS/  -  YES  ' ,  'OK'' . '  Y'  / 

:  5  .-=1.NL 
CMIT(L.il=. FALSE. 

FALSE. 

K;=MFR(L)/l()OtXi 

Cl=iMFR(L>-10000^K'l)/10iX) 

T=55.0 

IF(K1.E(1.4)  T=45.0 
IF(K1.EQ.5)  T=40.0 
IF(K1.GT.5)  T=30.0 
SPDLNK(L)=T 
i;FR(L)=D3TLNKIL)/T 
IF(CFR(L).LE.O.O)  CfR(L)=0.01 
5  CONTINUE 
KGE=0 
KGR=0 
K1=0 

200  WRITE(K0UT1,611) 

611.  FORMAT!/, IK, 'INPUT  EVAOJATIGN  PLAN  (YES  OR  NO)?') 

READ():iN1.511,ERR=200)  AYN 
511  FORMAT (A3) 

IF(AYN.NE.YES.AND.AYN.NE.Y3)  GO  TO  220 
)rf?ITE(K0UTl,612) 

612  FORMAT! IX, 'FOR  EACH  PAIR  INPUT:  EVACUATION  03UNTY  NO.,  REOTTION  0 
UXNTY  NO.,  NUMBER  OF  PEOPLE') 

NPR^l 

NE=0 

205  NPR=NPR-1 

IF'NPR.LT.O)  NPRO 
210  CONTINUE 

READ(KINl.t,ERR=210)  Il.Jl.NUM 

IF(I1.EC.999)  GO  TO  125 

IFdl.LE.O)  GO  TO  205 

NPR=NPR+1 

S(NPR)=0.0 

IT=Il/tOOO 

IC=(Il-l000»ITd)/2 

I  JPR(  I ,  NPR)  sLOOK  IT,  IC,  0,  NOLOC.NODES) 

I1=IJPR(1,NPR) 

IT=J1/1000 

Ii2:(Jl-1000»IT+li/2 

IJPR(2.NPR)=LOCK(IT,IC.O,NOLX,NODES) 

J1=IJPR(2,NPR) 

D£NPR(NPR)=NUM 

IFdl.EB.O.OR.Jl.EQ.O)  GO  TO  215 

HRITE(K0UT1,605)  NF'R.DENPR(PPR).(A(«MES(K.11).K=1,7),  1ANAMES(K,J11 
1,K=1,7) 

DO  212  IE=1,NE 
IF(UE(IE).EQ.I1)  GO  TO  214 
212  COHTDJJE 
NE=NE+1 
jE(NE)=I1 

CALL  PRICE( ll.NN,Y( l.NEl.LFR. IFR,11FR.CFR,CI1IT,N0DES, lA.Sl 
IE=NE 

214  S(NPR)=Y!J1,IE) 

00  TO  210 

215  HRITE(K0UT1,651) 

651  FORMAT!  IK, 'COUNTY  WOT  IN  'JETWORK') 

GO  TO  205 

220  HP.ITE(K0UT1,613) 

613  FORMAT! IX, 'SPECIFY  E/AD.lATION  COUNTIES  (YES  OR  NO)?') 
READ(KIN1,511,ERR=220)  AYN 
IFCAYN.NE.YES.AND.AVN.NE.YS)  GO  TO  240 
WRITE(K0!JT1,614) 

614  FORMAT! IX, 'FOR  EACH  EVADJATION  COUNTY  WIT:  COUNTY  NO.,  NUMBER' 

1,  OF  PEOPLE') 
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NE-1 

f''.'jE=  1 

;;5  NE=NE-l 

:f(ne.lt.oi  h£=::i 

230  CONTJMJE 

sehd«I'u.»,err=:.?oi  ii.'iui 
IFill.Ew. GO  TO  2;;-.0 
TFdi.LE.Oi  GO  TO  22*1 
'iE=NE*-l 
■■-1'  ■ 

;>^::-ivi:.)-‘iT+!i '2_ 

-■-i-.LF  uTi  .'.i'.  i;'',l‘i_  jCi^uOcoI 

..E;SE'=J 

2'-0F.E'=-;JLiM 

SUN=NLII1 

IF(J.EQ.O)  GO  TO  235 

wRITEiKClJTl.605)  NE.SVN.  (ftNAllE  ;!K,  J) ,K=1, 7) 

CALL  PRICE( J. NN> Y i  1 . NE)  .LFR.  IFR,«FR, CFR,:jt1IT.  NiODES.  Ift.S) 
ijij  TO  230 

235  K'FITE<ixOUri,451) 

GO  TO  225 

240  WRITEiKQlJTl,6!5) 

615  FGRMRTdX,  ' INPUT  VULNERABILITY  ijROLIP  CODE  (1  TO  4)  ') 
READ(KIN1,»,ERR=240)  KV 
IFIK'V.LT.l.OR.LV.GT.Al  GO  TO  240 

245  nRITE(K0UT1,6411 

641  fOfWAT', IX. 'INPUT  NlNI.tJM  PERCENT  l‘0  K  EVACUATED') 
READ(I':[N1.»,ERR=245)  KPV 
IFiKPV.LT.O.CR.KPV.OT.lOO)  GO  TO  245 
PEV=FLOAT( KPV)/ 100.0 
250  MRITE(KCiLITl.il6l 

616  FORMAT' IX. 'EVACUATE  COUNTOTFORCE  TARGETS  (YES  OR  N0)>') 
F:£AD(KINl,51l.ERR=250)  AYN 

IF ( AYN. £&. YES. OR. AYN. EO. YS)  KF=1 
NE=0 

DO  270  .J=1.NN 

:F(XLL0C(2.J).LE.0.0)  go  to  270 

JW=LL0C':.J)/10 

JF=LL0C(J)-10»JW 

IFiUU.EQ.l.AUD.KV.EQ.l)  GO  TO  260 
IFiJW.GT.l.AND.JW.LE.KV)  W  TO  260 
IFlJF.GT.O.AND.KF.EO.l)  GO  TO  260 
GO  TO  270 
260  NE=NEH 
JEINE)=J 

SV(NE)=XLLOC(l..J)*PEV 

IF(JF.GT.O.ANO.KF.ej.l)  SV(NE)=-SV(NEI 

CALL  PR  ICE  ( J ,  NN ,  Y  ( 1 ,  NE ) .  LFR.  IFR,  MFR,  CFR,  Om  T.  NoOE  j.  I A .  S ) 

270  OUNTINUE 
EMAXd.OElO 
390  NRITEIKOUT 1,626) 

626  FORMAT (IX. 'INPUT  POST  EVACUATION  MAXI.MUI  DENSITY  IN  EVACUATED  COUN 
ITIES  IN  PEOPLE/SQ,  ilILF') 

READ(KIN1,»,ERR=3901  MAXE 
EMAK^MAXE 

2f)0  WRITE(K0UTI,617) 

617  FORMATdX. 'SPECIFY  REYMUrOiN  OT.M''TiS  (YES  OR  NO)'’') 
P,EADa(IN1.5d,ERR=230)  AYN 
IF(AYN.NE.YES.AND,AYH.NE.VS)  GO  TO  'BIO 
ARITE().0IJT1.613) 

61S  FORMATdX.  i^OR  EACH  PECEPTUM  ..dUNTY  IWT:  caWY  NO..  NUMBER'. 

1  OF  PEOPLE' ) 

KGR=1 

NR-1 

290  NR=NR-l 

IF(NR.LT.Oi  NR^i 
300  CONTINUE 

READ':KINl.t,ERR=3<>0>  Jl.NUM 

IF':J1.E0.999)  GO  TO  350 

IF!Jl,LE.i)l  GO  TO  290 

NR=NR+1 

JT=.J1/1000 

JC=(Jl-lOOO*.JT+l)/2 

j-Ll'CK  1 2T ,  oC '  J. NOL'.'U-.  N00“  ■ ' 

DV(NR)=NIJM 

JR(NR)=J 

!;f(j.eq.o)  go  to  m 
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iJR :  TE  ( i(CUT ; .  iOS)  -IR ,  O'J  ■  ^NA'IEi- ;< ,  J !  •  K=1 ,  1 

.JViNR)=l 

:-o  TO  SCO 

:05  WRITE^KQUTT-fcSli 
C’l]  TO  2'?') 

210  >iRITE;((0;JTl,21?i 

?CR?1AT(1<.' INPUT  JIINlUUtl.  -iftXIMUtl  PRE-Ev'aCiATICN  'EC&'T.^'N  A-£o  TE 
1;'::-ITY  IN  Fi-OPL" ' II  P  ■ 

FEAD^KINb't.BiJ-^- 1  TiNQ.MAjiD 

?'n 

220  iRITEv'l4iTl,i20i 

■220  FCR1TAT;i<, 'INPUT  :1ININIJ«.  ilACMU-l  TRAVEL  TINE  IN  HOURS  <  100  I 
F;EAD(t:iNl, »,ERH=32;V!  NlNr-NAXT 
TMIN=FLOAr:(IINn/tOO.O 
TMAX=FLOAT(t1RXT)/100.0 
330  i4RITE(K0iJTl.fc21) 

621  FORMAT  ax. 'INPUT  MAXilfV  i-.U  POPULATION  MULTIPLE  X  100') 
READ(KIN1.»,ERR=330)  MP 
IF(MP.EQ.O)  MP=100000 
PM=FLOATlrtP)/100.0 
fiR=0 

BO  340  J=l,NN 

IF(XLLOC(2,J).EQ.O.O)  'Vi  '  i  340 
DEN=XLL0C(l..J)/XLi.0C(2..J) 

IF'CEN.LE.DMIN.OR.tf-N.GE.OMAX)  00  TO  34(i 
334  DO  336  1=1. NE 

IFlJ.EQ.JEd))  GO  TO  340 
336  CONTINUE 
NR=NR+1 
B'J(NR)=l.OElO 
jR(f«)=J 
340  CONTINUE 
350  CONTINUE 

IFINE.EQ.O.OR.NR.EQ.OI  ij]  TO  4l0 
395  RMAX=1.0E10 

IFlKGR.EQ.l)  GO  TO  405 
400  NRITETKOUTl  oai 

C.27  FORMAT  ax. 'INPUT  =^031  EVACUAri.H  lAXMUM  DENSITY  IN  PECEl'TlON  '.OUN 
ITIES  IN  PE0PCE/3Q.  MIlE'  I 
REAB(KINl.».ERR=4(Xi)  MAXR 
RMAX=NAXR 
405  DO  23  I=1,NE 
J=JE(I) 

IF(SV(I).GE.O.O)  GO  TO  22 

sva)=-sva> 

GO  TO  23 

22  STEST=XLLOr.a-  ji-emax>!'l.  V  ji 
IFlSTEST.GT.SVd))  ■5Vd)=STEST 

23  CONTINUE 
IF(KGR.EQ.l)  GO  TO  29 
DO  27  1=1. NR 
J=.)Fld) 

26  DEN=XLLnC(1.0).'XL.O'<2.  Jl 
IFTRMAX.LE.OEN)  DVd'=j.  i 
PMAX=(RMAX-DEN)MLLXl2.-J) 

JVd)=2 

IFIPMAX.LE.  T  "  >0  TO  27 
IFiPMtDEN.OE.RMA;',  50  '’3  25 
PMAX=(  PM-1 .0 )  ♦0£N»XLLrii' ( 2.  •)) 

IFlPMAX.LT.DVdll  Wd^-/. 

25  IF(FfttX,LT.DU{P  ■  .r  .)=a'a 

27  CONTINUE 

29  fl=NE>l 
N=NR 

4;10  SS=0.0 
DS=0.0 

IFINE.EB.O)  GO  TO  420 
DO  30  1=1. NE 

30  SS=SS+^3Vd  ) 

420  IF(NR.EQ.O)  01  TO  430 
00  40  J=I.NR 
40  DS=DS+DV1J) 

SV(M)=0.0 
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430  WRITE! KOLiri,,);'!.  IE  : 

625  FCRHATI/lX.'NUMBei  ::i-  EVACift'E)  C0Ui'|T:-:?',[4.i(.  E'rtOiiifW  ^OPULAT 
IIQN  •.F12.0/U,'NLIMSER  OF  REC-PflOM  DWKTIEi' .  (-(  t  '-ECr  fluN  A 
2FIEA  CAPACITY', F12.0I 
IFiKGE.BlV  GO  TO  42 
PEV=100.0»PE'J 
URITE(K0IJT1>552'  ov.c?; 

IF(KF.EQ.i)  MR  I.':;  'li 
42  IFiKGR.E'j.l)  GO  TO  44 

WRITE(KGUT1.654l  OH!  I,  jHA)!.  TOIM.  r«<,  ->1 
44  IF(KGE.EO.O)  URirHK:0lin.655>  EHAX 
IF(KGR.EQ.O)  WRlTEiKOUTl.tSi)  RMAX 
IF(I03UT2.EQ.6)  Cfl  TO  44 
CALL  HEADER(l(0UT2,l.O) 

WRITE(K0UT2.625)  NE.SS.HR.DS 
IF  KC€.EQ.l)  ijO  TO  52 
WRITE<K0UT2,652)  KV.PEV 
IF(KF.EQ.l)  HRITEIK0UT2,45;3) 

52  IFlKGR.EQ.n  GO  TO  54 

WRITE(K0LIT2.654)  I]|1IN,DMAX.Tf1IN.T11AX,Pi1 
54  IF'KGE.EQ.O)  HRlTE(KafT2,655)  E(WX 
IF(K:GR.EQ.0)  URIT£<K0in'2.656l  RIWX 
4fa  URITE<K0IJT1.457) 

READ(KIN1,512,ERR=44)  AYN 
IF(AYN.!E.iDK)  RETURN 
512  FORMAT! A 2) 

b52  FORMAT! IX, 'VILTERABILITY  GROUP  CODE:  M2/ 

IIX,  MINIMUM  EVACOUATED  POPULATION:  ',F4.0,'  PERL'ENT') 

653  FORMAT! IX, 'COUNTERFORCE  TARGETS  Af^  EVACUATED') 

654  FORMAT!'  MINIMUM  RECEPTION  AREA  DENSITY:  ',F6.0.'  PEOFIE/SQ  'ILE'/ 
IIX, 'MAXIMUM  RECEPTION  AREA  DENSITY:  ',F6.0,'  PEijPLE/95  MILE'/ 

21X, 'MINIMUM  TRAVEL  TIME  TO  RECEPTION  AREAS:  '.F6.:'-'  HOURS'/ 

31X, 'MAXIMUM  TRAVEL  TIi'E  TO  RECEPTION  AREAS:  ',F6.2,'  HajRS'/ 

41X, 'MAXIMUM  RESIDEr>T  PGPULATION  MULTIPLE:  ',F6.2) 

655  FORMAT! IX, 'POST  EVACUATION  MAXIMUM  DENSITY  IN  EVACUATED  AREAS:  ', 
1F6.0) 

656  FORMAT! IX, 'POST  EVACUATION  MAXIMUM  DENSITY  IN  RECEPTION  AREAS:  ', 
1F6.0) 

()57  FORMAT! /IX, 'TYPE  "OK"  TO  COMPUTE  EVACUATION  PLAN') 

K1=0 

K2=0 

<6=0 

IF!NE.EQ.0.CR.II?.EC!.0,0R.S3.GT.0S)  GO  TO  140 
IFIKCIUT2.NE.6)  CALI  HE(3DER!K0UT2,1.U 
WR:ITE1K0UT2,631) 

Di'J  450  IE=l,NE 
I=UE!IE) 

REC=0.0 
DO  440  J=l,NR 
.IY=JR!J) 

II-1Y(JY,IE).C€.T111N.AND.Y!.JY.IE).LE.TMAX)  REC=RB>DV1.J) 
DEN=XLL0C<1,I)/XLLXI2,I) 

440  CONTINUE 

WR’.TEIK0UT2, 62S)  ! ANAMES  K,  1 ) ,K=5,7)  ,XLLOC!  1,1), XLL0C!2, 1 )  ,[EN, SVl 
1IE),REC. 

450  CONTINUE 

631  FORMAT  I  /  IX ,  'EVAiCuATION  COUNTY ' ,  4X .  'RESIDENTS' ,  bX , '  AREA ' , 

13X,  DENSITY', 3X,  AT  RISK  Fijp. ' -2X-,  REC.  Cftf'ACITY') 

623  FORMAT' lX,3A4,SX.2F10.0,Fl0.2.2Fi5.0) 

629  FORMAT(/IX,  'RECEPTION  CaiNTY'.SX.'RESIIENTS'  ,6X,'AREA', 

13X,' DENSITY ',6X,  CRITERION', 2X,'REC.  r^AiCITY') 

WRITE!K0UT2,629) 

DO  460  J=1,NR 
wY=JR!J) 

DEN=XLL0C(I,.JY)/XLUX(2,.jY) 

IF!..'ViJ).EQ.l)  WR1TE1K0UT2,63:2)  !ANAI€3IK,JY),K=5,7:i,XLL0C!l,JY).X 
lLLCCi2,JY),DEN,IiVI.J) 

IFiJV!.J) .E0,2)  WRITE!KCUT2.693)  !ANAHE3IK. JY) ,K=5. 7) , XLLOC!  1, JY) , X 
1LLGC!2,JY),DEN,DV!.J) 

1F!.JV!.J).EQ.3)  WR!TEIK0Lrr2,634l  (ANAMESIK,JY),K=5,7),XLLCC(1,JY),X 
1LL0C!2,JY),DEN,9V1J) 

632  F0RMAT!lX,3A4,SX,2F10.i)rF!0.2,10X,  'GIVEN',F15.0) 

633  F0RMAT!lX,3A4,3X,2F10.0,F10.2..:X,’I€NSirY',Fl5.>i) 

634  FORMA"  1<>3A4,8X,2F10.0-F10.2.7X, 'MlJLTlFLE  ,F15.0) 

460  iXNTINUE 
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IFiSS.LT.Ki 
DO  50  1=1,  :€ 

:0  50  J=1,NR 

.:y=.R(.!) 

A(cA)=y(..v,I) 

:F■A^Al.LT.•:^-!IN.OR,A(JA).C.T.T^■A:;;  Ai.JA)=lCOO.O*Ai JA) 

5(  c;nt;'..e 

70  A(JA)=0.0 
i::=i 

call  TRHFL)( I :i .K1  ,K2>K.£„XO,ii,M.B, a, LV.DV, 2iP,R,3,0(l01 : ,DFR, 

IIV.jV) 

IF'Xl.NE.l)  GO  ■"0  lAO 

N0=(1+l*-l 

NPR=0 

DC  120  J=l.NO 
IF(Z(0).LE.1.0)  GO  TO  120 
IE=IV(.J) 

IFdE.EQ.M)  GO  '0  120 
.  NPR=NPR+1 

IJPR(1.NPR)=JE(IE) 

IR=JV(J) 

:JPRi2.1PR)=,-!R(IRl 

1)EHPR(NPR)=Z(.J) 

C  «UTE(fc.667)  J,lE,IR.NPR,IJPR(l,NPR),IJPR(2,rff>R),[0iPR(NPR) 
JA=M*(IR-1)+IE 
S(NPR)=A(JA) 

467  FORMAT! n, 614, F12.1) 

120  CONTINUE 

135  IF  NPR.LE.Ol  GO  TO  UO 
1NF£=0 

WRITE(Km!Tl,6351  i'lPR 

635  FORMATI/lX.  'EVACUATIijN  PROBLEM  IS  FEASIBLE'/IX, 
r NUMBER  OF  E'YADJATICN/F'ECEPTIiON  COUNTY  PAIRSM4) 

IF(K0UT2.ME.6)  CALL  HEA[eR(K0LIT2,1.2) 

WRITE  (iMjUTL,  404) 

604  F0RMAT(/1a,  'EvAi:LATI0N/R£CEPTACN  CajNTY  FAIRS'/IX.  FAIR  PECFLE 
1  ORIGINMCX,  'OX.VTYMOX,  '[■cSTimiCN',5i(,  ODUNTYMIX,  ' HOURS' 
2,7X, 'MILES', 5X, 'DENSITY') 

DO  130  J=l,i'FR 
IF(J.EQ.KPR)  GO  TO  127 
JS=J+1 

DO  126  K=JS,NPR 

IF!S  JI.GE.SiKJ)  GO  TO  126 

Il=IdFR(l,J) 

Jl=IviPR(2,J) 

S1=S(J) 

P1=DENPR(J) 

IJPR(l,J)=IUPR(l,K) 

IJPR(2,J)=IJPR(2,I<1 

S(J)=S(K) 

DENPR!J)=DENPR(K) 

IUPR(1,K)=I1 

iOPRi I.  =v,l 

S(K)=S1 

DFtfR(;0=Pl 

126  CONTINUE 

127  Il=LlFR(l,.i) 
j1=I.JPR(2,.J1 
DO  12S  IE=1,NE 
IF(JE(IE).EQ.[1)  m  TO  12'? 

128  CONTINUE 

12?  CALL  ROLiTE(Ji,Il,IFNO,l,IJR,CR,Y(l,IE),A.UR,IFR,MFR,CFR, 
10MIT(l,U,NCiDE51 
D3T=<).0 
)#=IoR(l,l)+l 
CO  135  J0=3,NF 
L=IJR(.jO,l) 

135  DST=DST<-D£T,.NK(L) 

DEN=XLLCC(1,J1) 

DO  150  K=1,NFR 

IF ( I DPR \ 2, K ) . EQ . J 1 )  DEN=DEN+DENPR ( K ) 

150  CONTINUE 
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DEN=BtJJ/XLLCC!2..il) 

,  *^7 f '■'’ *'•”’■  ( K.  1 1 ) ,  K=! ,  - ! ,  ( ANAMES i K,  J1 

i  t  /  /  1  0  U  M  Do  I  )  C  t  J 

;y!il  F6kl-AfaX,I3,2i3,F3.1) 

_  RETURN 

lA'J  wRITE'T-T.'T  1,602)  KliiQ-K*, 

o02  FORrATdX,  EVACJATICN  FRCSLEM  IS  ^CT  FEASTELZdTIc; 

end""" 


SUBROUTINE  LOADL(L,NRT.HRS.P'X,XL,TIMES,NX.NTOT) 

C  LOAD  LINK  L 

DIMENSION  XL'21,11,TIMES(15).XN(15),NX11) 

INCLUDE  SOURCE. INSERT 
DO  10  KT=l,14 
10  XN(KT)=0.0 
K  50  JRT=1,NRT 
IF (DR(URT).L£. 100.0)  GO  TO  50 
1JF=IJR!1.URT)+1 
DO  50  0=3. NF 
Jl=NF-J+3 

IF(IJRIJI.ORT).NE.L)  GO  TO  50 
00=01-1 

TYM=TIMER(0RT,PV,HR3.Tl,T2) 

TYM=2.0»T2 

DEN=DP(JRT)/(TYM»PU) 

DO  45  KT=1.14 
T=TIMEi(KT) 

XAD=0.0 

IF\XL(J0,0RT).EQ.XL(01.JRT)1  GO  TO  40 
25  XAO=OEN 

IF(T.C,T.XL(00,JRT).At'10.T.LT.XL;01,ORT)+TYM)  GO  TO  40 
XAD=0.0 

IF(T.LT.rLIOI.JRT>.OR,T,GT.XLI^lO,.,RT)+TYM)  GO  TO  40 
IFIT.iiE.XLIJl.ORTl+Trtl)  CO  TO  30 
XAD=UT-XL(01.JRTI)/(XL(0i;i,..I?T)-XL(Jl,.JRT))>»iei 
GO  TO  40 

30  XAD=((XL(OO..JRT)+rffl-Tl/(XL(00,ORT)-XL(01i,.TiTl  ))*DEN 
40  CONTINUE 

XN(KT)=>fNlKTl+XAD 

C  WRITE(6,SOO)  ORT.0.01.T.XAD.XL(JO,..fiT),XL(01.0RT) 

45  continue 
50  continue 
NTijT=0 

DO  60  KT=1,14 
X=100.0»XN(KT) 

X=X,'BFR(L)«'.4'?9'R 

NXIKT)=X 

60  NT0T=NT0T+NX(KT1 
liO  P0RMATiU.TI5.4F 1 2. 4) 

RETURN 

END 


INCLUDE  SOURCE. Insert 

’  JfO'  15.  l ) .  XS(  15. 1 ) ,  NX !  15 ) 

Ro^?=l.  (f  -  24. .  32. ,  40. ,  4:3.  ,60. ,  72. .  34 

DO  70  .J=  1,500 
90  XD(I.J)=0.0 
X3{I,1)=0.0 
XSiI,2)=0.0 
100  X3(l,3)=<).0 
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NE=0 

NR=100 

DO  200  i-:RT=l,l'IRT 
KF'R=IuR'-!2>i-?'Tl 
JR=I,.'PRiLKf'S) 
I5=:jFR(2,KPR) 
IFvIlE.EQ.Q)  Gij  TO  120 
30  _ 


Tr-i'MEjsJR 

.. 

130  IF  (SR.  1:5. 100)  DO  TO  iTO 
DO  140  .1=101.  NR 
IF(I5.EQ.I3S(J))  00  TO  160 
140  CONTINUE 
150  HR=NR+l 
IS3(NR)=1S 
J=NR 

160  DO  190  KT=1,14 
T=TIt'',ES(KTl 

T3=TIMER(KRT.PV,HRS,T1,T2I 

FD=1.0 

IF(T.GE,.2.0*T2)  DO  TO  170 
FD=T/(2.0*T2) 

170  FA=1.0 

IF(T.GE,.Tl+2.0tT2)  GO  TO  130 
Ffi=0.0 

IF(T.LE.Tl)  GO  TO  180 
FA=(T-T:i)/i2.0i*T2) 

ISO  XDIKT. I ;)=XD{KT, I)<-Fr:»DRlKR7) 

XD ( KT , J ) =X0 (KT . 0) +FA»DR (KR7) 

X3(K7,l)=XS(KT,l)+(1.0-FD)*a9(Fm) 

X3  ( KT.  2  )=XS(  ter,  2 ) + ( FD-FA )  tCRCKRT) 

X3(KT.3)=XS(KT.3)+FA*DR(KRT) 

190  CONTINUE 

XD(15,I)=X0(15,I)+DRIKRT) 

XD(lS.J;=XDil5..J)»-DR(KRT) 

X':.a5.1)=XS(15,l)*DR(KRT) 

200  CONTIWE 

WRITE I KOUT 1,611)  XS(15,1) 

IF(K0UT2.£Q.6)  GO  TO  210 
CALL  HEADERIK0UT2.3.0) 

URITE(K0UT2.6li;i  XS(15,1) 

210  OO  230  K=1.3 
DO  220  r;T=l,14 
<=100.0«XS(KT.K)/XS(15,i) 

220  NX(KT)=X 

IF(K.EQ.l)  HRITE(::0UT1.612)  (NX(k:T).KT=1,  14) 

IFiK.EO.2)  WRIT, EO'OUTl. 613)  (NX'KTl .KT=1. 14) 

IFIK.EQ.O)  URITE1K0UT1,614)  (NX(KT),KT=1,14) 

IF(K0UT2.EQ.6)  GO  TO  230 

IFIK.EQ.I)  URITE(K0UT2.612)  tNX(KT),KT=l. 14) 

IF(K.EQ.2)  WRITE(nOUT2,613)  (NX(i'J),KT=l.  14) 

IF(K.E0.31  WRITE(KCUT2.614)  (NX(KT),KT=1, 14) 

230  CCNTlNuL 

till  FORMATI/IX, 'EVACUEES:  ^F10.0,5X,■■ElJAPSED  HOURS:  2  4-  , 

r  6  8  12  16  20  24  32  40  48  60  72  34  ) 

612  FORMATflX, 'PERCEv'TAGE  OF  EVACUEES  AT  RI3KM0X,14I5) 

613  FORMAT  I IX. 'PERCENTAGE  I'jF  EVAiacES  EN  RiME'.9X.  1415) 

614  PCRMAT( IX, 'PERCENTAGE  OF  EVACUEES  AT  DESTINATION'. 3X, 1415) 
:F';kout;!.,ne.-‘''  call  header(kout2,3.i) 

WRITE(KCUT2  ol5) 

DO  250  1=1,  r€ 

DO  240  );T=1,14 
X=100.0»XD(KT.I),'XD(15,I) 

240  NX(KT)=X 
.R=ISS(I1 

250  WRITE(K0UT2.616)  iAi'JAMES(K,.jR),K=5.7).XD(15.I) .  (NX(KT) .KT=1.  14) 
I^.r-).T2..NE.6)  CALL  iHEAItR(KCirr2i3.,t) 

,RITEiK0LrT2.617) 

DO  230  J=101,SR 
CO  270  i;t=;,i4 
X=IC0.0*XD(KT,J))XD(;5.J) 

270  NX(KTt=X 
IS=IS3IJ) 
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230  '.iRITE(K()IJT:,t.li)  (ANfiMES(K,IS).K=5,7)..)(D(!5,J),  (NX(:<T)-KT=1,;4: 
615  FORMAT (.'IX,  EVACUATION  DEPARTURE  SUMMARY'/ 

1!.X,  EVACUATION  COUNTY  EVACUEESM9X.  PERCeiTAGE  DEPARTED  AT 
iLAF'UED  HQUR’/SOX,  ’  2  4  6', 


16  20  24 


40 


60 


84') 


616  FORMATdX, 3A4,5X,F10.0.2X,  14151 

617  F:RMAT(,'1.(.' RECEPTION  ARRIVAL  SUMMARY'/ 

11“,  kEUeiFTIuN  ■.^^.‘NT''  EVACUEES' ,,^OXi ■  rE,''.LEiVTA'jE  nSRIVtD 


IF',K!jUT2,fiE.o)  '-ALL  HEADER (KuLlT2i 
NRITE(KOUT2,600) 

LRITElkGlJT2,601) 


40^ 

2,S1 


84' ) 


DO  30  .1=1,  NN 
L3=LFR(I) 

LF=LFR(I+1)-1 
CO  20  L=LS,LF 

IFiDFR(Ll.EQ.1.0E5)  GO  TO  20 
IF(DFR(L)/(P'/*eFR(LI).LT.HR31  C-O  TO  20 
J=IFRi:L) 

NH=ffR<L)-1000*(l1=R(L)/1000) 

KH=MFR(Ll/1000i>H 

CALL  LOADL(L,NRT,HF3,PV.XL.TI(1ES,NX,NTOT) 

IF(NTOT.LT.l)  CHj  TO  20 

IF(KH.EQ.l)  WRITE(KOUT2,604)  (ANAfC3(K,J).K=l,4),(ANAMES(K,I),K=l, 
14),NH,D3TLNK(L),SPDLNK(Ll,BFR(L),lNX(KT),l':T=l,14) 

IF'IKH.EQ.D)  WRITEIK0UT2.6JJ3)  (i»IAME3(K,J),K=l,4!,(ANAMES(K,I),K=l, 
14),NH,DETLM<(L),SPDLNK(L),BFR(L),':NX(KT),KT=1,1^) 

IF(KH.EQ.3)  WRITEiK0UT2,602)  (ANAMES(K,J),K=1,41, iANAMES(K,I),K=l, 
14),NH,DSTLNK(L),SPDLNK(L),EFR(L),(NX(KT),KT=1,14) 

604  F0RMAT(1X,4A4,:JX,4A4,3X,'INTER3TATE',I4,2^7.1,F8.D,14I4) 

603  FORMAT!  1.X,4A4,3X,4A4,3X,'U.S.  ROUTE', I4,2F7.1,F3. 0,1414) 

602  F0RMAT(1X,4A4,3X,4A4,3X, 'STATE  ROAD'. I4.2F7.1.F3.0, 1414) 

20  aMTINUE 
30  CONTINUE 

600  FORMAT! /IX, '■^ETWORK  LINK  CHARACTERISTICS  AND  LOADINGS') 

601  FORMATdX, 'FROM', 15X, 'TO', 17X, 'ROUTE  NPf1E'.4X,'  DIST.  SPEED', 

1'  CAP.',6X.'PERCS4TAGE  OF  CAPAiIITY  IN  USE  AT  ELAPSED  HOUR', 
2/75X,'  2  4  6  3  12  16  20  24  32  40  4?  60  72  34') 

RETURN 


END 


SUBROUTINE  LOADER (NRT,NL,KTCC) 
C  LOAD  NETWORK  LINKS 

INCLUDE  SOUWIE. INSERT 
DO  5  L=1,NL 
DFR(L)=>).0 
5  CONTINUE 
DO  130  JRT=1,NRT 
IF(DR(JRT).LE.1.0)  GiO  TO  130 
NF=IJR(l,cRT)+l 
DO  120  J=5,NF 
L=IJR(J,.RT) 

LO=JTOILi 

110  DFR(L)=DFR(L)+DR(.JRT) 
IF(KTCC.1£.3)  iiO  TO  120 
115  DFR(LO)=DFR1LO)+DR(URT) 

120  CONTINUE 
130  CONTlNLiE 
RETURN 
END 
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^ .  SRvi' ' :  ?L :  :  f.  COT  1 .  ciiT2 ,  an  ,  ,  KTc: ,  F V ,  HR': .  T 1 ,  N I T ,  ss :  1 
lELETS  CAPACITATED  i.pr5 
INCLUDE  SC'JRCE.  I^wERT 
Tl=0.i') 

1  :':in 


"2. LE. T! .1  vi.i  lO  O.' 


-'>L 

:0  CONTINUE 

;1F(LO.E0.O)  RETURN 
IFINIT.EQ.l)  URITElKOUTl.bOl) 

IF(NIT.Eu.l.AND.K0UT2.NE.&)  NR!TE't(0UT2,i01) 

T3=0.?0*T1 
DO  50  1=1, NN 
LS=LFRiI) 

LF=LFR(M)-1 
00  50  L=LS,LF 
IF(0«ITIL,1))  GO  TO  50 
T2=dfri:l)/(pv*bfr(L)! 

IFa2.LT.T3)  00  TO  50 
IF(CFR(L'.EQ.l.0E5)  00  TO  20 
J=IFR(L) 

NH=tlFR(L)-1000*(MFR(L)/iO0O) 

f;H=«FR(Ll/ICCOO+l 

T0=r2/2.0 

VCLS=DFRIL)/PV 

IF'LO.GT.O.AND.KH.EO.l)  UR[TEiK0UTl,602)  NIT,NP.T,BFR(LI ,VCLS, TO, 
1'.:ANAI1ES(K,  J)  .K=l,7l ,  (ANAMEOiK,  I)  ,K=1,7)  .NH 
IF''  0.GT.0.AND.KH.EQ.2)  yRITE(K0lJTl,603)  NIT-.NRT.BFRiD.VCLS.TO, 
l(ANAnES(K,J),K=l,7),(A,NAMES(K,I),K=l,7),NH 
IF(LO. GT. 0. AND.KH. EQ. 3)  «RITE(K0UT1 . E04)  NIT, NRT, BFRd. )  > VCLS, TO, 
1(ANA(1ES(K,JI,K.=1,7),(ANA«ES(K,I).K=1,7I,NH 
IFdO.En.O.ANO.KH.EO.I)  WRITE(K0IJT1,6'05)  BFR(L),VCLS.TO, 

1 ,  A  NAMES  ( ,  J  > ,  r;  =  1 , 7 ) ,  i  AN  AMES  ( K ,  I ; ,  K= 1 , 7 ) ,  NH 
I F 1 10 .  E2 . 0 .  AND .  KH ,  EQ .  2 '!  mR  I TE  ( 1  OUT  I ,  fc06 )  BFR  '•  L  i ,  VCLb ,  TO , 
2(ANAMEStK,J),K=l,7),iANAMES(K,I),K=l,7),NH 
IF(L0.E{1.0.AND.KH.EQ.3)  WRITE(K0UTl,fe07)  BFR(L),VCLS,TO, 
IIANW1ES(K,J),K=1,7),(ANAHES(K.I).K=1.7),NH 
IF1KOUT2.E0.6)  GO  TO  30 

IF(L0.GT.0.AND,K'K.Ei3.i;t  MRrTE(K0lfr2,602)  NIT,NRT,BFR(L),VCLS,TO, 
l(ANAMES(K,J),K=l,7),(ANArtES(K,I),K=l,7),NH 
IF(L0.G7.f).AND.KH.E0.2)  WRITE(K0UT2,603)  NIT,NRT,BFR(L),VCLS.TO, 
1(ANAMES(K,J),K=1,7),(ANAMES(K,I),K=1,7),NH 
IF(L0.GT.0.AND.KH.B3.3)  WRITE(K0UT2,604)  NIT, NRT, BFRd). VCIS, TO, 

I  ( ANA,rES(K.  JI,K=1,7),  ( A-WESiK, I ) ,K=1 ,7)  ,NH 
IFlLO.Ea.D.AlC.KH.EQ.l)  WRlTE(K0UT2,fc05)  BFR(L),XLS,TO, 

l(ANAMES(K.J).K=l.7),(ANAME3(K,I),K.=l,7l,rJH 
IF(LO.EQ,O.AND.KH.EQ.2)  WR[TE(K0UT2,606)  BFRID.VCIS.TO, 
2(ANAlMES(fc,J).K;=l,7),(ANAf1ES(KM),K=l,7),NH 
IF(L0.EQ,0.AND.KH,EQ,3)  '4RITE(K.0UT2,«)7)  BFR(L),VCL3,T0, 

I I  iNAilES  ( K ,  J ) ,  K=1 ,  .  i  AiMES  ( K .  I ) ,  K=  1 . 7 ; .  NH 

30  ijMlT(L.l)=.TRl)E. 

L0=0 

IF ( KTCC. EQ. 2. AND, T1 . GT. HRS )  BFR ( L ) =2. 0»BFR( L) 

40  K=,JTO(L) 

CMIT(K.2)=.TRUE. 

IF(KTCC.EQ.2.hN11,T1.GT.HRS)  BFR(K)=1,0E-5 
IF(KTCC,Ea.2)  0MIT(K,1>=.TRUE, 

50  CONTINUE 

401  F0RMAT(/1X, 'PASS  ROUTES  CAPACITY  VEHICLES  DELAY  FROM', 

I12X, 'COUNTY', lOX, 'TO'. 14X. 'OOUNTY'.lOX, 'ROUTE  NAME') 

402  FORMATI IX, 13. 17, IX. 2F10.0,F7, .,3X,7A4. 4X,7A4, 4X. ' INTERSTATE  , 14) 
‘.03  F0RMAr<lX,I3,I7,lX.2Fl0.0,F7.2.3X.7A4,4X,7A4,4X,'U.S.  ROUTE',  14) 
i04  F0RMAT(IX,13.I7,1X,2F10.0,F7.2,3X,7A4,4X.7A4,4X.'3TATE  ROAD', 14) 
.;05  FORMAT(L2X,2F10.0,F7.2..5X.7A4,4X,7A4,4X,  INTERSTATE',  14) 

4U4  FriiMAT(12X,2F10.0,F7.2,3X.7A4,4X,7A4.4X,'U.S.  ROUTE', 14) 

407  FC««ATll2X,2F10,0.F7.2,3X.7A4,4X,7A4,4X.'STArE  ROAD', 14) 
IFITl.LT.HRS)  T1=HRS 
RETURN 
END 
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ACTION  IR:AMEL:RT,;'IRT, L:R) 

LlflENSICfi  kiR(21.1) 

IFSA«E=0 

IF(MRT.LE.O)  RETUR!>1 
N..=IJR(l,wRT)+l 
DO  20  t..RT=l,f®T 
NK=IJR(1.::;RT)+1 
IF:>iJ.SE.fK)  GO  TO  20 
DO  Ij  I-D.'ij 

^F. I..R( D1..FT! .,‘DE. I JR(I i .'RTl )  00  TO  DC 

IFSAI1E=1 
RETUrri 
dO  DlNTIM'E 
RETURN 
END 


nmiGN  TIMER(KRT,  PV. HRS, T1 .  T2) 
C  COI^>UTE  LONGEST  TIME  ON  ROUTE 
INCLUDE  SOIKCE. INSERT 
T1=CR(KRT) 

30  T2=HRS;2.0 
NF=1JR(1,KRT)+1 
DO  40  J=.],NF 
1=IJR(J,KRT1 

T3=DFR(L)/(2.0»PV»BFR(L)) 
IF(T2.GE.T3)  GO  TO  40 
T2=T3 

40  CONTINUE 
TIMEF;=Ti+T2 
F£TURN 
END 


3J3R0UTINE  HEADER IIFILE.XNUM) 
C  ATTACH  HEADIER 

WRITE! IFILE, 401)  XNUM 
601  Fi:iRI1AT(//XlHl. 'TABLE  '.F3.1) 
FETLIRN 
END 


SUBROUTin  TRNPLX(Il,Kl,K2.K6,X0.n,N,B,C.30.D0,X.P,R,S,D,Y,I,J) 
MAIN  CALLING  SEQUENCE 
LOGICAL  B 

REAL*8  X,F,R,X0.R«,A1,DI 

CIMESSICN  3\l),Clll,3ill,D(l),X(l),P(l),R(l),S0':i;.D0(I),Y'vI),; 
1,J(1) 

L=M+N-1 

T2=1.0E-6 

T3=-1.0E-6 

K2=0 

i;7=M»N 

R0=-1.0E19 

TC=0.0 

DO  5  10=1, H 

T0=T0vS0(I0) 

DO  5  ,.0=1,  N 
KO=M*(.JO-1)+IO 
IF(RO.lT.C(KO))  R0=C(K0) 

5  CONTINUE 
PO=1.01»R0 
T2=T2*T0 
T3=T3»R0 
KF‘?=0 


GO  10(10.201.11 

10  CALL  1EWIM.M.L.RO.T2.T3.B.C.3.D.SO.DO.X.P.R) 
20  CALL  RESETlM.M.L.Kb.XO. B.C. iJ.D.'sO.DO. I.J. XI 
30  CALL  XCK(J1,M.L.T2.I.J.X) 

IF(l-:}:?.EQ.l)  GO  TO  50 
‘  r^=l 

GO  10(50.40). _I1^  _ 

40  iF'Jl.E'j.O)  00  TO  40 

'(■)  :f;4;,.lt.ji)  on  to  20 

■■  I-O-i.EO.ol)  00  TO  ') 


LnLl.  mIL.  17.  uT.  Ill  .u'i.LO.  i  2.  u. 
iF  .IT.OT.O)  GO)  TO  SO 
’5 

oO  CALL  C£LTA(  I7i  J7»L>L1 1 1 >  Ji 
CALL  RCW(K9..A1.L1.T2,X.Y) 

IF(K9.GT.O)  GO  TO  90 

RETIJRN 

90  IFdC.GE.Il?)  Gj  TO  100 

CALL  PI VOT I H. KF , 17 . 07 . LI . A1 . XO, B, C. I . J. X . Y ) 

«5  K'2=K2+1 
GO  TO  30 
100  Kl=4 
RETURN 
END 


SUBROUTINE  RESETUI. N. L.K6. XO. B. C. 3.D. SO, BO. I . J. X ) 

C  RESET  S.D.I.U  AND  X 
LOGICAL  e 
REAL*8  X.X0.S7 

DIMENSION  S(  1). C(l). S(l). Dili, soil). DO(l).  1(1). JlD.Xd) 

CALL  RESDIM.N.S.D.SO.DOI 

K.i.=l 

X0=0.0 

):=l 

CO  20  I1=1.M 
CO  20  ^;i=i,N 
i\7=M*(Jl-l)+Il 
IF(B(K7))  GO  TO  10 
100=11 
■J(K)=v.ll 
K=K+1 

10  B(K7)=.TRUE. 

20  CONTINUE 
25  L1=0 

CO  70  K=1.L 
I7=Ii);:) 

J7=J(K) 

1.7=M*(.J7-l)+I7 
:F(.N0T.B(L7)i  m  TO  70 
S7=S(I7l 
DO  30  INd.L 

:f(I(ihi.ne.i7.cr.:h.eq.):)  so  to  so 

N>l=M*(JilM)-l)+i7 
IFIBOIH))  GO  TO  40 
37=S7-X(IH> 

30  CONTINUE 
00  TO  40 
40  37=D(J7) 

DO  50  IH=1.L 

IFsJdHl.NE.JT.CR.IH.EQ.K)  <X)  TO  50 
),H=M*(.J7-1)+I(IH) 

IF(B(|)}I)1  GO  TO  45 
S7=S7-X(IH) 

50  CONTINUE 

tf)  B(K7)=.N0T.B(K7) 

X(K)=S7 
CO  TO  70 
45  Ll=l 
70  CCNTIiXI.E 

IFdl.EO.  1 )  C4]  TO  25 
DO  SO  i';=l.L 
K7=nKU(K)-D+I(K) 

SO  :<o=xo+C(k:7)«x(k) 
return 
END 
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:3jBR0UTINE  XCK( J1 . tl.L. 72. 1 , J,  X i 
C  RESET  <  m  CHECK  FOR  IMFEASIBILITY 
REALMS  X 

Dit'.ENsicix  iaj,jai,x(U 

.ji=0 

DO  10  K=1,L 
i;7=I(K)+tJ{Ki-l)ifM 
IF  i  DABS  ( X  U'.) ;  .LT.Ti)  Xi“- j=<,i.O 
IFiXiKj.LT.~Ti)  ol=l 
•C  TCNTTM'E 
RET  LRU 
END 


SUBROUTINE  StMP(K,Ll.I.J.X) 
C  SWAP  LOCATIONS  IN  I,J  fWD  X 
REAL»8  X.X9 

DinENSION  I(l)iJ(l).X(l) 
19=1 (LI) 

I(L1)=I(K) 

HK)=I9 

I9=J(L1) 

J(L1)=J(K) 

J(K)=I9 

X9=X(L1) 

X(L1)=X(K) 

X(K)=X9 

RETURN 

END 


SUBROUTINE  RESD(M.N.S,D.SO>00) 
C  RESET  S  k  D  FROM  SO  k  DO 

DlfENSION  SM.),D(1),30(1),D0(1I 
DO  10  1='-*^ 

10  S(I)=SO' 

DO  20 

20  [)(J)=DO(J) 

RETURN 

END 


SUBROUTINE  RCIW(K9,B1,L1,T2>X,Y) 

C  ROW  SELECTION 
REAL *8  X.Al 
DINENSION  X(ll,Y(l) 

K9=0 

A1=1.0E20 
00  10  K=l.Ll 

IF(X(K).LT.-T2.0R.Y(KI.LT.0.0)  GO  TO  10 

IF(X(K).GE.A1)  GO  TO  10 

K9=K 

A1=X(K) 

10  CONTINUE 
DO  20  K=l,Ll 

IF(X(KI.LT.-T2.AND.fiK,).r.E,-Al.AND.Y(K!.LT.0.0)  C>0  TO  SO 
20  rjDNTINUE 
RETURN 
30  K')=X 
A1=-X(KT 
RETURN 
END 
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oUBROUT INE  3ET ; I1irJ)LiP.Oir'.C'i 7_iCi  1 1 oi Ai r iR , 
C  GET  PRICES 

REAL^y  !(i  Pi  RiRUi  i.'.i 

uIMENSION  KlhJdl.iKDiPdl.RlD.Xd) 

"0  IC  I1=1,I1 
10  ?dl;=0.0 

:b  ;c  ji=i,n 

CO  -dl 


:•]  50  f=\,i 

ic=i;r;) 

C0=C(K7) 

IF(K6.EQ.O)  as  TO  30 
!F(X':K).GE.-T2)  C0=0.0 
iF(X(K).LT.-T2)  CO=-RO+I.O 
30  IF(PdO).GT.O.O.AND.R(.JO).GT.O.O)  00  TO  50 
IF(P(I0).EQ.0.C.A.ND.R(J0).EQ.0.0)  SO  TO  40 
IF(P(IO).EQ.O.O)  P(I0)=R(J0l-C0 
IF(R(JO).EC!.0.0)  R(JO)=P(IO)+CO 
GO  TO  50 
40  Ll=l 
50  CONTINUE 

IFiLl.GT.O)  GO  TO  CS 

RETURN 

END 


SUBROUTINE  DELTAd?, J7.L,Lli I, J,X,Y) 

C  COh.PiJTE  PER  UNIT  CHANGES  IN  BASIS  VARIAMS 
R£Al.*3  X 

DIMENSION  Id),Jd)-Xd).Yd) 

10=17 

JCI=U7 

L0=0 

Ll=l 

L2=L 

17=0 

5  DO  10  K=L1,L2 

IF(K7.EQ.O.AND.I(K).EQ.IO)  TO  15 
IF(K7.Eil.l.AND.JIK).EQ.J0)  OH  TO  15 
10  CONTINUE 

IF(LO.EQ.O)  GO  TO  20 

Ll=l.O 

K=L2 

lO=lO-1 

L2=L2-1 

CALL  SHAP(K,L1.I..J.X) 

GO  TO  20 

15  IF(K.GT.Ll)  CALL  SWAP(K.L1.I.J,X) 
IFOd.EC'.O.AND.JlLl !.E'3.J7!  GO  TO  40 
IFiK7.EQ.l.AND.IiLl).EQ.I7)  GO  TO  40 
LO=LOd 

l:=li+i 

20  >7=1-K7 

IF'LC.EG.O)  GO  TO  30 
10=1 iLO) 

JO=U(LO) 

GO  TO  5 
30  10=17 
jO=J7 
GO  TO  5 
40  DO  50  k:=i,l 
50  Y(X)=0.0 
37=1.0 
00  faO  K=1,L1 
Y(K:)=S7 
60  S7=-S7 
RETURN 
END 


Ri;2  THE  IDA/BPT  CRISIS  RELOCATION  PLANNING  MODEL: 

AD  A126  652  DOCUMENTATION  A..(UI  INSTITUTE  EOR  DEFENSE 

ANALYSES  ALEXANDRIA  VA  PROGRAM  ANALYSIS. . 

UNCLASSIFIED  E  S  PEARSALL  ET  AL.  22  DEC  82  F/G  9/2 


NL 
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SUSROUrM  HIN(I7.J7,Dl,D3.K6,T3.f1.N.B.C.P,R) 
C  illN  :t-J,  SELECTS  VARIABLE  TO  ENTER  BASIS 
LOGICAL  B 
REAL»8  P,R,I!1,D2 
DIMENSION  B(l)-C(l),?(l),R(l) 

17=0 

01=03 

DO  DO  1=1  mM 
CO  10  0=1 iN 
K=MKJ-1)+I 

IF(.NOT.BiK))  GO  TO  10 
D2=C(iO 

IF(K6.EQ.l)  D2=0.0 
D2=D2+P(I)-R(.J) 

IF(DABS(D2).LT.-T3)  D2=0.0 

IF(D2.GE.D1)  GO  TO  10 

D1=D2 

17=1 

J7=J 

10  CONTINUE 
20  CONTINUE 
RETURN 
END 


SUBROUTINE  N£W(H,N,L,R0,T2.T3.B.C,S.D,S0.D0.X,P,R) 
C  START  PHASE  1 
LOGICAL  B 
REAL*8  X,P,R,R0,D1 

DIMENSION  B(1).3<1).D(1),S0(11,D0(1),C(1) 

DIMENSION  Xdl.Pdl.Rd) 

CALL  RESDdl.N.S.D.SO.DO) 

U)  5  10=1, M 
DO  5  JO=I,N 

t:=M»(Jo-i)+io 

B(K)=.TRUE.  • 

5  CONTINUE 
DO  10  10=1, H 
10  P(10)=-R0 
DO  20  J0=1,N 
20  R(JO)=RO 
DO  30  10=1, H 
DO  30  J0=1,N 
K=M*(J0-1)+I0 

IF(C(K).LT.-P(I0))  P(I0)=-C(K) 

30  CONTINUE 
DO  40  J0=1,N 
DO  40  10=1, M 
K=M«<J0-1)+I0 

IF(C!K)+PaO).LT.R(.JO))  R(J0)=C(H;)+P(I0) 

40  CONTINUE 
D3=1.0E20 
i;=i 
K6=0 

45  CALL  MIN(I7,J7,D1,D3,K6.T3,I1.N.B.C.P.R) 
i-.;7=M*(J7-l)+I7 

IF(S(I7).LT.D(J7)d2)  P(I7)=R0 
IF(S(I7).GT.D(J7)-T2I  R(J7)=-R0 
B(K7)=. FALSE. 

IF(Sa7).GT.D(J7))  ijO  TO  50 
D(J7)=DIU7)-S(I7) 

S(I7)=0.0 
GO  TO  60 

50  S(I7)=S(I7)-D(J7) 

0(J7)=0.0 
60  K=K+1 

IF(K.LE.L)  GO  TO  45 

RETURN 

END 
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SUBROUTINE  PIV0T(I1,1':9,I7.  J7,L1.A1.X0,3.C.I. J.X.V) 
PIVOT 

LOGICflU  B 
REAU3  X,XO,Al 

DIMENSION  B(l).C(l).I(ll.U(l).X(l).Y(l) 

J9=0(K9) 

K=tH(U9-l)+I9 

;50=xo-C';k)»x(K9) 

s(k;=.tf.ue. 

K7=l1»(U7-l)ti7 

B(K7)=.FAL3E. 

I(K9)=!7 

J(K9)=J7 

Y(K9)=-1.0 

X(K9)=0.0 

DO  10  K=1,L1 

K7=tt»(J(K)-l)+IlK) 

X0=X0-ftl»V(K)»CIK71 

X(K)=X(X.)-ftUY(K) 

10  CONTINUE 
RETURN 
ENU 


I 


4.5  Data  Files 

MAP  -  Map  File 
Field  Format 


2  16 

3  16 

4  16 

5  16 

6  2X,4A4 

USND  -  Nodes  File 
Field  Format 

1  12 

2  13 

3  4X,I1 

4  11X,11,4A4 

5  8X,2F9.4 

6  3A4 

7  2X,F8.0 

8  F8.0 

9  4X,I2 

USNL  -  Links  File 
Field  Format 

1  12 

2  13 

3  IX, II 

4  215 

5  IX, 12 

6  13 

7  IX, 11 

8  215 

9  IX, 12 

10  15 


Contents 

State  FIPS  number 
USND  starting  record  number 
USND  ending  record  number 
USLD  starting  record  number 
USLD  ending  record  number 
State  name 


Contents 

State  FIPS  number 
County  FIPS  number 
Sequence  number 

Urban/Rural  flag.  County  center  name 
Latitude,  Longitude  in  degrees 
County  name 
1980  Population  (or  0) 

County  area  (or  0) 

FQIA  risk  code 


Contents 

Origin  state  FIPS  code 
Origin  county  FIPS  code 
Origin  sequence  number 
Origin  northing,  easting 
Destination  state  FIPS  code 
Destination  county  FIPS  code 
Destination  sequence  number 
Destination  northing,  easting 
Route  characteristics 
Highway  number 
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MAP 


1 

1 

72 

i 

200 

0 

.1 

f\ 

■i 

0 

0 

0 

0 

73 

li'X> 

201 

101 

m 

251 

412 

i 

177 

413 

613 

7 

0 

0 

0 

fi 

3 

287 

359 

619 

770 

9 

360 

396 

771 

370 

10 

397 

404 

371 

394 

11 

405 

405 

395 

907 

12 

406 

484 

908 

1116 

13 

485 

646 

1117 

1463 

14 

0 

0 

0 

0 

15 

0 

0 

0 

0 

16 

647 

695 

1469 

1555 

17 

696 

808 

1556 

1899 

18 

309 

907 

1900 

2177 

19 

903 

1025 

2173 

2452 

20 

1026 

1142 

2453 

2698 

21 

1143 

1265 

2699 

2988 

22 

1266 

1331 

2939 

3123 

23 

1332 

i;556 

3129 

3185 

24 

1357 

im 

3186 

3273 

25 

1391 

1429 

3274 

3385 

2b 

1430 

1524 

3386 

3572 

27 

1525 

1613 

3573 

3730 

28 

1619 

1702 

3781 

3954 

29 

1703 

1820 

39^ 

4213 

30 

1821 

1682 

4214 

4319 

;u 

1333 

1973 

4320 

4525 

32 

1979 

2006 

4526 

4581 

33 

2007 

2039 

4582 

4673 

34 

2040 

2104 

4674 

4364 

35 

2105 

2144 

4865 

4950 

36 

2145 

2282 

4951 

5306 

37 

2283 

2383 

5307 

5545 

38 

2384 

2437 

5546 

5638 

39 

2438 

2556 

5639 

5982 

40 

2557 

2636 

5983 

6184 

41 

2637 

2687 

6185 

6282 

42 

2688 

2784 

6233 

6550 

43 

0 

0 

0 

0 

44 

2735 

2799 

6551 

6588 

45 

2300 

2846 

6589 

6750 

4b 

2847 

2^n;5 

6751 

6381 

47 

2916 

3012 

6882 

7131 

4:3 

3013 

3297 

7132 

7751 

49 

3298 

:333l 

7752 

7314 

50 

3332 

3352 

7315 

7373 

51 

3353 

3497 

7374 

3200 

52 

0 

0 

0 

0 

53 

3498 

3549 

8201 

3292 

54 

3550 

360? 

3293 

342? 

55 

3610 

3637 

3430 

3606 

5b 

5688 

3713 

8607 

3675 

57 

0 

0 

0 

0 

53 

0 

0 

0 

0 

59 

0 

0 

0 

0 

60 

0 

0 

0 

0 

ARIZONA 

ARKANSAS 

CALIFORNIA 

CCiLORADO 

CQNtCCTlCUT 

[CLAWARE 

DIST  OF  COLUMBIA 

FLORIDA 

C£CRG1A 


IDM! 

ILLINOIS 

INDIANA 

I0«A 

KANSAS 

KENTUCKY 

LCUISIANA 

MAINE 

MARYLM) 

MASSACHUSETTS 

MICHIGW 

MINNESOTA 

MISSISSIPPI 

MIS:50URI 

MONTANA 

NEBRASKA 

NEVADA 

NEW  HAI5>SHIRE 
NEW  vERSEY 
NEW  (€XIC0 
NEW  YORK 
»3RTH  CAROLINA 
^aRTH  DAKOTA 
CWIO 
(KIWIDMA 
OFECm 

PENNSYLVANIA 

RHODE  ISLAND 

SOUTH  CARCLINA 

SOUTH  DWQTA 

TENNESEE 

TEXAS 

UTAH 

VERMCWT 

VIRGINIA 

WAa<INGTCN 
WEST  VIRGINIA 
WISCONSIN 
WYOMING 


USND  -  Partial  Dump 


24001 

0 

114 

66  0Cl'J''BERLANL 

73.7556Ai.LEGHNY 

30548. 

“26. 

•■A 

2^003 

0 

A5 

187  C'A!vi\AFuLIS 

33.9477 

76.5142ANNE  ARUNDEL 

:-70775, 

417. 

0 

103 

130  0CCCfEV3VILLE 

3  • .  -riC"? 

7i.6;':LEA.':MCF.E 

.  -1^/3  • 

;  1 . 

^  .• 

24009 

0 

33 

134  oPRInCE  FREDRICK 

33.5500 

76. 5600C  AL'.’ERT 

"19. 

24011 

0 

58 

224  ODENTON 

38.3500 

75.8300CAR0LINE 

23143. 

320. 

•A.'*, 

24013 

0 

lilS 

161  CliESTHINISTER 

39.568? 

77.C000CARRGLL 

?6-j56. 

45  i. 

.)() 

24015 

0 

112 

226  OELKTON 

39.62(10 

75.7326CECIL 

60430. 

352. 

CO 

24017 

0 

43 

165  OWALDORF 

38.6312 

76.9156CHARLES 

72751. 

458. 

40 

24019 

0 

40 

212  0CN1BRI2GE 

38.5834 

76.0553D0RCHESTER 

30623. 

580. 

CO 

24021 

0 

97 

137  OFREDRICK 

39.4122 

77.4472FREDERia< 

114263. 

664, 

:>o 

24023 

0 

118 

44  OGRANTSVILLE 

39.6334 

79.169SGARREn 

2649E:. 

.663. 

00 

24025 

0 

107 

197  OBEL  AIR 

39.5600 

76.3200HARFOTD 

145930. 

i48. 

00 

24027 

0 

39 

169  OCOLUMBIA 

39.2”?23 

76.S514H0WARD 

116572. 

251. 

00 

24029 

0 

82 

211  OCHESTERTON 

39.1922 

76.0600KENT 

16695. 

2S4. 

00 

24031 

0 

31 

145  OGERMANTOWN 

39.18:34 

77.2902110NTG0MERY 

579053. 

494. 

20 

24033 

0 

54 

177  OUPPER  MARLBORO 

38.7811 

76.7C12PRINCE  GEORGES  665071. 

485. 

20 

24035 

0 

63 

210  OCENTERVILLE 

33.9922 

76.0900QUEEN  ANNES 

.55^)8. 

373, 

CO 

24037 

0 

17 

190  OLEXINGTON  PARK 

38.2466 

76.4584ST  MARYS 

59895. 

567, 

00 

2403? 

0 

13 

231  OPRIfCESS  ANNE 

33. 1967 

75.6956SGMERSET 

19133. 

332. 

no 

24041 

0 

52 

211  OEASTON 

33.7600 

76.0742TALB0T 

15604. 

279. 

00 

24043 

0 

113 

120  OHAGERSTOHN 

39.6389 

77.7586UASHINGT0N 

113086. 

462. 

30 

24045 

0 

25 

236  03ALISBURY 

•38.3700 

75.6000WIDjMIC0 

64540. 

3£0. 

CO 

24047 

0 

3 

239  OPOCOMOKE  CITY 

38.0521 

75.5431W0RCESTER 

30889. 

483. 

00 

24510 

0 

90 

181  IBALTIMORE 

:39.3100 

76. 621 IBALTIMORE. MD  P  786775. 

7<j 

K) 

24003 

1 

61 

174  ODAVIDSONVILLE 

38.  SS-?? 

76.7559ANNE  ARUNDEL 

0. 

0. 

00 

24015 

1 

103 

215  OPERRYVILLE 

39.5622 

75.9872CECIL 

0. 

0. 

00 

24021 

1 

117 

142  OEMIITSOWG 

39.7012 

77.3542FREDER1CK 

0. 

0. 

00 

24023 

1 

38 

28  ORBDHOUSE 

39.2738 

79.4674CARREn 

0. 

0, 

00 

24025 

1 

105 

207  OABERDEEN 

39.5311 

76.1341HARF0RD 

0. 

0. 

00 

24035 

1 

77 

217  OY 

39.1223 

75.959aiUEEN  ANNES 

0. 

0. 

00 

24043 

1 

113 

98  OHANCOCK 

39.7112 

78.1677MASHINGT0N 

0. 

0. 

00 

24047 

1 

23 

262  OOCEWl  CITY 

38.3412 

TS.llSSWCffCESTER 

0. 

0. 

00 

24015 

2 

115 

208  OLIBERTY  GROVE 

39.6634 

76.1174CECIL 

0. 

0. 

00 

24021 

2 

93 

130  OJEFFERSQN 

39.3544 

77.5774FREDERICK 

0. 

0. 

00 

25001 

0 

51 

178  OSAG/mS 

41.7389 

70.5233B<«NSTABLE  . 

147925. 

399. 

00 

25003 

0 

99 

39  IPITTSFIELD 

42.4434 

73.2300BERK3HIRE 

145110. 

942. 

40 

25005 

0 

61 

149  OTAINTON 

41.8878 

71.1017BRI3TOL 

474641. 

556. 

40 

25007 

0 

29 

100  oedgwtohn 

41.4200 

70.4933DUKES 

8942. 

106. 

00 

25009 

0 

106 

15^.  0IWNVER;5 

42.5455 

70.9994ESSEX 

633632. 

500. 

40 

25011 

0 

109 

73  OGREENFIELD 

42.5800 

72.5700FRANKLIN 

64317. 

7-22. 

00 

25013 

0 

76 

76  ISPRINGFIELD 

42.1100 

72.5200HAMPDEN 

443013. 

621. 

31 

25015 

0 

94 

74  ONORTHHAfPTQN 

42.3678 

72.5417HAMPSHIRE 

138313. 

537. 

00 

25017 

0 

107 

127  OLITTLETON 

42.5523 

71,5172MIDDLESEX 

1367034. 

829. 

40 

25019 

0 

17 

201  ONANTUCKEH 

41.2600 

70.0900NANTUCKET 

5087. 

46. 

40 

25021 

0 

75 

l;32  OFRAMCLIN 

42.0911 

71.4261N0RF0LK 

606587. 

398, 

00 

25023 

0 

65 

172  OPLYMCtJTH 

41.9422 

70.6450PLYnOUTH 

405437. 

664. 

40 

25025 

0 

91 

150  IB3ST0N 

42.3292 

71.0300SLFFOU<  PSD 

650142. 

55. 

20 

25027 

0 

87 

114  OUORCHESTER 

42.2666 

71.7806U0RCE5TER 

646352. 

1532. 

40 

25003 

1 

37 

36  OLEE 

42.2699 

73.2S33BERK3HIRE 

(1. 

0. 

no 

25005 

1 

48 

146  im  BEDFORD 

41.6999 

71.1t00BRIST0L 

6. 

0. 

00 

USNL  -  Partial  Dump 


24031  0 

31 

145  24021  0 

97 

137  21 

701 

24021  0 

01 

137  24043  0 

113 

120  21 

70 

24043  0 

113 

120  24043  1 

113 

93  21 

70 

24510  0 

90 

131  24027  0 

16V  il 

702 

24027  0 

35 

169  24021  0 

97 

137  21 

702 

24510  0 

90 

181  24005  0 

103 

180  21 

81 

24021  0 

97 

137  24043  0 

113 

120  52 

40 

24043  0 

113 

120  24043  1 

118 

98  52 

40 

24043  1 

118 

98  24001  0 

114 

66  52 

40 

24001  0 

114 

66  24023  0 

116 

44  52 

40 

24001  0 

114 

66  24023  0 

116 

44  32 

48 

24003  1 

61 

174  24003  0 

65 

187  32 

50 

24003  0 

65 

187  240.35  0 

68 

210  32 

50 

24035  0 

68 

210  24041  0 

52 

211  42 

50 

24041  0 

52 

211  24019  0 

40 

212  42 

50 

24019  0 

40 

212  24045  0 

25 

236  42 

50 

24045  0 

25 

236  24047  1 

23 

262  42 

50 

24017  0 

43 

165  24037  0 

17 

190  43 

235 

24033  0 

54 

177  24009  0 

38 

184  43 

4 

24510  0 

90 

181  24003  1 

61 

174  43 

3 

24003  1 

61 

174  24033  0 

54 

177  42 

301 

24033  0 

54 

177  24017  0 

43 

165  42 

301 

24510  0 

90 

181  24003  0 

.65 

137  43 

505 

24021  1 

117 

142  24021  0 

97 

137  52 

15 

24021  0 

97 

137  24021  2 

93 

131)  42 

15 

24021  1 

117 

142  24013  0 

108 

161  53 

97 

24013  0 

108 

161  24510  0 

■90 

ISl  43 

97 

24510  0 

90 

181  24025  0 

107 

197  52 

1 

24025  0 

107 
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The  IDA/BPT  Crlali  Belocetlon  Planning  Kxlal:  Daacrlptton, 

Docunentatlon  and  Uaers'  Guide  to  the  Conputer  Prograa.  (tOa 
Record  Docunent  I>-11)  by  Edward  S.  Paaraall,  IPT,  Robert  C. 

Bushnell,  BPT  under  aubcontract  for  the  Inatitute  for  Dafenac 
Analyaea,  (Contract  FEHA  E1W-C-0749C) 

Abstract 

This  report  describes  work  perforaed  by  Bushnell,  Pearsall  and 
Trozzo,  Inc.,  under  aubcontract  with  the  Institute  f<'r  Defense  Analyses 
on  Task  A-l  of  IDA  Contract  No.  £MU-C*074B  with  the  Federal  Eiaergency 
Managenent  Agency.  Task  A-1  calls  for  the  developaent  of  "a  isodel  to 
simulate  population  aovement  during  an  evacuation  from  the  risk  area 
to  the  various  host  areas  over  a  transportat  Ion  network.*' 

This  report  describes,  documents  and  provides  a  user's  guide  to  a 
system  of  computer  routines  which  perform  the  various  computations  re¬ 
quired  to  apply  a  crisis  relocation  model  developed  Jointly  by  IDA  and 
BPT,  Inc.  The  computer  routines  together  comprise  an  Interactl^  system 
resident  on  the  FCMA  Univsc  1106  facility.  The  sndel  and  its  attached 
national  data  baae  can  be  used  to  analyze  in  detail  the  evacuation  of 
risk  areas  anywhere  in  the  continental  (/nlced  States  under  a  wide  range 
of  different  assumptions  regarding  the  assignment  of  reception  areas  and 
the  performance  of  the  transportation  system  during  the  evacuation. 


The  IDA/BPT  Crisis  Relocation  Planning  Model:  Description, 

Documentation  and  I'scrs*  Guide  to  the  Cosiputer  Program.  (IDA 
Record  Document  >11)  by  Edward  S.  Pearsall,  BPT,  Robert  C. 

Bushnell,  BPT  under  subcontract  for  the  Institute  for  Defense 
Analyses,  (Contract  FEMA  EfW-C-0749C) 

Abstract 

This  report  describes  work  performed  by  Bushnell,  Pearsall  and 
Trozzo,  Inc.,  under  subcontract  with  the  Institute  for  Defense  Analyses 
on  Task  A't  of  IDA  Contract  No.  CMW-C''0749  with  Che  Federal  Emergency 
Management  Agency.  Task  A-l  calls  for  the  development  of  "a  model  to 
simulate  population  movement  during  an  evacuation  from  Che  risk  area 
to  the  various  host  areas  over  a  transportation  network." 

This  report  describes,  documents  and  provides  a  user's  guide  to  a 
system  of  computer  routines  which  perform  the  various  computations  re¬ 
quired  to  apply  a  crisis  relocation  model  developed  Jointly  by  IDA  and 
BPT.  Inc.  The  cosiputer  routines  together  comprise  an  interactive  system 
resident  on  the  FEMA  Univsc  1100  facility.  The  model  and  its  attached 
national  data  base  can  be  used  to  analyte  in  detail  the  evacuation  of 
risk  areas  anywhere  In  the  continental  United  States  under  a  wide  range 
of  different  assumptions  regarding  the  assignsient  of  reception  areas  and 
the  perforsiance  of  the  transportation  system  during  the  evacuation. 
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Bushnell,  BPT  under  subcontract  for  tha  Inatitute  for  Defense 
Analyses,  (Contract  FEMA  EW-C-0749C) 

Abstract 

This  report  describes  work  perforswd  by  Bushnell,  Pearsell  and 
Trozzo,  Inc.,  under  subcontract  with  the  Institute  for  Defense  Analyiq 
on  Task  A^l  of  IDA  Contract  No.  EMW-C-0749  with  the  Federal  Eaergencf 
.Management  Agency.  Task  A-l  calls  for  the  development  of  "a  model  to 
simulate  population  movement  during  an  evacuation  from  the  risk  area 
to  the  various  host  areas  over  a  transportation  network." 

This  report  describes,  documents  and  provides  a  user's  guide  to  • 
system  of  computer  routines  which  perform  the  various  computations  rmH 
qulred  to  apply  a  crisis  relocation  model  developed  Jointly  by  IDA  and 
BPT,  Inc.  The  computer  routines  together  comprise  an  Interactive  syM 
resident  on  the  FEMA  Univsc  U08  facility.  The  model  and  its  attacbM 
national  data  base  can  be  used  to  analyze  in  detail  the  evacuation  of] 
risk  areas  anywhere  In  the  continental  United  States  under  a  wide  raM 
of  different  aseumptlons  regarding  the  assignment  of  reception  areas  I 
the  perforaunce  of  the  transportation  system  during  the  evacuation.  1 
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